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k1 general  matrix  formulation  of  the  aero-servo-elastic  stability  problem 
is  presented  based  on  the  British  method  of  flutter  analysis.  The  analysis 
applies  to  the  closed-loop  system,  and  the  digital  computer  program  (MPASES) 
developed  to  perform  the  analysis  is  described.  An  application  is  made 
to  a typical  air-to-air  missile  with  three  feedback  loops. 
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SUMMARY 


A general  matrix  formulation  of  the  aero-servo-elastic  stability 
problem  for  a closed-loop  system  is  presented.  It  is  based  on  a classi- 
cal method  for  determining  servomechanism  stability  (using  linear 
differential  equations)  and  the  British  method  of  flutter  analysis  which 
represents  the  aerodynamic  forces  as  frequency  dependent  springs  and 
dampers.  The  combination  of  the  aerodynamic  forces  with  the  mechanical 
springs  and  dampers  in  the  equations  of  motion  and  then  with  their 
electro-mechanical  equivalents  in  the  servo  system  leads  to  a consistent 
formulation  of  the  linear  equations  of  motion  for  the  closed-loop  aero- 
servo-elastic  system.  An  iterative  real  eigenvalue  solution  accounts 
for  the  frequency  dependence  of  the  aerodynamic  forces  which  is  only 
secondary. 


The  digital  computer  program  MPASES  (modified  program  for  aero- 
servo-elastic  stability)  developed  to  perform  the  analysis  is  described. 
Dynamic  storage  allocation  is  utilized  throughout  MPASES  to  provide  the 
most  efficient  use  of  core  storage  through  the  variable  dimensioning  of 


all  arrays.  Program  limitations  are  minimal  allowing  stability  results 
in  one  computer  run  to  be  obtained  for  up  to  ten  velocities  at  each  of 
five  altitudes.  In  addition,  as  many  as  twenty-five  servo  element  gain 
variation#  are  permitted.  fTj  ' i-, — r. 
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SECTION  I 


INTRODUCTION 

This  report  documents  a modification  to  Ref.  1 which  presented  a digital 
computer  program  for  the  analysis  of  aero-servo-elastic  system  stability. 
Reference  1 formulated  the  problem  by  combining  one  of  the  classical  methods 
for  determining  servomechanism  system  stability  with  the  American  method  of 
flutter  analysis.  Since  the  American  method  of  flutter  analysis  represents 
the  aerodynamic  forces  as  frequency  dependent  masses  and  utilizes  the  mathe- 
matical concept  of  an  artificial  structural  damping  to  determine  the  reduced 
frequency  at  flutter,  the  combined  aero-servo-elastic  stability  analysis  does 
not  obtain  physically  meaningful  frequencies  and  dampings  at  flight  conditions 
other  than  neutrally  stable  ones.  On  the  other  hand,  the  British  method  of 
flutter  analysis  represents  the  aerodynamic  forces  more  realistically  as 
frequency  dependent  springs  and  dampers.  This  is  only  an  approximation  for 
transient  motion  since  the  aerodynamic  forces  are  generally  known  only  for 
harmonic  motion  and  it  is  not  completely  correct  to  identify  the  forces  in 
phase  with  the  displacements  as  aerodynamic  springs  and  the  forces  in  phase 
with  the  velocities  as  aerodynamic  dampers.  However,  comparisons  of  results 
from  the  various  methods  for  flutter  analysis  have  been  made  by  Jocelyn 
Lawrence  and  Jackson  in  Ref.  2,  and  it  has  been  found  that  the  British 
method  of  flutter  analysis  leads  to  reasonable  predictions  of  transient  aero- 
elastic  behavior  when  the  damping  levels  are  low. 

The  British  method  of  treating  the  aerodynamic  loads  as  springs  and 
dampers  permits  combining  the  aerodynamic  forces  with  the  mechanical  springs 
and  dampers  in  the  equations  of  motion  and  then  with  their  electro-mechanical 
equivalents  in  the  servo  system.  This  results  in  a consistent  formulation  of 
the  equations  of  motion  for  the  closed-loop  aero-servo-elastic  system  because 
the  frequency  dependence  of  the  aerodynamic  forces  is  only  secondary. 


Although  the  frequency  dependence  is  secondary,  however,  it  is  important  and 
an  iteration  is  necessary  to  "line-up"  the  reduced  frequency  for  which  the 
aerodynamics  are  determined  and  the  frequency  determined  by  the  equations 
of  motion;  one  iteration,  however,  is  sufficient  to  achieve  sufficient 
accuracy  between  the  aerodynamic  forces  and  the  equations  of  motion. 

The  iteration  begins  at  a particular  velocity  by  assuming  the  reduced 
frequency  k*0.  (The  aerodynamic  damping  can  only  be  found  from  a low  fre- 
quency k«0.0  but  an  interpolation  scheme  will  estimate  the  damping  at  k«0.) 
Solving  the  eigenvalue  problem  from  the  equations  of  motion  yields  either 
separate  and  real  roots  or  complex  conjugate  pairs.  Any  real  roots  found 
for  k»0  are  correct,  e.g.,  the  rigid  body  roll-damping  root  or  a structural 
divergence  root.  However,  the  complex  conjugate  roots  must  be  lined  up  and 
the  lowest  frequency  roots  permit  the  aerodynamic  forces  to  be  revised  (by 
interpolation)  to  the  corresponding  reduced  frequency.  With  the  revised 
aerodynamics,  the  eigenvalue  problem  is  solved  again  and  the  lowest  oscilla- 
tory root  is  obtained  along  with  an  estimate  of  the  next  higher  frequency. 
The  aerodynamic  forces  are  then  revised  again  to  correspond  to  the  estimate 
of  the  reduced  frequency  of  the  second  oscillatory  mode  and  the  new  eigen- 
value problem  is  solved.  This  results  in  the  second  oscillatory  root  and  an 
estimate  of  the  third  frequency.  The  process  of  revising  the  aerodynamics 
and  solving  new  eigenvalue  problems  continues  until  the  desired  number  of 
roots  have  been  obtained  for  the  particular  velocity.  The  whole  procedure 
is  repeated  for  the  next  higher  velocity  and  is  continued  until  the  velocity 
range  of  interest  has  been  covered.  From  the  solutions  for  frequency  and 
damping  of  each  mode,  root  locus  charts  can  be  drawn  for  use  in  redesign  of 
servo  components,  e.g.,  amplifier  gains,  or  mechanical  components,  e.g., 
mass  balancing. 


2 


The  Modified  jgrogram  for  aero-servo-elastic  stability  (MPASES)  requires 
four  revisions  to  the  program  PASES  of  Ref.  1.  The  first  is  the  inclusion 
of  the  aerodyne*,  forces  as  spring  and  damping  terns  in  the  equations  of 
notion  rather  than  as  complex  inertial  terns.  The  second  is  the  iterative 
eigenvalue  solution  required  to  line  up  the  frequencies  between  the  aero- 
dynamic forces  and  the  equations  of  motion  as  discussed  above.  The  third 
problem  is  the  interpolation  of  the  aerodynamic  forces  necessary  in  the 
iteration  to  minimize  aerodynamic  computational  expense.  The  last  problem 
is  the  result  of  the  new  formulation  dealing  only  with  real  matrices;  a 
specialized  eigenvalue  extraction  method  that  analyzes  real  matrices  which 
have  either  real  or  complex  conjugate  roots  may  be  utilized  for  computational 
efficiency.  Each  of  these  modifications  is  discussed  in  the  following  sections. 


SECTION  II 

CLOSED-LOOP  AERO-SERVO-ELASTIC  STABILITY  ANALYSIS 


New  Equations  of  Motion 

The  theoretical  derivation  of  the  equations  of  motion  for  the  PASES 
computer  program  is  taken  from  Ref.  1 and  is  reproduced  in  this  report  in 
the  Appendix  for  ease  of  reference.  It  is  only  necessary  here  to  rederive 
the  aeroelastic  equations  of  motion  since  the  servo  equations  are  not 
affected  by  the  change  in  representation  of  the  aerodynamic  forces  from 
complex  masses  to  real  springs  and  dampers. 

Our  new  definitions  of  aerodynamic  influence  coefficients  (AIC's)  are 
taken  from  the  survey  of  unsteady  AIC's  in  Ref.  3.  The  unsteady  force  is 
given  in  terms  of  the  deflections  and  their  velocities  by 


<F}  - (qS/c)  ([Cjl5]{h}  ♦ [C^Khc/V}) 


The  AIC's  [Cj^]  and  [C^]  ®ay  be  regarded  as  static  and  dynamic  stability 
derivatives  but  are,  in  general,  functions  of  the  reduced  frequency 


k » ub  /V 
r 


where 


Assuming  the  AIC's  to  be  constants  independent  of  k simplifies  the  analysis 
and  is  adequate  for  low  frequencies,  but  it  is  inaccurate  at  higher  fre- 
quencies and  the  reduced  frequencies  should  be  lined  up  with  the  frequencies 
obtained  from  the  equations  of  motion  as  discussed  in  the  Introduction. 
However,  the  choice  of  lining  up  the  frequencies  or  not  is  left  to  the 
user  in  his  choice  of  frequency  dependence- of  the  AICs. 
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A definition  of  complex  oscillatory  AIC's  was  also  given  in  Ref.  3 as 


{F}  - Pa^bJsCCbKh}  (4) 

for  use  in  the  American  method  of  flutter  analysis.  Equations  (1)  and  (4) 
must  be  identical  for  harmonic  motion  and  we  have 

(qs/c)([chs]  ♦ iCttc/vjtCn^th)  » Pa)2brs(cchR]  + 1Cchi]) Ch}  (S) 

where  [C^R]  and 
Identifying  the  real  and  imaginary  parts  of  Eq.  (5)  leads  to 

[C^]  - 2k2(£s/S)[ChR]  (6) 

and 

t°hDh]  * k Ccs/S)  [ChI]  (7) 

Equations  (6)  and  (7)  determine  the  aerodynamic  stiffness  and  damping  AIC's 
when  the  oscillatory  AIC's  are  given  in  the  format  of  Eq.  (4),  as  is  the 
case  with  a number  of  available  computer  programs,  e.g..  Refs.  4 and  5. 

In  terms  of  the  stiffness  and  damping  AIC's,  the  equations  of  motion  of 
the  aeroelastic  system  including  the  rotations  of  control  surfaces  appear  as 
(cf.  Appendix  A,  Eq.  (2-15)) 

[M] {h}  ♦ (C]{h)  + [K]({h>  - (h6]{6»  - {F}  (8a) 

- CqS/c)([Ch3]{h}  ♦ [^Hhc/V})  (8b) 

[M]{h>  ♦ (C](h>  ♦ [R]{h}  - [K][h5]{6>  - 0 (9) 


[C^j]  are  the  real  and  imaginary  parts  of  [C^] , respectively. 


or 


where 


[C]  - [c]  - Jjpvstc^] 

(10) 

[R]  - [K]  - (qS/c)  [C^SJ 

(11) 

The  bo del  solution  proceeds  as  before  (see  Appendix  A) . The  series  for 
the  deflections  is 

{h>  - [hpjfap}  * ♦ [h5](6}  (12) 

T 

Substituting  Eq.  (12)  into  Eq.  (9)  and  premultiplying  by  [hp)  leads  to 
the  modal  equations  for  the  flexible  degrees  of  freedom. 

c*Ficy  ♦ tw^ioy  + c*FJ]«> 

♦ (CpHij)  * (fFRKljl  + [fF4)tS>  (13) 

♦ [*FK«p>  ♦ [Sf,]!*,)  + t*p4K«>  * o 

where 


[Mp]  * 

[hF]T(M][hF] 

(14) 

[«pr)  ■ 

(hF]T[M] (hR] 

(15) 

t*p4]  ■ 

[hp3T(M] [h6] 

(16) 

(Spl  - 

(hpftCUhyl 

(17) 

[hF3T[CHhR] 

(18a) 

■ 

-l»VSthpj't[ChDh]  [hR] 

(18b) 

■ -*vsthFlXDhithsl 

(19b) 

[*f]  • [Hp]T[R][hF] 

(20) 

(iFR]  - [hF]T[R][hR] 

(21a) 

. -US/c)thF]T[ChJ][hR] 

(21b) 

t*F4]  ■ li.F]Tt*H»4] 

(22a) 

• -CqS/£)  [hp^tC^j]  [ha] 

(22b) 

Equation  (13)  may  be  compared  to  Eq.  (2-20)  in  Appendix  A.  We  h ave  utilized 
the  fact  above  that  the  rigid  body  displacements  cause  no  internal  damping  or 
structural  forces.  We  have  assumed  the  vibration  modes  are  either  free-free 
modes  for  the  entire  system  or  restrained  modes  for  individual  components, 
but  may  not  be  arbitrarily  chosen  modes.  The  limited  orthogonality  of  the 
free- free  or  restrained  modes  (i.e.,  [Mp]  is  not  a diagonal  matrix  unless 
all  modes  are  free-free  modes)  leads  to  a diagonal  form  for  the  generalized 
stiffness  matrix  TXpJ . If  we  denote  the  diagonal  elements  of  [Mp]  by  , 
then  the  generalized  stiffness  matrix  is 


f*pJ  * [hp]T[K] [hFJ 


(23a) 


■ fwpjfWpj 


(23b) 


and 


[ftp]  • f*Fi  - (qS/c) Chp]T[Chs] [hpl 


(24) 


i dots  the  generalized  stiffness,  but  is  assumed  so  as  an  approximation  that  is 

justified  by  low  levels  of  structural  damping.  The  approximate  form  of  the 

j 

equivalent  viscous  structural  damping  is 

fCpJ  - fgpAupJfKpJ  C2S*) 

• r«Fj  fu*pj  r«Fj  (25b) 

so  that 

[Cf]  - \C?1  - %pVS(hp)T[ChDh][hp)  (26) 

T 

Next,  substituting  Eq.  (12)  into  Eq.  (9)  and  premultiplying  by  [hR] 
leads  to  the  modal  equations  for  the  rigid  body  degrees  of  freedom. 

tiVRF1{aF}  * ♦ ttfR6H6} 

♦ (O^lOp)  ♦ (27) 

* f*RpK*p}  * ■ 0 

where 


[Hyp]  m [«pRl 

(28) 

(Mr]  - (hR]T[M][hR) 

(29) 

- (hR]T(M](h6] 

(30) 

- CHR]T(CHhF] 

(31a) 

• - 4FpWSClkR]TCCh0h]  (hp) 

(31b) 
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(32a) 


tPR]  » thR]TtC][hR] 

» - %pVS[hR3T[ChDh] [hR]  (32b) 

t2R6J  ■ thR]T[C](h5]  (33a) 

- - %pVS[hR]T[ChDh][h6]  (33b) 

t%3  « (Kr3T£R]  thp]  (34a) 

- -CqS/c)(hR]T[ChsHhF]  (34b) 

(Sr  1 - [hR]T[K][hR]  (35a) 

- -(qS/c) [hR]T[Chs] [hR]  (3Sb) 

(*R61  * (36a) 

- -(qS/c) (^R]T[Chsl [hg]  (36b) 


and  we  have  again  noted  that  rigid  body  displacements  produce  no  internal 
damping  or  structural  forces.  Equation  (27)  should  be  compared  to  Eq.  (2-29) 
in  Appendix  A. 

The  new  matrix  partitions  in  the  aero-servo-elastic  equations  of  motion 
now  become 


(xx3 {x} 


>p] 

t*W 

0 

( 

[*rJ 

t*R«l 

0 

0 

0 

rccj 

[CS2] 

] * ( 

(FSAJ 

IRSAJ 

[SC] 

[SS2] 

( ®2 

(37) 


l 

"(tp] 

tfF»l 

0 

K] 

(Xi](i>  • 

ttw] 

tfR) 

0 

111 

K 

A 

0 

0 

fCC) 

[CS2] 

) \ 

JFSGI 

IRSGj 

[SC] 

(SS2]_ 

\ *2 1 

"(*pl 

(^F61 

0 " 

1* f\ 

tXJ(x)  - 

‘V 

0 

JM 

A 

0 

0 

rccj 

[CS2] 

J 6 f 

0 

0 

[SC] 

[SS2] 

I e2  j 

These  may  be  compared  with  Eqs. 

(2-43) 

. (2-44) 

and  (2-4S)  in 

(38) 


(39) 


The  coefficient  matrices  [X>]  and  [X^]  each  have  five  new  nonzero  partitions 
as  expected  from  moving  the  aerodynamic  terms  from  the  mass  matrix  to  the 
stiffness  and  damping  matrices. 

The  Eigenvalue  Problem 

The  new  representation  of  the  aerodynamic  forces  changes  the  eigenvalue 
problem  only  to  the  extent  that  the  matrices  are  real  now  rather  than  complex. 
The  equation  to  be  solved  is  still 


(Y(A]  ♦ [8]){V}  • 0 

where  the  amplitudes  of  motion,  (V),  are  defined  by 

(v)  ■ (V)exp(Yt) 


(40) 


(41) 


and  instability  occurs  when  the  airspeed  and/or  control  system  gains  are 
such  that  the  real  part  of  y is  positive. 


1 
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Although  [A]  and  (B ] art  now  real  matrices,  they  still  nay  be  singular 
and  obtaining  the  canonical  form  of  the  eigenvalue  problem  by  a shift  in 
eigenvalues  is  still  appropriate.  We  let 

Y - Yc  - 1A  (42) 

where  y is  an  arbitrarily  chosen  real  number,  and  then  the  new  eigenvalue  is 


X - 1/(Y0  " Y) 


and  the  new  eigenvalue  problem  is 


X{V)  - Cy0(a]  ♦ Cb|)“1[ai(v> 


The  shift  value  y0  is  arbitrary  to  the  extant  that  it  must  be  chosen  so  the 
linear  combination  yg[A]  + [B]  Is  nonsingular.  A value  which  scales  [a]  to 
be  the  same  order  of  magnitude  as  [B]  and  of  the  same  sign  Is  recommended. 

The  eigenvalues  of  Eq.  (44)  axe  either  real  or  complex  conjugates.  A 

* 

subroutine  ALLMAT  (Ref.  6)  for  complex  matrices  was  used  in  Ref.  1. 

A more  recent  development  for  the  real  case  of  Eq.  (44)  is  the  sub- 
routine El  GW  given  in  the  International  Mathematical  and  Statistical  Library 
(IMSL,  Ref.  7).  Subroutine  EIGBT  calls  IMSL  routine  EBALAF  to  balance  the 

matrix.  Then  IMSL  routine  EHESSF  reduces  the  balanced  matrix  to  an  upper 

\ 

Hessenberg  form  and  routine  EQRH3F  computes  all  of  the  real  and/or  complex 
conjugate  pairs  of  eigenvalues  of  the  Hessenberg  matrix. 

The  IMSL  Package  is  universally  used  and  is  usually  Incorporated  Into 
the  scientific  libraries  of  major  computer  systems.  When  requested,  the 
eigenvectors  are  found  in  two  iterations  by  the  Inverse  Power  Method  with 
Shifts  (Ref.  8,  pp.  323,  626-628)  in  subroutine  EGNVCT  (Ref.  9).  Subroutine 
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[U  - AIJ{u}  ■ 0 (4S) 

! 


The  subroutine  is  used  in  the  present  development  by  setting 

A ■ 0 (46) 

end 

(U1  ■ y[A]  ♦ [B]  (47) 

in  Eq.  (45) . If  eigenvectors  ere  requested,  ell  reel  eigenvelues  ere  used 

in  Eq.  (47) , but  only  the  coaplex  conjugete  eigenvelues  with  positive 
isuiginery  parts  (positive  frequencies)  ere  used,  since  the  eigenvectors  ere 

elso  coaplex  conjugete  peirs. 

H 

Writing  the  coaplex  eigenvalue  es 

Y ■ u ♦ iw  (48) 

where  p is  the  decay  rate  and  w is  the  daaping  frequency,  we  find  the  cyclic 

frequency  to  be 

i 

I 

f ■ uj/Zir  (49) 

end  the  fraction  of  critical  damping  c to  be 

5 ■ -y//p4  ♦ w*  (50) 

For  comparison  to  structural  damping  levels,  twice  the  deaping  ratio  Z is  I 

I 

a preferable  output  quantity  since 


for  a structurally-damped  single  degree  of  freedom  oscillator.  For  a non* 
oscillatory  root  a different  definition  of  damping  ratio  is  necessary  and 
we  chose  the  tine  to  half  amplitude 

T^-in2/(-u)  (52) 

If  the  motion  is  unstable,  i.e.,  p > 0,  Eq.  (52)  gives  the  (negative)  time 
to  double  amplitude.  The  stability  of  the  oscillatory  roots  can  also  be 
compared  using  Eq.  (52)  and  this  will  be  an  additional  output  quantity. 

Lining  Up  the  Reduced  Frequency 

The  need  for  lining  up  the  reduced  frequency  for  a specified  mode  of 
motion  with  the  frequency  determined  by  the  eigenvalue  problem  for  that  mode 
was  discussed  briefly  in  the  Introduction  and  at  some  length  in  Ref.  2.  The 
necessary  equations  for  the  iteration  are  given  in  this  section. 

It  is  possible  to  begin  the  iteration  with  any  value  of  k.  However,  a 
finite  value  of  k may  be  representative  of  an  oscillatory  mode  but  we  are 
equally  interested  in  static  modes.  Therefore,  we  begin  with  k-0  and  any 
real  roots  will  be  determined,  e.g.,  the  roll-damping  root,  a static  struc- 
tural divergence  root,  and  any  over-damped  roots  from  control  system  compon- 
ents. The  first  oscillatory  root  for  a free  vehicle  may  be  either  the  short 
period  mode  or  the  Dutch  roll  mode  and  the  choice  of  k»0  provides  a good 
estimate  of  that. 

Let  the  complex  conjugate  pairs  of  roots  be  denoted  by 

Y ■ p t iut  (53) 

Trs  wrs  rs  v ' 

where  r denotes  the  oscillatory  mode  number  ordered  by  frequency, 

<“L  < u>2j  < ...),  and  s denotes  the  number  of  the  mode  under  investigation. 


The  reduced  frequency  is 


k • wc/2V  (54) 

and  the  aerodynamic  forces  are  determined  for  k?  which  should  be 

k ■ w c/2V  (55) 

9 33 

However,  <*»ss  is  not  known  at  the  outsets  From  the  initial  solution  with  kQ«0, 
we  estimate  the  nonzero  frequencies  (w^,  w2q,  oj^q,  ...)•  Refining  the  aero* 
dynamics  by  interpolation  for  k1  - w10c/2V  and  repeating  the  eigenvalue  solu* 
tion,  we  find  a new  sat  of  frequencies  (w^,  ^2i  • w3i>  •••)•  The  value  of 
is  taken  as  the  correct  value  for  the  first  mode  and  its  damping 
determines  its  stability.  The  eigenvectors  for  the  first  root  may  then 
be  calculated  if  the  mode  shapes  are  desired. 

The  aerodynamics  are  next  determined  by  interpolation  for  k2  * U21c/2V 
and  the  next  eigenvalue  solution  yields  the  frequencies  (w12,  si)j2»  • ••)• 
The  frequency  b>22  is  now  taken  to  be  correct  and  its  damping  a22  measures 
the  second  oscillatory  mode  stability.  We  continue  with  kj  * ta,32br/V* 
finding  the  aerodynamics  for  kj,  and  then  the  eigenvalues;  w32  and  a32  are 
assumed  to  be  correct.  The  process  continues  until  all  roots  of  interest, 
and  all  eigenvectors,  if  requested,  have  been  found. 

Care  must  be  exercised  in  tracking  the  oscillatory  roots  because,  as 
the  reduced  frequency  is  increased  the  roots  which  were  real  at  k-0  may 
become  complex  at  higher  values  of  k or  roots  which  were  complex  at  k-0 
may  become  real  at  higher  k's.  A suitable  algorithm  is  to  choose  wJ$  as 
the  closest  value  to  u.  _ , , and  then  to  choose  u. . , _ as  the  closest  higher 

3 |Sa 1 3^1 |3 

value  to  u>ss. 


t 
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Linear  Spline  Interpolation  of  Aerodynamic  Terns 


In  the  modal  matrix  equations  of  motion,  twelve  of  the  partitions 
depend  on  the  aerodynamic  influence  coefficients  (AIC's)  which  are  functions 
of  the  reduced  frequency  k.  As  the  reduced  frequency  changes  when  the  roots 
are  tracked  in  the  frequency  lining  up  process,  the  AIC's  will  change  and 
interpolation  on  k is  necessary.  It  is  computationally  simpler  to  inter- 
polate the  partitions  rather  than  the  AIC's.  Two  of  the  partitions  depend 
on  structural  parameters  also  but  these  are  independent  of  k and  can  be 
carried  along  in  the  interpolation.  A linear  spline  is  chosen  for  the 
interpolation  because  it  offers  cubic  accuracy  and  continuity  throughout 
the  regions  of  interpolation  and  extrapolation. 

A linear  spline  is  a mathematical  device  for  interpolating  for  a 
function  y(x)  for  all  points  x,  when  y is  known  for  a discrete  set  of 
points,  y ^ » y(xA).  The  spline  passes  through  all  of  the  known  points. 

The  mfthematical  spline  takes  its  name  from  the  plastic  spline  used  by 

A 

draftsmen  for  drawing  curves  through  specified  points.  If  the  plastic 
spline  may  be  regarded  as  a uniform  beam,  the  linear  spline  representing 
it  mathematically  is  a solution  to  the  uniform  beam  differential  equation. 
The  following  derivation  for  the  linear  spline  is  taken  from  Ref.  10,  App.  E, 
by  R.  L.  Harder. 

Ne  wish  to  determine  the  deflection  curve  of  a continuous  beam  over 
multiple  supports.  Consider  the  fundamental  solution  to  the  deflection 
equation 


EI^« 

dx* 


(56a) 


(56b) 


that  is  syaaetrical  about  a support  at  the  origin  x « 0.  The  solution  is 


y ■ A,  ■ B j x | ♦ Cx*  ♦ D | x | 1 


Continuity  of  slope  requires  B ■ 0,  and  equilibrium  with  the  support  load  P 


requires 


c/2 

,ta  / wdx_?  li®  El 
c-0  -e/2  e+0  -c 


t/i 

f *^dx 

h **' 


(58a) 


d»y 

lia  El  =-*- 
tr*C  ***  -e/2 


(S8b) 


Froa  Eq.  (57)  the  third  derivative  is  found  to  be 


(S8c) 


6D  sgn  x 


and  Eq.  (S8b)  becoaes 


(12EI)  D - P 


Combining  Eqs.  (57)  and  (60)  and  generalizing  for  all  supports  at  x * xi 


leads  to 


y(*)  ■ ♦ ct(x  - xt)*  ♦ qjx  - xi|*) 
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y(x)  ■ x*sgn 


a large  distance  from  the  supports,  the  deflection 
For  large  x ♦ t»,  Eq.  (61)  behaves  like 

x I ♦ x2^(Ci  - SQ^sgn  x)  * 0(x)  (62) 


and  the  required  linear  behavior  necessitates 


l Qi  - 0 (63) 

I Ct  - 0 (64) 

l QiX.  • 0 . (65) 

Equations  (63)  and  (65)  are  equilibrium  equations,  and  Eq.  (64)  permits 
writing 

I (a i ♦ Ct(x  - xi)  2)  » aQ  ♦ a^x  (66) 

in  Eq.  (61),  where  aQ  and  a^  are  new  constants,  so  that  Eq.  (61)  becomes 

y(x)  - aQ  ♦ axx  ♦ l Qi|x  - xjs  (67) 


In  the  case  where  the  spline  is  symmetrical  about  the  origin,  as  is 
the  case  for  the  AICs  [Cj^OO]  and  ICj^OO]  which  are  symmetrical  functions  of 
the  reduced  frequency  k,  we  may  use  the  method  of  images  and  Eq.  (67)  becomes 

y(x)  « aQ  ♦ l Q^lx  - xii*  ♦ |x  ♦ xil*)  (68) 


and  only  Eq.  (63)  needs  to  be  satisfied  in  addition  to  Eq.  (68).  The  matrix 
form  for  Eqs.  (63)  and  (68)  in  terms  of  the  specified  values  of  y(x^)  is 


t ' 
ii 


I 
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where  « | - xj*  ♦ |x^  ♦ xi|* 

The  left  hand  side  of  Eq.  (69)  may  be  written 


o 

0 0 0 ...  0 

*1 

1 0 0 ...  0 

>• 

0 1 0 ...  0 

] 

• 

• 

0 0 1 ...  0 

• • • • 

• 

• • • • 

0 0 0 ...  1 

The  aatrix  form  for  an  interpolated  value  of  y (x^) 


11  hi  *k2  ••• 


where 


hi  ' l*k  * xjl’  * l*k  * xjl‘ 


Combining  Eqs.  (69),  (71)  and  (72)  leads  to  the  desired  interpolation  coeffi- 


1 


cients  1^  defined  by 


where  {y^}  denotes 


and 


rN . 


(74) 
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0 

0 

0 

0 


(75) 


i 

< 


Equation  (74)  provides  the  linear  spline  interpolation  scheme  needed  to 

obtain  the  AICs  in  the  frequency  lining  up  process.  The  inverse  required 

in  Eq.  (75)  is  ill-conditioned  for  a large  number  N of  interpolated  points, 

say  N*50;  however,  a large  number  are  not  needed  in  the  flutter  analysis  and 

N will  be  limited  to  N £ 10  in  the  present  applications.* 

* 

A lower  limit,  say  N s 5 is  adequate  for  most  applications  unless  an 
extremely  wide  frequency  spectrum  must  be  covered. 
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SECTION  III 


AERO-SERVO-ELASTIC  STABILITY  ANALYSIS  PROGRAM:  MPASES 

V i 

Program  Description 

MPASES  is  a general  purpose  digital  computer  program  for  the  analysis  of  | 

the  closed-loop  stability  problem.  This  program,  which  is  a modification  of  ! 

PASES  (Ref.  1),  formulates  the  problem  by  combining  a classical  method  for  j 

•i 

determining  servomechanism  system  stability  with  the  more  realistic  British  | 

method  of  flutter  analysis  in  lieu  of  the  American  method  used  in  PASES.  With 
the  input  of  an  arbitrary  number  of  elastic  degrees  of  freedom,  aerodynamic 
influence  coefficients  and  the  control  system  description,  the  stability 

of  a missile  or  aircraft  configuration  as  an  aero-servo-elastic  system  can 

l j 

be  investigated.  When  no  aerodynamics  are  involved,  the  servo-elastic  system  t 

I f 

| stability  can  be  determined.  j 

Dynamic  storage  allocation  utilized  throughout  MPASES  provides  the 

i | j ; 

most  efficient  use  of  core  storage  through  the  variable  dimensioning  of 
all  arrays.  Because  the  program  is  based  on  the  theory  that  permits  the 
inclusion  of  the  aerodynamic  forces  as  spring  and  damper  terms  in  the 

; i 

equations  of  motion  rather  than  as  complex  inertial  terms,  only  real  !' 

matrices  are  used  in  the  eigenvalue  solution.  This  makes  for  additional  j 

economy  in  computer  usage.  j 


MPASES  uses  the  method  of  revising  the  aerodynamics  by  lining  up  the 
frequencies  and  performing  the  spline  interpolation  of  the  aerodynamic 
terms  necessary  in  the  iterative  eigenvalue  solution.  This  process  is 
repeated  until  the  required  number  of  modes  and  eigenvectors  are  obtained 
for  a particular  velocity.  Moreover,  if  the  analysis  Mach  number  (velocity 
determining  factor)  differs  from  the  AIC  Mach  numbers  by  more  than  a speci- 
fied ’deviation'  value  (input  by  the  submitter),  spline  interpolation  of  the 
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L \ 


aerodynamic  forces  is  performed  also.  In  this  case,  the  program  interpolates 
for  Mach  number  before  the  reduced  frequency  interpolation  for  each  mode  is 
carried  out. 


f 


The  variables  in  the  stability  analysis  of  an  aero-servo-elastic  system 
are  servo-gains,  velocities  (through  Mach  number  and  speed  of  sound  input) 
and  altitude  (air  density  input).  An  option  is  provided  to  vary  the  gain 
of  a single  servo  component  or  control  surface  with  the  density  and  Mach 
number  held  constant.  When  this  option  is  executed,  the  coefficients  of  the 
output  of  the  desired  servo  element  or  control  surface  are  divided  by  the 
gain  factor,  K.  This  method  facilitates  programming,  since  there  may  be 
more  than  one  input. 

Since  dynamic  dimensioning  is  used,  the  program  limitations  are  minimal. 
Stability  results  can  be  obtained  for  up  to  five  (S)  altitudes  and  for  each 
density,  as  many  as  tan  (10)  different  velocities  can  be  analyzed.  In 
addition,  twenty-five  (25)  gain  variationa  are  permitted.  The  number  of  AXC 
matrices  that  can  be  input  la  limited  to  ten  (10)  reduced  frequencies  for  each 
of  five  (5)  Mach  numbers.  Thera  are  no  size  restrictions  for  any  of  the  other 
parameters  existing  in  the  program. 

The  total  number  of  memory  units  (words  or  bytes)  required  to  execute 
MPASES  is  cooq? lately  dependent  upon  input.  The  size  of  the  program  is 
reflected  in  the  length  of  the  blank  common  block  found  in  the  'MAIN' 
section  of  the  program.  The  length  can  readily  be  altered  to  accommodate 
different  analyses.  . Likewise,  dimensions  of  the  arrays  limiting  the  pro- 
gram as  stipulated  above  can  be  changed  to  the  desired  size. 
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Subroutine  Description 

MAIN 

Mein  program  for  MPASES. 

Reeds  end  prints  basic  input  data. 

Allocates  dynamic  storage  for  all  arrays. 

Calling  program  for  subroutines  MASS,  SERV0,  F0RMAB,  AIC,  MPSE 
and  MPAS. 

AIC 

Reads  and  prints  complex  AIC  matrices. 

Separates  AIC's  into  real  and  imaginary  parts. 

Formulates  generalized  aerodynamic  forces  and  stores  on  temporary  j 

File  1.  * 

I 

Calls  subroutine  GEMMAT. 

EGNVCT 

Finds  the  eigenvector  from  a real  matrix  for  which  the  complex  ! 

eigenvalue  la  known. 

Uses  the  Inverse  power  method  with  shifts.  L 

The  eigenvector  obtained  Is  complex  with  the  largest  value  | 

normalized  to  unity.  | 

A singular  matrix  returns  an  eigenvector  of  zeros.  fi 

EIGRF 

Computes  eigenvalues  of  a real,  triangular  matrix. 

Calls  routines  from  the  IMSL  Package  (see  Ref.  7). 

Called  by  subroutine  SOLV. 

F0RMAB 

Generates  the  A and  B matrices  which  constitute  the  eigenvalue  p 

problem.  i 

1 i 

Stores  [A]  and  [B]  on  temporary  File  2. 

GENMAT 

Formulates  the  generalized  mass  and  aerodynamic  force  matrices  { 

utilizing  rigid  body,  control  surface  and  vibration  modes.  ij 

Calls  subroutine  MAIMPL.  j 

Called  by  subroutines  AIC  and  MASS.  j 

' ! 

, i 

'i 

INTERP  Generates  the  constant  portion  of  the  linear  spline  interpolation 
coefficients  matrix. 

Prints  message  if  spline  interpolation  matrix  is  singular. 

Calls  subroutine  INVERS. 

INVERS  Finds  the  inverse  of  a real,  square  matrix. 

Called  by  subroutines  INTERP  and  S0LV. 

MASS  Reads  and  prints  the  weight  matrix  in  pounds.  Converts  weight 
to  mass  units  (slugs) . 

Reads  and  prints  rigid  body  control  surface  modes,  free-free 
vibration  modes  and  restrained  (rigid  body)  modes. 

Stores  modal  data  on  temporary  File  4. 

Reads  and  prints  damping  coefficients  and  vibration  frequencies 
(Hertz).  Converts  frequencies  to  rps  units. 

Formulates  generalized  mass  matrix. 

Generates  generalized  stiffness  and  damping  matrices  if  flexible 
modes  are  input. 

Calls  subroutine  GENMAT. 

MATMPL  Multiplies  real,  two-dimensional  matrices. 

If  requested,  transposes  post-multiplier  matrix  to  perform  as  a 
pre-multiplier. 

Called  by  subroutines  INTERP,  SERV0 , GENMAT  and  S0LV. 


Obtains  stability  results  from  aero-servo-elastic  analysis  for 
altitudes  and  velocities  requested. 

Reads  Mach  numbers  at  each  altitude  for  which  analysis  is  to  be  made. 

Calculates  velocity  (from  Mach  number  and  speed  of  sound)  and 
constants  for  the  aerodynamic  forces. 

Performs  Mach  number  interpolation  if  necessary. 

Reads  File  2 - brings  (A)  and  [B]  into  core. 

Executes  gain  option  if  requested;  stores  A and  B matrices  re- 
flecting gain  factor  on  tempoerary  File  3. 


MPAS  Carries  out  spline  interpolation  of  the  aerodynamic  forces  for 
(Cont.)  reduced  frequency  for  each  mode. 

Generates  new  A and  B matrices  with  changed  aerodynamics  for  the 
particular  mode. 

Calls  subroutine  S0LV  for  eigenvalue  solution  of  each  mode. 

Prints  stability  results. 

Prints  eigenvectors,  if  requested,  for  the  analysis. 

Reads  File  4. 

Computes  and  prints  structural  deflections  at  system  mass  points 
for  each  mode. 

Calls  subroutines  INTERP,  SPLINE  and  S0LV. 

MPSE  Obtains  stability  results  from  servo-elastic  analysis  for  number 
of  modes  requested. 

Reads  File  2 - places  [A]  and  [B]  into  core. 

Executes  gain  option  if  requested. 

Calls  subroutine  S0LV  for  eigenvalue  solution  of  required  modes. 
Prints  stability  results. 

Prints  eigenvectors  for  each  mode. 

Note  that  the  eigenvectors  are  always  computed  and  printed  for 
a servo-elastic  analysis. 

Reads  File  4. 

Computes  and  prints  structural  deflections  at  system  mass  points 
for  each  mode. 

Calls  subroutine  S0LV. 

RESULT  Forms  the  array  where  stability  results  for  each  mode  are  stored. 

Calls  subroutine  EGNVCT  to  compute  the  eigenvectors  (only  upon 
request  for  aero-servo-elastic  analysis,  always  for  servo-elastic 
analysis) . 

Called  by  subroutine  S0LV. 

SERV0  Reads  coefficients  from  servo  differential  equations. 

Prints  control  system  description. 

Computes  input  to  rate  gyros  and  accelerometers  from  differentia- 
tion or  interpolation  matrices. 

Determines  maximum  order  of  servo  elements. 

Prints  pertinent  information  concerning  all  servo  elements. 
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SERV0 

(Cont.) 


Generates  [X2],  [XI],  [XO],  [Yl],  [YO] . [ZO] . 
Calls  subroutine  MATWPL. 


S0LV  Prints  the  A and  B matrices. 

Generates  the  dynamic  matrix  (eigenvalue  problem) . 

Calls  subroutine  EIGRF  for  eigenvalue  solution. 

"Lines -up  the  frequencies"  in  aero-servo-elastic  analysis. 
Determines  real  and  oscillatory  roots  for  each  mode. 

Estimates  new  reduced  frequencies. 

Calls  subroutine  RESULT  to  compute  stability  results  and  eigen- 
vector for  each  mode. 

Calls  subroutines  MAIMPL,  INVERS , EIGRF  (from  IMSL  Package), 
and  RESULT. 

Called  by  subroutines  MPAS  and  MPSE. 


SPLINE  Generates  the  linear  spline  interpolation  matrix. 

Determines  the  Lagrangian  coefficients. 

Used  in  the  frequency  lining-up  process. 

Called  by  subroutine  MPAS. 

Programming  Symbols  - A partial  list  of  F0RTRAN  symbols  used  in  MPASES  is 
presented.  Additional  symbols  are  defined  in  the  Input  Instructions. 

1.  Integer  Variables  !<jj 


ID 

Density  index 

IE 

Mode  counter 

IG 

Gain  index 

IK 

Reduced  frequency  index 

IM 

Mach  number  index 

IM0DE 

Actual  numbeT  of  modes  calculated  for  stability 

analysis 

<1 


25 


: *■  i it  i "VfryiiTinit-iru.  > 


mAtmeemarntt 


J 


IOB  First  real  root  (eigenvalue  with  positive  real 

part  and  zero  imaginary  part) 

IOE  Last  real  root 

ISING  Singular  matrix  control.  If  ISING  * 1,  matrix  is 

non-singular;  return  from  subroutine  INVERS. 

J1  BLANK  COMMON  length  required  to  call  subroutine  MASS 

J2  BLANK  COMMON  length  required  to  call  subroutine  SERV0 

J3  BLANK  COMMON  length  required  to  call  subroutine  F0RMAB 

J4  BLANK  COMMON  length  required  to  call  subroutine  MPSE 

or  subroutine  AIC 

J5  BLANK  CONMON  length  required  to  call  subroutine  MPAS 

JT0T  Largest  values  of  Jl,  J2,  J3  and  J4  for  servo- 

elastic  analysis;  or  largest  value  of  Jl,  J2,  J3, 

J4  and  JS  for  aero-servo-elastic  analysis.  JT0T 
is  the  minimum  BLANK  COMMON  length  required  to 
execute  MPASES. 

K0P  Column  assignment  in  A and/or  B matrices  for  zero 

and  first  order  servos. 

K01  Column  assignment  in  A and/or  B matrices  for  second 

order  servo-velocity. 

K02  Column  assignment  in  A and/or  B matrices  for  second 

order  servo-displacement. 

KXM  Control  for  Mach  lumber  interpolation. 

K1  Order  of  servo  element  for  which  the  gain  varies. 

K2  Column  assignment  in  the  X,  Y and/or  Z arrays  for 

servo  whose  gain  varies. 

L Counter  for  coefficients  from  servo  differential 

equations  including  inputs  to  control  surfaces, 
outputs  from  servos,  and  inputs  to  servos. 

MDF  Number  of  mass  points;  used  as  a variable  dimension 

involving  flexible  modes.  (MDF  * 1 when  NFM  » 0) 

MDR  Number  of  mass  points;  used  as  a variable  dimension 

involving  rigid  body  modes.  (MDR  » 1 when  NRM  ■ 0) 

MFM  Number  of  flexible  modes;  used  as  a variable  dimen- 

sion. (MFM  ■ 1 when  NFM  ■ 0) 
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MRM 


Number  of  rigid  body  modes;  used  as  a variable 
dimension.  (MRM  ■ 1 when  NRM  ■ 0) 


Ml 


M2 


M3 


M31 

NFR 

NIO 

NK1 

NM0DE 

NM1 

NPART 

NPR 

NRO 

NSERV 


NT 


NTC 

NTT 

NVEC 

NO 

N1 

N2 


Number  of  flexible  modes  or  rigid  body  modes,  which- 
ever is  larger. 

Ml  or  the  number  of  control  surfaces,  whichever  is 
larger. 

Number  of  reduced  frequencies  or  Mach  numbers,  which- 
ever is  larger.  (AIC  input) 

M3  ♦ 1. 

Sum  of  the  number  of  flexible  modes  and  rigid  body 
modes. 

Counter  for  real  roots. 

NK  ♦ 1 

Sum  of  the  number  of  flexible  modes,  rigid  body 
modes  and  control  surfaces. 

NM  ♦ 1 

Number  of  partitions  in  AIC  matrix.  (NPART  • NA 

if  NA  > 1) 

* 

Counter  for  number  of  actual  modes  in  the  servo- 
elastic  analysis. 

Control  for  determining  the  number  of  real  roots. 

Total  number  of  coefficients  from  servo  differential 
equations  including  inputs  to  control  surfaces, 
outputs  from  servos,  and  inputs  to  servos;  used  as 
a variable  dimension. 

Sum  of  the  number  of  flexible  modes,  rigid  body 
modes,  control  surfaces  and  all  servo  elements. 

Sum  of  the  number  of  flexible  modes,  rigid  body 
modes,  control  surfaces  and  second  order  servo 
elements . 

Order  of  the  eigenvalue  problem,  2*NTC+N0+N1. 

Control  for  calculation  and  printing  of  eigenvectors. 
Number  of  zero  order  servos 
Number  of  first  order  servos 
Number  of  second  order  servos 
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2.  Real  Variables 


AN0RM 

C0NAI 

C0NAR 

C0NI 

C0NR 

DET 

EK 

EW 

VEL 


Tost  for  normalising  structural  deflections. 

Constant  for  generalised  aerodynamic  forces* 
imaginary  part,  -HpVS. 

Constant  for  generalised  aerodynamic  forces* 
real  part,  -HpV^S/c. 

Constant  for  imaginary  part  of  AIC  matrix,  k(cs/S) . 

Constant  for  real  part  of  AIC  matrix,  2kl(cs/S). 

Value  of  determinant  returned  from  subroutine  INVERS 
if  matrix  is  non*singular. 

Reduced  frequency,  bru ✓V,  used  in  the  iterative  pro- 
cedure; determined  by  the  estimated  frequency. w . 

Estimated  frequency  for  successive  iterations  of  the 
eigenvalue  solution;  determined  by  the  'lining-up 
the  frequencies'  process. 

Velocity  for  the  stability  analysis,  ft/sec. 


VELK  Velocity  for  the  stability  analysis,  knots. 


XN0RM  Normalising  factor  for  structural  deflections. 


3.  Integer  Arrays  (Variable  dimensions  indicated  in  parentheses) 
I ANA (NTT) 


INDEX (NTT, 3) 
K0(NSE, 2) 

LI (NTT) 

U(NTT) 

NL(NSERV,2) 


N#R 

(NSE.NSERV) 


Used  in  subroutine  SOLV  as  an  argument  when  subroutine 
RESULT  Is  called.  This  vector  is  eventually  used  in 
EGNVCT  as  Ll  described  below. 

Used  in  subroutine  INVERS  for  working  storage. 

Stores  information  concerning  the  maxim um  order  of 
all  servo  elements  and  their  column  assignments  in 
the  X,  Y and  Z arrays. 

Used  in  subroutine  EGNVCT  to  restore  order  of  the 
elements  in  the  eigenvector. 

Used  in  subroutine  EGNVCT  as  working  storage. 

Stores  information  concerning  the  number  of  the 
servo  element  from  which  there  is  input  and  output. 

Stores  information  concerning  the  order  of  the 
servo  element  coefficients. 
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4.  Real  Arrays  (Variable  dimensions  indicated  in  parentheses) 


A (NTT, NTT) 

API 

(NFR,NM0DE,NK) 

APR 

(NFR,NM0DE,NK) 

AMI 

(NFR,NM0DE,NM) 

AMR 

(NFR,NM0OE  ,NM) 
B (NTT, NTT) 


Matrix  in  the  eigenvalue  problem  formulation. 

Generalized  aerodynamics  used  for  reduced  frequency 
spline  interpolation  - imaginary  part. 

Generalized  aerodynamics  used  for  reduced  frequency 
spline  interpolation  - real  part. 

Generalized  aerodynamics  used  for  Mach  number  spline 
interpolation  - imaginary  part. 

Generalized  aerodynamics  used  for  Mach  number  spline 
interpolation  - real  part. 

Matrix  in  the  eigenvalue  problem  formulation. 


BINV(NTT.l)  Column  of  constants  used  in  subroutine  INVERS. 


BC (NTT, NTT) 
C (NTT, NTT) 
CF(MFM) 


CHDH(NDF,NDF) 


CHS(NDF.NDF) 


DFM(1,MFM) 


DRM(1  ,MRM) 


FRC(M2,M2) 


GAI  (NFR.NM0DE) 


GAR(NFR,NM0DE) 


Yq[A]  ♦ [B],  where  yq  is  the  eigenvalue  shift. 

Dynamic  matrix  used  for  the  eigenvalue  solution. 

Structural  damping  matrix  in  its  equivalent  viscous 
form  (diagonal) 

k(cs/S)  [Cjjj] , where  [C^]  is  the  imaginary  part  of 
the  AIC  matrix  [C^] . 

2K2 (cs/S) [Cj^] , where  [C^]  is  the  real  part  of  the 
AIC  matrix  [C^]. 

-K[D] (hp) , where  K is  the  rate  gyro  or  accelerometer 
gain,  D is  the  differention  or  interpolation  vector, 
and  hp  are  the  flexible  modes. 

-K(DJ [hR] , where  K is  the  rate  gyro  or  accelerometer 
gain,  D is  the  differentiation  or  interpolation 
vector,  and  hR  are  the  rigid  body  modes. 

Partition  of  generalized  mass  or  aerodynamic  matrix 
used  in  subroutine  GENMAT. 

Matrix  of  generalized  aerodynamics  as  output  from 
subroutine  GENMAT  - imaginary  part. 

Matrix  of  generalized  aerodynamics  as  output  from 
subroutine  GENMAT  - real  part. 


GENM 

(NFR.NM0OE) 

P(M31) 


Generalized  mass  matrix 


i. 


PSK(NK) 

PSN(NM) 

RI (NTT) 
RR(NTT) 

SAFI 

(NFR.NMgDE) 

SAFR 

(NFR.NMgDE) 
SI (M31 ,M3) 
SM(NM1,NM) 
SMK(M31,M31) 

SK(MFM) 

SPK(NK1,NK) 

STAB(M0DE,6) 

TO(M1,NDF) 

nang 

(NDF ,NM0DE) 
XO(NT,NT) 

XI (NT, NT) 


Row  vector  for  spline  interpolation  used  in  sub- 
routine SPLINE. 

Vector  of  LaGrangian  coefficients  from  reduced 
frequency  spline  interpolation 

Vector  of  LaGrangian  coefficients  from  Mach  num- 
ber spline  interpolation 

Imaginary  part  of  eigenvalue 

Real  part  of  eigenvalue 

Matrix  of  generalized  aerodynamic  forces  for  esti- 
mated reduced  frequency  obtained  from  spline  inter- 
polation - imaginary  part. 

Matrix  of  generalized  aerodynamic  forces  for  esti- 
mated reduced  frequency  obtained  from  spline  inter- 
polation - real  part. 

Matrix  of  'ones'  for  spline  interpolation  - used 
in  subroutine  INTERP. 

Constant  spline  interpolation  matrix  for  Mach  num- 
bers. 

Matrix  of  Mach  numbers  or  reduced  frequencies  in 
the  constant  portion  of  the  spline  interpolation 
matrix  - used  in  subroutine  INTERP. 

Generalized  stiffness  matrix  (diagonal). 

Constant  spline  interpolation  matrix  for  reduced 
frequencies . 

Stores  stability  results  for  each  mode. 

Intermediate  array  used  in  subroutine  GENMAT. 

Matrix  of  flexible,  rigid  body  and  control  surface 
modes. 

Coefficient  matrix  of  second  order  variables,  dis- 
placement. 

Coefficient  matrix  of  second  order  variables,  velo- 
city. 


X2(NT,NT)  Coefficient  matrix  of  second  order  variables,  ac- 

celeration. 


i 
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J 


J 


YO(NT,NSE)  Coefficient  matrix  of  first  order  variables,  dis- 
placement . 

Yl(NT,NSE)  Coefficient  matrix  of  first  order  variables,  velo 
city. 

Z0(NT,NSE)  Coefficient  matrix  of  zero  order  variables,  dis- 
placement. 


Z 


Dynamic  dimensioning  array  in  BLANK  C0W40N. 


S.  Complex  Arrays  (Variable  dimensions  indicated  in  parentheses) 


C2(NTT) 

C3(HTT) 

DETL  (NDF) 

E(NTT) 

U (NTT, NTT) 

V(l,l) 

VEC(M0DE,NTT) 

WK(NTT) 


Complex  working  storage  for  subroutine  EGNVCT. 

Complex  eigenvector  for  a particular  mode, 
computed  by  subroutine  EGNVCT. 

Structural  deflections  for  the  system  mass  points 
(determined  by  input  modal  data  and  eigenvectors) 

Complex  eigenvalues  resulting  from  the  eigenvalue 
solution  in  subroutine  EIGRF. 

YCA]  + [B],  where  y is  the  complex  eigenvalue  for 
a particular  mode;  used  in  subroutine  EGNVCT  to 
determine  eigenvector. 

Pauedo  complex  eigenvector  storage  in  subroutine 
EIGRF  (not  used). 

Scores  eigenvectors  for  each  mode  (complex). 

Work  area  in  subroutine  EIGRF . 


Processing  and  Programming  Considerations 
1.  Operation 

Standard  FORTRAN  IV  processor  system.  Operable  on  the  CDC  com- 
puter; model  6600,  Cyber  175-Scope  3.4.5  system. 

Note:  MPASES  can  be  made  operable  on  all  computer  systems  with 
minor  modifications.  Probable  necessary  changes  are 
listed  below: 


Deletion  of  PR0GRAM  statement  at  beginning  of  'MAIN* 
section  of  program. 

END-0F-FILE  (E0F)  statement  alterations. 

Alphanumeric  modifications  dependent  upon  the  number  of 
characters  per  word  in  the  operating  system  used. 

Single  precision  to  double  precision  accuracy.  (MPASES, 
as  presented  in  this  report,  has  single  precision 
accuracy. } 


2. 


Core  Storage 

Number  of  memory  units  (words,  bytes,  etc.)  required  to  execute 
is  dependent  upon  input  data  reflected  in  the  length  of  BLANK 


COftMON.  Estimate  of  BLANK  COKWON  length  can  be  accomplished  as 
follows : 


Servo-Elastic  Analysis: 

LENGTH  « 1 ♦ NTT(4NTT  ♦ 2M0DE  ♦ 9)  ♦ NDF (NM0DE  ♦ 2)  ♦ 6M0DE 


LENGTH  ■ 1 ♦ NDF (NDF+NM0DE  ♦ Ml)  ♦ NFR(NM0DE)  ♦ 4NFM  + (M2)2 
Use  the  larger  of  the  above  two  estimates. 

Aero-Servo-Elastic  Analysis: 

LENGTH  - 11  ♦ NTTC4NTT  ♦ 2M0DE  ♦ 9) 

♦ 2(NFR)(NM0DE)(NK  ♦ NM  + 1)  ♦ 2(M3  * 1) 2 

♦ NM(NM  ♦ 2)  ♦ NK(NK  ♦ 2)  ♦ NDF (NM0DE  ♦ 2) 

+ 6M0DE 


LENGTH  ■ 1 ♦ NDF (4NDF  ♦ NM0DE  + Ml)  ♦ 2 (NFR) (NM0DE)  ♦ (M2)2 

Use  the  larger  of  the  two.  If  a number  greater  than  NTT  was 

assigned  to  the  parameter  'MODE,'  use  MODE*NTT  in  the  above 

calculations.  Note  that  the  first  estimate  for  each  analysis 

is  eigenvalue  problem  size  oriented,  while  the  second  is 

dependent  upon  the  number  of  system  mass  points. 

Refer  to  this  section  under  Programming  Symbols  for  definitions. 

Actual  BLANK  COMMON  length  required  for  the  analysis  follows  the 
stability  results  in  the  program's  printed  output.  This  information 
will  enable  the  user  to  make  more  efficient  use  of  core  storage  in 
subsequent  analyses  by  recompiling  the  'MAIN'  section  with  a more 
realistic  BLANK  COMMON  length. 

Note.  MPASES,  as  presented  in  this  report,  contains  a BLANK  COMMON 
length  of  10.000.  This  length  was  more  than  adequate  to  per- 
form the  analyses  in  the  sample  problems  used  as  examples, 
(see  Section  4.0.) 

3.  Auxiliary  Files 

Standard  input  read  file  (5) . 

Standard  output  print  file  (6). 

Four  temporary  utility  files  (1,  2,  3,  and  4). 


i mi 


Note:  To  avoid  the  cumbersome  handling  of  data  cards,  it  is  sug- 
gested that  data  sets  be  generated  containing  the  AIC’s, 
nodal  data,  etc.  Consequently,  MPASES  can  be  modified  to 
accept  this  input  from  auxiliary  files  in  lieu  of  cards. 


Input  Instructions 

Units  - All  units  are  taken  in  the  pound- feet - second  system  with  the  ex- 
ception of  density,  which  must  be  in  slugs/cubic  feet.  The  weight,  which 
is  input  in  pounds,  is  converted  to  slugs  internally  in  the  program.  Like 
wise,  MPASES  converts  the  frequency  from  Hertz  to  ips  units  within  the 
program.  The  AIC's  must  be  non-dimensional  when  input  with  no  altitude 
consideration  (i.e.,  p*l). 

Data  Deck  Setup 


I 


I 


1.  Title  cards  (2). 

2.  Control  card  describing  data  input. 

3.  Eigenvalue  shift,  y . 

4.  Geometric  properties  of  air  vehicle  (s,  c,  S)  and  Mach  number 
deviation  (DM) . 

5.  Mach  numbers  and  reduced  frequencies,  i.e.,  for  AIC  input. 

6.  Density  and  speed  of  sound  for  each  altitude. 

7.  Number  of  velocities  to  be  analyzed  for  each  altitude. 

8.  Gain  factors,  if  gain  option  is  to  be  exercised. 

9.  Weight  matrix. 

10.  Rigid  body  control  surface  modes. 

11.  Structural  damping  coefficients  if  NFM  > 0. 

12.  Frequencies  of  vibration  (flexible)  modes  if  NFM  > 0. 

13.  Flexible  modes,  if  any. 
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14.  Restrained  (rigid  body)  modes,  if  any. 


1 

IS. 

Coefficients  of  differential  equations  describing  the  servo 
elements  in  the  control  system. 

16. 

AIC  matrices. 

17. 

Mach  numbers  to  determine  velocities  to  be  analyzed  at  each 
altitude. 

Items  4,  S,  6,  7,  16  and  17  are  included  for  aero-servo-elastic 
analysis  only. 

Detailed  description  of  input  data  follows. 


Input  Data  Description 

FORTRAN 


NO. 

CARD 

FORMAT 

COLUMNS 

NAME 

DESCRIPTION 

1. 

TITLE 

18A4 

1-72 

Any  alphanumeric  statement. 

2. 

TITLE 

18A4 

1-72 

Any  alphanumeric  statement. 

NOTE: 

Two  cards 

must  be  input.  May  be  blank. 

3. 

CONTROL 

1415 

1-5 

MC0DE 

» 1,  diagonal  weight  matrix 
input . 

■ 2,  coupled  weight  matrix 

6-10 

NDF 

input . 

Number  of  mass  points  in  the 

11-15 

NFM 

complete  system  - degrees  of 
freedom  (NDF  > 0) . 

Number  of  flexible  modes 

16-20 

NRM 

(NFM  >0).  If  NRM  • 0,  NFM 

> o).- 

Number  of  rigid  body  modes 

21-25 

NC 

(NRM  >0).  If  NFM  ■ 0,  NRM 
> 0).” 

Number  of  rigid  body  control 

26-30 

NSE 

surface  modes;  i.e.,  number 
of  control  surfaces  (NC  > 0) . 
Total  number  of  servo  elements 

in  the  control  system,  includ- 
ing rate  gyros  and  accelero- 
meters (NSE  > 0). 

I 

i 


| 

r 


i 

! 

i 
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FORTRAN 

NO.  CARD  FORMAT  COLUMNS  NAME  DESCRIPTION 

31-35  NA  ■ 0,  no  aerodynamics,  i.e., 

servo-elastic  analysis  only. 

> 0,  AIC  input,  i.e.,  aero- 
servo-elastic  analysis.  Input 
as  follows:  NA  * number  of 
partitions  contained  in  each 
AIC  matrix. 

NOTE:  NA  - 1,  AICs  derived  from  Doublet-Lattice, 
Mach  Box  theories,  etc. 

NA  > 1,  AICs  derived  from  Strip,  Piston 
Theories,  etc.  (NA  » number  of 
strips) . 

It  is  assumed  that  all  AICs  are  obtained 
from  the  same  method  for  any  one  aero- 
servo-elastic  analysis. 

36-40  NM  Number  of  Mach  numbers  for 

which  there  is  AIC  input 
(NM  < 5).  If  NA  - 0,  NM  « 0. 

41-43  NK  Number  of  reduced  frequencies 

(V)  for  each  Mach  number 
(NK  < 10).  If  NA  « 0,  NK  - 0. 


NOTE:  The  reduced  frequencies  must  be  the  same 
for  each  Mach  number. 

46-40  ND  Number  of  altitudes  (ND  < 5) . 

If  NA  - 0,  ND  - 0. 

51-35  MODE  Number  of  modes  requested  for  the 

stability  analysis.  Use  a large 
number  if  all  modes  Inherent  to  the 
analysis  are  to  be  considered. 


NOTE:  Due  to  the  nature  of  the  'lining  up  the 

frequencies 1 method  of  analysis,  the  eigen- 
value solution  may  result  in  less  modes 
than  requested. 


56-60  NG 


61-65  NGS 


■ 1.  no  gain  variation. 

> 1,  gain  option  executed. 
Input  as  follows:  NG  * number 
of  gain  variations  (K). 

(I  < NG  < 25) 

■ 0,  no  gain  variation. 

> 0,  NGS  ■ number  of  the  par- 
ticular servo  element  for 
which  the  gain  varies. 

< 0,  NGS  ■ number  of  the  con- 
trol surface  for  which  the 
gain  varies. 
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DESCRIPTION 


FORTRAN 

CARD  FORMAT  COLUMNS  NAME 


NOTE:  The  gain  factor  (K)  can  vary  for  only  one 
servo  element  or  control  surface  in  any 
one  stability  analysis. 

66-70  NAB  ■ 0,  A and  B matrices  not 

printed. 

■ l,  A and  B matrices  printed 
by  rows. 

NOTE:  If  NA  ■ 0,  [A]  and  [B]  as  generated  for 
the  analysis  are  printed. 

If  NA  > 0,  [A]  and  [B]  as  generated  for 
k*0  only  are  printed. 

SHIFT  E12.0  1-12  GAMMA  Shift  eigenvalue  Cyq) • 

Choice  of  y is  left  to  the 
o 

user. 

NOTE:  See  Section  2 (Theoretical  Development) 
and  Section  4 (Sample  Problems)  for 
details. 

NA  ■ 0,  OMIT  the  following  cards:  5,  6,  7,  8 and  9. 


AER0  4E12.0 

1-12 

SS 

Semi-span  (s),  feet. 

CONSTANTS 

13-24 

CBAR 

Mean  aerodynamic  chord  (c) , 
feet. 

2S-36 

S 

Surface  area  (S) , sq.ft. 

37-48 

DM 

Mach  number  deviation.  Choice 
of  DM  is  left  to  the  user. 

NOTE:  If  the  difference  between  the  Mach  number 
for  which  the  analysis  is  to  be  performed 
and  the  Mach  number  for  which  there  is  AIC 
input  is  greater  than  DM,  Mach  number  inter- 
polation of  the  AICs  will  be  made.  If 
NM  ■ 1,  no  interpolation. 


AIC 

DATA 

5E12.0 

XM(I) 

If  NA  - 0,  OMIT  Card  5. 

Mach  number;  I » 1,  NM. 

AIC 

DATA 

6E12.0 

XK(I) 

If  NA  « 0,  OMIT  CArd  6. 

Reduced  frequency,  k;  I»l, 

NK.  Continue  on  next  card 
if  necessary. 

NO. 


CARD 


FORTRAN 

COLUMNS 

NAME 

DESCRIPTION 

1-12 

DENS(l) 

Air  density, P , associated 
with  the  first  altitude  at 
which  the  analysis  is  to  be 
made;  slugs/cu.ft. 

13-24 

S0S(1) 

Speed  of  sound  associated 
with  the  first  altitude; 
ft/sec. 

25-36 

DENS (2) 

Density  for  second  altitude. 

37-48 

S0S(2) 

Speed  of  sound  for  second 
altitude. 

9.  AER0  DATA  SIS 


DENS (I) , S0S(I) ; I»l,  ND. 
Continue  on  next  card  if 
necessary. 

If  NA  * 0,  OMIT  Card  8. 

MD(I)  Number  of  velocities  (deter- 

mined by  Mach  number  input  - 
see  Card  26)  for  each  alti- 
tude; I - 1,  ND.  (MD  * 10) 

If  NA  - 0,  OMIT  Card  9. 


10.  GAIN  6E12.0 


11.  WEIGHT  6E12.0 


GAIN(I)  Gain  factor,  K.  I * 1,  NG. 

Continue  on  next  card  if 
necessary.  (NG  < 25) 

If  NGS  * 0,  OMIT  Card  10. 

WT(I,I)  If  MC0DE  * 1,  vector  of  di- 
agonal elements  of  weight 
matrix,  pounds.  1*1,  NDF. 
Continue  on  next  card  if 
necessary. 

WT(I, J)  If  MC0DE  • 2,  elements  of 
upper  triangle  of  coupled 
weight  matrix,  lbs.  Input 
by  rows;  each  row  starts  on 
a new  card.  1*1,  NDF, 

J - I,  NDF. 

Continue  on  successive  cards 
to  complete  each  row. 
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NO.  CARD 


FORMAT  COLUMNS 


FORTRAN 

NAME 


DESCRIPTION 


12.  M0DE  6E12.0 


13.  DAMPING  6E12.0 


14.  FREQUENCY  6E12.0 


IS.  M0DE  6E12.0 


16.  M0OE  6E12.0 


CM®DE(J,I)  Rigid  body  control  surface 
modes.  Each  mode  starts  on 
a new  card.  I ■ 1,  NC,  J ■ 1, 
NDF. 

Continue  on  successive  cards 
to  complete  each  mode. 

GF(I)  Structural  damping  coefficient 

corresponding  to  vibration 
mode  (I).  1*1,  NFM.  Con- 

tinue on  next  card  if  necessary. 


If  NFM  « 0,  OMIT  Card  13. 

FREQ (I ) Frequency  of  vibration  mode 

(I),  Hz . I • 1,  NFM.  Con- 
tinue on  next  card  if  necessary. 


If  NFM  « 0,  OMIT  Card  14. 

FM0DE(J,I)  Free-free  vibration  (flexible) 
modes.  Each  mode  starts  on  a 
new  card.  1*1,  NFM,  J * 1, 
NDF.  Continue  on  successive 
cards  to  complete  each  mode. 


If  NFM  - 0,  OMIT  Card  15. 

RM0DE(J,I)  Restrained  (rigid  body)  modes. 

Each  mode  starts  on  a new  card. 
1*1,  NRM,  J * 1,  NDF.  Con- 
tinue on  successive  cards  to 
complete  each  mode. 

If  NRM  • 0,  OMIT  Card  16. 


Input  description  of  the  servo  differential  equation  coefficients  follows: 

IMPORTANT:  All  servo  differential  equations  must  be  equated  to  zero, 
i.e.,  all  terms  should  be  on  one  side  of  the  equation,  be- 
fore attempting  to  input  the  coefficients. 

NOTE:  Reference  to  servo  elements  and  control  surfaces  is  by  num- 
ber; therefore,  numbers  should  be  assigned  to  each  servo 
element  (1  to  NSE)  and  control  surface  (1  to  NC) . Servo 
elements  also  include  the  rate  gyTos  and  accelerometers  in 
the  control  system. 


i 


CARD  FORMAT 


COLUMNS 


FORTRAN 

NAME 


DESCRIPTION 


REPEAT  the  following  cards  17,  18  and  19  for  each  control  surface,  i.e., 
Nd  times. 

17.  Output  from  Control  Surface.  FORMAT  (3E12.0) 

Columns  1-12  X0(I)  Zero  order  coefficient  for  Control  Sur- 


Columns  1-12  X0(I)  Zero  order  coefficient  for  Control  Sur- 
face (I) . 

13-24  XI (I)  First  order  coefficient  for  Control  Sur- 

face (I). 

25-36  X2(I)  Second  order  coefficient  for  Control 

Surface  (I). 

Note:  1-1,  NC. 

18.  Input  Control  Card.  FORMAT  (IS) 


Columns  1-S  INC 


Number  of  inputs  from  servo  elements  to 
control  surface  (I) . 


REPEAT  card  19  for  input  from  each  servo  element  to  Control  Surface  (I), 
i.e. , INC  times. 

19.  Input  from  Servo  Element.  FORMAT  (2IS,2X,3E12.0) 


Columns  1-S  K 
6-10  N0 
13-24  CO 
2S-36  Cl 
37-48  C2 


Servo  element  from  which  there  is  input. 
Order  of  servo  element  coefficients. 

Zero  order  coefficient. 

First  order  coefficient. 

Second  order  coefficient. 


If  N0  » 0,  zero  order;  only  CO  input. 

■ 1,  first  order;  Cl  must  be  input  (CO  may  be 
zero. 

• 2,  second  order;  C2  must  be  input  (CO  and 
Cl  may  be  zero) . 

REPEAT  the  above  card  19  INC  times. 

REPEAT  the  above  cards  17,  18  and  19  NC  times. 

REPEAT  the  following  cards  20,  21  and  22  for  each  servo  element,  i.e., 
tfSfe  times  (not  necessarily  in  sequential  order). 


FORTRAN 

NO.  CARD  FORMAT  COLUMNS  NAME  DESCRIPTION 


20.  Output  from  Servo  Element.  FORMAT  (2I5,2X,3E12.0) 


Columns  1-5 

I 

Servo  element  number  from  which  there  is 

6-10 

N0 

output . 

Order  of  servo  element  coefficients. 

13-24 

CO 

Zero  order  coefficient 

25-36 

Cl 

First  order  coefficient. 

37-48 

C2 

Second  order  coefficient. 

If  N0  - 0, 

only  CO  is  input. 

■ 1,  Cl  must  be  input  (CO  may  be  zero). 

■ 2 , C2  must  be  input  (CO  and  Cl  may  be  zero) . 

21.  Input  Control  Card.  FORMAT  (IS) 

Columns  1-5  INS  Number  of  inputs  from  servo  elements  (or 

control  surfaces)  to  servo  element  (I). 

If  there  are  no  inputs  to  Servo  Element 
(I),  INS  « 0. 

REPEAT  card  22  for  input  from  each  servo  element  or  control  surface  to 
Servo  Element  (I),  i.e.,  INS  times. 

If  INS  • 0,  OMIT  Card  22. 


22.  Input  from  Servo  Element  or  Control  Surface.  FORMAT  (2IS,2X,3E12.0) 
Columns  1-5  K 

-K 

6-10  N0 

£ 

| 

N0  ■ -1,  input  from  body  angular  rate  to  rate  gyro. 
■ -2,  input  from  body  acceleration  to  accelero- 
meter. 

NOTE:  If  N0  < 0,  K ■ 0. 

13-24  CO  Zero  order  coefficient. 

25-36  Cl  First  order  coefficient. 

37-48  C2  Second  order  coefficient. 


Servo  element  number  from  which  there  is 
input . 

Control  surface  number  from  which  there 
is  input. 

If  K > 0,  N0  ■ order  of  servo  element  co- 
efficients. 

If  K < 0,  N0  ■ order  of  control  surface 
coefficients. 


j 

\ 

i 
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NO.  CARD 


FORMAT  COLUMNS 


FORTRAN 

NAME 


DESCRIPTION 


If  N0  ■ 0,  only  CO  input. 

• 1,  Cl  must  b«  input  (CO  may  be  zero) 

» 2,  C2  must  be  input  (CO  and  Cl  may  be  zero) 

■ -1,  only  CO  input,  where 

CO  ■ -Kg  (rate  gyro  gain) 

■ -2,  only  CO  input,  where 

CO  ■ -K  (accelerometer  gain) 


NOTE:  When  N0  < 0,  Card  22  must  be  followed  by  Card  23. 
23.  If  N0  • -1, 

Differentiation  Row  Vector  used  to  determine  angular  velocity  for 
the  rate  gyro. 

If  m « -2, 

Interpolation  Row  Vector  used  to  describe  body  acceleration  for  the 
accelerometer. 

FORMAT  6E12.0  D(I)  Element  of  differentiation  or  interpolation  row 

matrix.  I • 1,  NDF. 

Continue  on  next  card  if  necessary. 

NOTE:  See  Ref.  1 for  description. 

REPEAT  the  above  Card  22  (and  Card  23,  if  applicable)  INS  times. 

REPEAT  the  above  Cards  20,  21  and  22  NSE  times. 


Input  instructions  for  the  AIC  matrices  [C^]  follow: 

NOTE:  It  is  assumed  that  the  AlC's  are  dimensionless  and  do  not 
reflect  density  (i.e.,  p ■ 1). 

If  NA  • 0,  OMIT  Cards  24  and  2S . 

If  NA  « 1,  input  only  Card  25. 

If  NA  >1,  input  Cards  24  and  25. 

The  AICs  are  input  for  all  the  reduced  frequencies,  i.e.,  for  XK(I), 
where  I • 1,  NK  for  each  Mach  number  XM(J),  where  J « 1,  NM. 

NOTE:  The  reduced  frequencies,  k,  must  be  the  same  for  each  Mach 
number. 

If  NA  ■ 1,  REPEAT  Card  25  for  each  AIC  matrix,  i.e.,  (NK*NM)  times. 


NO. 


CARD 


FORMAT  COLUMNS 


FORTRAN 

NAME 


DESCRIPTION 


T 


1 


( 


If  NA  > 1,  REPEAT  Cards  24  and  2S  NA  times  for  each  AIC  matrix. 
See  Card  3.  NA  ■ Number  of  partitions  (NPART)  in  the  AIC  matrix. 
All  the  input  AIC's  must  have  the  same  number  of  partitions. 

REPEAT  until  all  AIC's  are  input,  i.e.,  (NK*NM)  times. 


24.  Partition  Control  Card.  FORMAT  (2IS) 

Columns  1-5  MS  Size  of  partition. 

NOTE:  The  sum  of  the  sizes  (MS)  of  all  the 
partitions  in  each  AIC  matrix  must  equal  NDF. 

6-10  NZER0  ■ 1,  all  elements  of  partitions  are  equal 

to  zero.  Do  not  input.  OMIT  Card 
25. 


■ 0,  non- zero  partition.  Input  Card  25. 


25.  AIC  Partition  or  Matrix.  FORMAT  (6E12.0) 


CH(I,J)  If  NA  » 1,  complete  AIC  matrix  (complex) . 

I » 1,  NDF,  J * 1,  NDF. 

If  NA  > l,  AIC  partition  (complex). 

I » 1,  MS,  J ■ 1,  MS. 


Columns 


Input  by  rows.  Each  row  starts  on  a new  card.  Continue 
on  successive  cards  to  complete  each  row. 

The  imaginary  part  of  each  element  follows  the  real  part. 


e.g. : 

1-12 

CH(l.l) 

Real 

13-24 

CH(1,1) 

Imaginary 

25-36 

CH(1 ,2) 

Real 

37-48 

CH(1,2) 

Imaginary 

OMIT  the  following  Card  26,  if  NA  ■ 0. 

REPEAT  Card  26  for  each  altitude,  i.e.,  ND  times. 


h 

! 

(• 
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FORTRAN 

NO.  CARO  FORMAT  COLUMNS  NAME  DESCRIPTION 


26.  Mach  Numbers  for  Analysis.  FORMAT  (6E12.0) 

XMACH(I)  Mach  numbers  selected  for  stability  analysis  at  a 
particular  altitude.  I » 1,  NMACH,  where  NMACH  is 
the  number  of  velocities  (determined  by  the  Mach 
number,  XMACH(I)  and  the  speed  of  sound,  S0S(J)  for 
altitude  J),  i.e.,  NMACH«MD(J) . Refer  to  Card  9. 
(fWACH  £ 10) . Continue  on  next  card  if  necessary. 
No  eigenvectors  are  obtained  or  printed  with  the 
stability  analysis  when  XMACH(I)  > 0. 

-XMACH(I)  Same  as  above. 

Complex  eigenvectors  are  calculated  and  printed  for 
each  mode  when  XMACH(l) < 0. 

NOTES:  XMACH(I)  may  be  the  same  or  different  for  each 
altitude.  Do  not  confuse  XMACH  with  XM  which 
are  the  Mach  numbers  for  the  input  AIC's. 


REPEAT  above  Card  26  ND  times. 


Program  Output  Description 


( 


Input  Data 

1.  Upper  triangle  of  weight  matrix  (lbs.). 

2.  Damping  coefficients  and  frequencies  (if  free-free  vibration 
modes  are  present) . 

3.  Mode  shapes  - flexible,  rigid  body  and  control  surface. 

4.  Control  system  description  - coefficients  from  servo  differ- 
ential equations,  rate  gyro  gains  and/or  accelerometer  gains. 

5.  Maximum  order  of  each  servo  element  (determined  internally 
in  the  program) . 

6.  Rows  and  columns  (assigned  by  the  program)  in  the  A and/or 
B matrices  for  the  coefficients  of  each  servo  element. 

7.  Eigenvector  element  associated  with  each  servo  element  (velo- 
city and/or  displacement). 

8.  AIC  matrix  (by  rows)  for  each  reduced  frequency  and  Mach 
number,  if  aero-servo-elastic  analysis. 

A and  B Matrices  (if  requested) 

1.  For  an  aero-servo-elastic  analysis,  only  [A]  and  [B]  for 
k»0  are  printed. 

2.  For  a servo-elastic  analysis,  [A]  and  [B]  as  generated  for 
the  analysis  are  printed. 

Stability  Analysis  Results 

The  results  are  identified  and  tabulated  for  each  gain  factor 

variation  (if  any).  In  addition,  when  aerodynamics  are  included, 

the  results  are  printed  for  the  velocities  analyzed  at  each 

altitude. 

1.  For  each  mode  calculated,  the  following  is  printed  In  tabulated 
form: 
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Complex  eigenvalue,  y - real  part,  m>,  1/aac. 

Coop  lax  eigenvalue,  y - imaginary  part,  u>,  rad/sac. 
Damped  frequency,  Hz. 

Undamped  fraquancy,  Hz. 

Raducad  fraquancy,  k,  uaad  to  datarmina  aiganvalua 
problam  for  particular  aoda. 

fraction  of  critical  dashing,  Q . 

Tima  to  half  amplitude,  sac. 

Complex  eigenvectors  corresponding  to  each  aiganvalua,  if 
requested  for  the  aaro-aarvo-alaatic  analysis  of  a particular 
velocity.  The  eigenvectors  are  always  printed  with  the  sta- 
bility results  of  a servo-elastic  analysis. 

Structural  deflections  (complex)  representing  the  system  mass 
points  for  each  mode  are  printed  whenever  eigenvectors  are 
calculated . 


The  minimum  BLANK  C0MM0N  length  required  to  execute  the  program. 
This  length  is  dependent  upon  the  input  data. 


u 


2. 


3. 


! 

I 

SECTION  IV 
SAMPLE  PROBLEMS 

As  examples  of  Che  program  capabilities  we  consider  a simplified  I 

i 

missile  and  servo  configuration  first.  Case  1,  as  a rigid  body  in  vacuo, 
second.  Case  2,  as  a flexible  body  in  vacuo,  and  third.  Case  3,  as  a flex 
ible  body  in  an  incompressible  airstream  at  sea  level.  The  missile  is 
idealized  as  a uniform  beam  and  control  surface  as  shown  in  Fig.  1. 

k 

h 

£ | 

I 

y 

i 

Figure  1.  Uniform  free-free  missile. 

The  modal  characteristics  and  differentiation  matrix  for  this  configuration  ; 

, I 

are  given  in  Ref.  11.  For  the  numeridkl  work  we  assume  the  missile  to  weigh 
10m  * 1000  lbs.,  to  be  ISO  inches  long,  and  to  have  a fundamental  frequency 
of  45  Hz.  From  the  data  in  Ref.  11.  this  results  in  higher  frequencies  of 

I 

{2  * 12S. 4 Hz  and  fj  ■ 248.2  Hz;  only  three  modes  will  be  considered.  The 
structural  dampings  in  the  three  modes  are  assumed  to  be  gj  « 0.03,  gj  * 0.05,  ' 

and  gj  ■ 0.08. 

A simple  servo  system  is  considered  whose  block  diagram  is  shown  in 
Fig.  2. 
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Figure  2.  Servo  block  diagram. 

We  assume  the  transfer  functions  appear  as  follows. 
For  the  command  error, 

e.  * e.  ■ e^ 

1 4 c 

for  the  servoposition  error, 

•l  “ ®2  • *3  ‘ 0 

for  the  actuator, 

e3/6  - Kas/(Tas  + 1) 
for  the  position  potentiometer, 

.j/S  • Kp 

and  for  the  rate  gyro. 


•4/h'  * Kg/(32/u)g2  ♦ 2^gs/wg  * 


The  differential  equations  corresponding  to  Eqs.  (78)  and  (80)  are 

Ta*3  + e3  ' V “ 0 

and 


(X/“g2)  «4  * (25g/wg)  * e4 


i 


r 


PL 


The  actuator  is  seen  to  be  a first  order  servo  element  and  the  rate  gyro  is 
a second  order  element.  Therefore  the  second  order  variables  are 


(83) 


the  first  order  variable  is 

(y)  " {*3} 

and  the  zero  order  variables  are 


(84) 


M - CM 


(8S) 


For  the  numerical  work  we  assume  Ka  » 1/6  deg.  per  deg. /sec.,  Tft  ■ 0.01  sec., 
Kp  - 1 deg. /deg.,  oig  » 376.991  rad/sec.  (60  Hz),  Cg  ■ 0.70,  and  Kg  - 0.3  deg. 
per  deg. /sec.  The  differentiation  matrix  for  the  rate  gyro  from  Ref.  11  is 

[D]  - (1/241)  (0  0 +1  -27  +27  -1  0 0 0 0 ] (86) 

where  1 ■ 15  inches. 

The  case  descriptions  below  present  the  variables  in  the  eigenvectors 
and  discuss  the  solutions.  The  data  input  code  sheets  and  the  program 
printed  output  for  the  three  sample  problems  appear  following  the  case  de- 
scriptions. 


u 

[• 
» JS 
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Case  1 Description 

The  first  case  is  the  rigid  body  in  vacuo.  There  are  two  rigid  body 
modes  of  plunging  and  pitching.  The  eigenvalue  calculation  is  made  choosing 


LJ 
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y * -1.0.  The  11—  order  eigenvector  that  appears  is  in  the  following 
o 

order: 

{V,T  ■ l**l  *»2  1 *4  *R1  *R2  4 *4  *3  *1  <87) 

There  are  four  zero  eigenvalues  from  the  two  rigid  body  modes  and  four 
infinite  eigenvalues  from  the  command  error,  the  servoposition  error,  the 
position  potentiometer,  and  the  actuator.  The  three  non-trivial  eigen- 
values consist  of  the  real  actuator  damping  and  the  complex  conjugate  rate 
gyro  frequency  and  damping.  However,  since  only  the  eigenvalues  with 
positive  frequencies  in  the  conjugate  pairs  are  printed  in  the  stability 
results,  only  nine  solutions  are  presented  in  the  printed  output.  The 
structural  deflections  at  the  ten  system  mass  points  for  each  mode  are  also 
printed. 


Case  2 Description 

The  second  case  is  the  flexible  body  in  vacuo . The  addition  of  three 
flexible  modes  to  Case  1 results  in  a 17—  order  eigenvector  which  has  the 
variables  printed  out  in  the  following  order: 

{Y>T  « lapi  ip2  ^F3  ^R1  *R2  ^ ®4  ®F1  *F2  *F3  aRl  *R2  3 *4  ®3  ®1  ®2J  (88) 

The  eigenvalues  are  obtained  by  again  choosing  y • -1.0.  Nine  non-trivial 

l ° 

eigenvalues  are  obtained  in  this  case.  Three  are  those  obtained  in  Case  1 
and  the  additional  three  complex  conjugate  pairs  correspond  to  the  three 
flexible  modes.  The  stability  results  present  the  modes  with  the  frequencies 
in  ascending  order,  the  negative  values  being  omitted.  Therefore,  only 
thirteen  solutions  are  printed;  the  structural  deflections  for  the  thirteen 
modes  are  also  printed. 


0 


r 
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An  unsatisfactory  design  is  seen  to  exist  from  the  negative  damping 
in  the  133.0  Hz  mode;  the  missile,  will  "buzz"  in  this  mode. 

Case  3 Description 

The  third  case  is  the  flexible  body  in  an  incompressible  flow  at  sea 
level  with  a density  p • 0.00237692  slugs/cu.  ft.  and  velocity  V « 500  fps. 

It  is  the  same  as  Case  2 with  the  addition  of  aerodynamic  loads.  The 
aerodynamic  loads  are  assumed  to  act  only  on  the  control  surface  and  are 
derived  from  the  incompressible  strip  theory  presented  in  Ref.  15.  The 
control  surface  has  an  exposed  span  of  20  inches  and  a semichord  of  15  inches. 
Five  reduced  frequencies  are  chosen  for  the  aerodynamic  interpolation, 
k • 0.0S,  0.10,  0.20,  0.50,  and  1.00;  the  minimum  value  of  k * 0.05  is 
chosen  because  of  the  singularity  in  the  aerodynamic  damping  at  zero  fre- 
quency when  it  is  computed  from  Theodorsen’s  function. 

Although  there  are  17  degrees  of  freedom  in  the  eigenvalue  problems 
solved,  only  nine  modes  can  be  obtained.  This  results  from  the  frequency 
"lining-up"  process.  The  reduced  frequencies  determined  by  the  estimated 
frequency  in  each  eigenvalue  solution  are  shown  along  with  the  stability 
results  for  each  mode.  In  addition,  the  time  to  half  the  amplitude  as  well 
as  the  dating  ratio  is  indicated  to  facilitate  stability  evaluation. 

The  solutions  for  Case  3 are  seen  not  to  be  significantly  different 
from  those  of  Case  2.  However,  we  note  that  the  short  period  mode  shows 
up  in  place  of  the  rigid  body  pitching  mode.  The  reduced  frequency  of 
the  short  period  mode  is  k « 0.03346  and  its  damping  ratio  is  C ■ 0.1401. 

The  aerodynamic  loads  do  not  cause  a large  disturbance  for  the  (impractical) 
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configuration  chosen  as  the  example.  The  servoelastic  "buzz"  observed 
in  Case  2 is  still  unstable  here.  More  practical  configurations  should 
be  studied  that  have  more  critical  hinge  line  locations  and  mass  balancing 
on  the  control  surface. 
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.44  1 N9U4  0 
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lna|  W«H|  f *U 

i<*a|0...n  *» 

i 1 W*lU.Nl I HA 

1  1 t«4  1U*-«H  «A 

n<«ni-*rr  ha 

llJ«U2*Mt»»  ha 

l0tPall3.*UP»2  „ ?* 

l»al«lutP  •nhf  “**  ha 

lvalH»i«<.n  «* 

(/(na. •.*..«!. }*i  10  Ja" »a 

ha 

t»LL  ha 

1 2(  II VI  .£(1  111  .21  1UI  . UAJW  .Hll.HUUfc.HHUUC.MPn.QAHHA.HO/  .2(1131  t HA 

i2(U.tH.Ul^l)i<(i»)l  “ ■ “ “ha 

jTuTa»A*0(J|.O2«U.J.,jal  HA 

aNlTfc  (V.olOl  JlUl HA 

W lu  I HA 

3a  11a«IIJ»(«hMI-*h  "A 

l I3«l  l »•  (ilA  » 1 I "rrf HA 

H JanAAVIf.rl.rtAl  HA 

1 !••!  !*♦  »Hj«  I ) »h3  ha 

1 1 /»1 1*1*  (O*  l.l  • (H  J.  l ) HA 

“iivaiiwu  «» 

Huai  ln.nP>«*MH(JUfc»Mt  HA 

lju«l  l .»l«p  •.•NrtUufc'HA HA 

l2la|2V.MpHANHUUfc»MM  HA 

122*12 1 ••.PHANrtulifc  Ann  ■ HA 

___12j«l<2*NP«.Afinuu€. HA 

U»*lil*W'K,Af1W(t  HA 

i/aal2A«*»«  ha 

12a*l2"»*M<  ..  _ ha 

12J«l2n*HJM  »A 

10ePal2'/».'.uP»2  "A 

___l»Sl«ll»tP*HO»»HHuufe  __  __  _____  HA 

i»«ih>i.nh  ' — ha 

wS»I»*2  "A 

CALL  -HAa  (2 1 1 11  .ZIIO  .4(1  J1  «2U»t  .2  dal  .2  llal  .A  1 1 /)  .2(  lal  .2  (IV)  . ha 

1 2(  I1U1  .21 II  It  >2(1121  .41 11  J)  .2(11*1  .filial  .21  llol  .2  (11  I)  .2(ll«l  . HA 

2 2 ( I1H1.2  IUVI  .211211  .2( 1221.2  1 12J)  .21 12» I . Z ( 1 25)  . 2 U 201  .HIT, n/H.  ha 

3.  HMOOt  .MH  .»<H»  1 .«A  • HA  • 1 • *Odfc  (H3«H  J*  | , OAHHh .HOP.  (1.1  1 »2(  IQEPl.  HA 

*’((IH*1!.2I1V>)  ....  - HA 

jlufaAU»(.ilia(i<Jji<J«<g)l  ha 

anlttin.siai.y.lul HA 

(HI  to  I ha 

atl  Pohnaiuhahi  ha 

>«»«  fohhai  (i*i>t ha 

3VJ  FOHHAf (At  12.01  HA 

*•1  FvlHHA  f ( (H|  • ASA  . 1 HA4//2VA • 1 AAA// 1 HA 

••2  POhhAHIhU  «.tA.j/H4fc..U-afcKVU-fcLAaUC  aiAXlLlTT  ANALYSIS  //>  HA 

•03  PuHHAiiimi  otiMHUPwi  tLAsric  mahIliit  analysis  //  i ha 

P0n»A|(1hV  AHA.  12  • l*H  UfcO*lt»  UP  PHttuUN/AH*.  12,  13H  PLtMBLC  HOOHA 

. IfcS/AuA. I*. I /N  Hlttlu  (ftlllf  HUUtS/AHA.  12.  I /"  LONlHUL  >QH>  ACt>/  A (HU  HA 

21  J.  Ian  Se .HAW  IttHtA  fi// A<A.  25HAH  IP  » t IUCN*  *LUl  ( QAMHA  1 * IHCll.J  HA 

3 //AHA. 13* I HH  HUUtS  »EWUt»ltU  //I  HA 

•oa  /UrtHa I llHV.AUA.il. 2/H  ACllIUOt  ♦AH(AriUPIS//3aA.l?HA(C  HAfHlCta  PUHHA 

1 .11. Inn  HALH  AVHIHPJ  Anil  . I2.2VH  HtOULtU  /HtOUENCltS//*/*.  HA 

2 /2HHtPtnt'<P  atH|-i.nUH)>  a lHtl2.3./n  r l/*3«. | |n*tH|-»rtAM  a HA] 

j iHtla.a.jn  ii/viaiiaamnmci  ant  a a mti2.a.HH_a«  y//*j»i hai 

A 24HHALH  WuHNpH  Ufc/IAUUN  a >|<.l  //I  HA] 

l»«HH»r  lln0.3/A.12.A|n  VAIN  P AC  1 oh  vahIAMumS  POW  afcHVO  CLCHCHT  13  ha] 

1//1  . ..  _ ...  _ _ _ HAI 

•(II  PUH*A|  (lHV«j!A.l2'AJn  l»«iN  PAClOH  VAHlAliuNH  PUN  LUNIHUL  aUN/ACk~H4| 

112  //I  HAI 

.•It  POHHA 1 ( 1HI  .2/  . JnnHlnlHUH  VLAN*  CUHrtun  LfcNvTH  HtUUlHlO  a 1».  1H./3A.HA1 

i Ajn«A»feu  u».  inhoi  uai*  mv  Anavrata  ntwwtartu.i  ~ ' "hai 

bNU  HAI 
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~rn«»  ti«r»u!in  «M 
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1 N*0  __ 
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00  14  J*l,N0P 
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i* 
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00  14  J4I.N0* 
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•**n  -con  4 no  point 

■ • 1 Tf  f *,* 1 4 1 

00  *4  !4|.NC  ■ ' 

• 1 *0(4,19))  (C*nOt (J, I )• J*l .NO*) 

*4  -•|Tf(N,*||)  T.(tNQ0C;J.n,j4),N0*) 

t((N»-.f.».9)  GO  TO  T9' 

•1*0(1,19))  (0P<l).t«).N*N) 

•(40(1.101)  (*M0tn,t4),N*N) 

~ ■•!)({*. *1)1 

0«)  *1  1 4 ) , N*N 

•7  *0(4,101)  (*-n0f(J.I).j4),N0*)  

— *«m  (n.*u)  i .*•(•< t ).&*(().  (tnoo*  t j7 i )T3M7^o() 
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TO  tr(N4N,|0.4)  60  TO  44 

N« IT* (4,4)4) 
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■*4444)4 
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4*444)40 
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4*544)90 
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»*4*01«1 
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44)50)00 
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4*4)4)90 
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4*4501)0 
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•1)  »0«4*T(///  l*,l*4PL(*t4UI  400(5  //}  

*1  J M1*4*f(|N0  )*,*44P(1*  ,t),14M  • PMOUfNCT  4 PI),), *4  CPt/tll,  ' 
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00  1 Jat.NT 
-X9(t,4>«0,9 

I X2(t.J)«0.0  

- 00  k JJM.NJf 
y«d.jj)«o,o 
y»<t.jj)«o.o 

rmr,jj»«o,a 

f eoatiNut 
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60  TO  (3 

I NOB(«K,L)aMQ 

i court Mut 
«UL.I)<t 

'HLU.21«* 

» eoartNor 

l CONTIMuF 

NtNC»L 

00  113  Jat.MJC 

“ »e*o(3,3o«j  •T.wj7eo(a;cr(U7C2(n 

aotrc(k.*o3)  t,t.c2(u.cuu.eo(u> 

Hta.nat 
""NL(L»2)al 
00  13  xxat.Nsr 
I9<XX.fO,t)  60  TO  19 

-*>■(**. L>«0  
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( COntinuC 
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00  M9  "•( • tNG 

0(40(3.399)  k.nO.CKC)  <Cl  (L)  >C2(L)  

>0(0, (.T. 9)  60  TO  (99 
K(*0.lT,9>  GO  TO  39 

_ »atTC(k.499>  X.C2(L).C)(L)aC9(L) 

OO  as  <yai .ait 

tKnx.fO.X)  60  TO  99 

_N0*(*«,U)»9 

GO  TO  as 

I NO*(ax.l )*Nfl 

( eoNttNur  

N*,(l.  t )at 

«l(l.2)aa 

. 60  TO  U0  
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HNyooio 

ICPV0020 

ijayooio 

,uavo9«o 

leavooso 

ITMVOOaO 

learooro 

navoioo 

icrvouo 

iravoico 

ICRvOt 19 
•CavOlO 
icavoiso 
tiavotao 
•ravoKo 

9(NV0 t 99 

uavotao 
Ifay9200 
ICRv02t0 
iravoixo 
K avo2io 

9(mV0240 

9tNy9230 

seavocao 

STNV0270 

9CRv02'S 

irav3290 

»t»V029« 

icavoioo 
lc»yg)io 
J£»vO 120 
9C4V3  1)0 
sravoioo 
uavoiso 
IMvoiaO 
ITuy9)70 
icayoiNO 
SI-V0190 
»|avoa«0 
icayoato 
9d«0>» 
sravoato 
ICavOdaO 
ieavo«39 
•KavOaoO 
IC»vO«?0 
■ravoaao 

' sravasoo 
iravgsto 

IERV0S29 

'irav0939 

UavoS39 

~tcav93*9' 

•e*yo3To 

If MV9**0 

t(ayoa99 

I(*«9al9 

G(NV9a20 

icavoaio 

irav9a«9 

irwyoaso 

ir*voaa9 

9tav9aro 


KavoToo 

uayona 

9(ay0T29 

KNV0T19 

IfavOTO 

»avoT39 

ItayOTaO 

ICNyOTTO 
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M aUl.liaO 

HUliilM  

t»<«o.ia,<*m  to  ro  *» 
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60  TO  at  

i%  ••i»i«*.*o*i  c«u> 

ii  pc  a 
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00  ••  

so  to  n 

»4u».jji««4u».  JJnocu.jj) 

60  Tfl  M 

66  CO*TI*Ut 

_ 66  !*(***, d.O)  60  TO  li* 

c»u  •»t'*»uo,*«ooc.o»".i.»o«,ui,»«6i.,iio*,i) 
00  «1  JJa|,«a« 

_ icaacjj  

t*{Mo,(6.(-tn  oo  to  *a 
lac  ta.  ternon*.  ic>  •own  ,jj> 

60  TO  *6 

*4  *i ( t*« ic  >»»  i mUIci  »o*an • jj) 

•9  coarimif 

60  TO  III  _ 

i ih»ut  «o*  cont»ql  iu»>C( 

166  11 (t. I la# 

NCa.IT** 

■•itJ t»* , c 1 1 u . fUU.coTU 

00  l«*  *«■!.»>( 
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I*aN«nor*I 

_______________ 

•«(!».  ie»»c#tu 
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_ »*(T»,ie>»cj(u)  

116  eONTINUt 
m eo«<Tt«iUf 

I»(tNO,r«.«)  60  TO  116 

com*  i 
ac»t*o  i 

•atrm.tm 

— Wacom  i.ionTJ-OTWlTHOfi 
t’lrormi  in. in 

III  eOMTINU*  

t»*m<*,»it)  l,<0ii733)7JJai;aoai 
00  TO  116 

II*  a«»t*o  i 

no  COnti no* 

C 0«t*»-|n*  oaoca  or  navo  wmi 

oo  i«a  aai.nsr 

oo  l>a  UM.l 

IP(NOa(N,UI.*0.<)  to  TO  111 

too  eoofiNur  

»o  149  u*»7i 

l*(NO*(t.UI.I0.H  60  TO  1)6 

1H  CONTCu*  

a(l(a.  I lat 
Man.) 

<C<«.»)as6 

60  TO  )«« 

116  KO«.U*l 

____  Ntwial  

«0<*,4nNl  ' 

60  to  | at 
IH  aOia.na> 

1 alaat.i  ~ 

•Oti.liwl 

too  eONTjNu*  ____________ 

t * TOW  I,  ».  I aiTliell 
OO  IT*  tat, I 


It  ayg  M 
l(*«9T« 

IIMIII 

IM«9ll 

StavOli 

Him* 

Slay#*) 

mom 

tfavoia 

lt«y|la 


loayiati 

lt«vl«)l 

lC'«t«)l 

K*vO*a| 

ttay0«*0 

l*a«««ao 

tiaaoara 

t(yy«a«o 
IImvi  jga 
If ay  to  to 
U»V1040 

ua.toio 
ttayioa# 
ii ay  1 0*4 
K«V|0k4 

irayioto 
Uay  i««o 
Itavioag 
ttavuot 
Ifavilio 
Itiati4i 
iriyiut 

»|ay I l«4 
***y>  no 
ITayiiaO 
Stay!  WO 

icavtti* 

6tay|  i«o 
iiayuoo 
uaviiio 
*r«yi*it 
ir ayiiii 
icavtai) 
liaviita 
il*ym* 
iravwta 
«*yillT 
Iftyiail 
Ofay  lit* 
ttay tifo 
If *yi44l 
If  ay  1444 
»fay 1 4 JO 
Of ayi4«0 
Kay  t4*0 
I(tvi4*4 
Ifavtin 
M*vi4l* 
K*yt4ao 
Kaviioi 
tiaymo 
Stay t)44 
II  ay  Dio 
uay  i no 
liay|i*a 
Itay  t la* 
Ifayiiro 

tiavt loo 
**ay| l*o 
liaaitoo 
Kaytato 
IMviall 
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APPENDIX  A 


THEORETICAL  DERIVATION  FROM  REFERENCE  I 

2.0  THE  CLOSED-LOOP  AERO-SERVO-ELASTIC 
STABILITY  PROBLEM 

2.1  NOMENCLATURE 

A Element  of  coefficient  matrix  of  |v|  in  Equati  * (2-4) 

AIC  Aerodynamic  influence  coefficient 

a Amplitude  of  generalized  displacement  coordinate 

B Element  of  coefficient  matrix  of  {v  | in  Equation  (2-4)  | 

Reference  dimension 

C Element  of  discrete  damping  matrix 

C Element  of  generalized  damping  matrix 

D Element  of  differentiation  matrix 

e Servo  signal  - 

I* 

F Element  of  forcing  function  matrix  £ 

i 

gp  Modal  structural  damping  coefficient 

H Element  of  interpolation  matrix 

h Deflection  of  aeroelastic  system 

hg  Deflection  due  to  unit  shaft  rotation 

I Element  of  unit  matrix 

i Imaginary  unit 

K Servo  system  gain  constant;  element  of  stiffness  matrix 

it  Element  of  generalized  stiffness  matrix 

k Reduced  frequency  (Strouhal  number) 

t 

M Element  of  discrete  mass  matrix  ! 

i 

i 

m Element  of  generalized  mass  matrix 

I 

Q Element  of  generalized  force  matrix 

Qa  Element  of  AIC  matrix 


J 
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Re(  ) Denotes  real  part  of  ( ) 


s 

T 

u 

V 
v 
X 
x 

V 

y 

z 

z 

V 

^o 

6 

S 

\ 

\ 

o 

« 

C) 

Subscripts 

a 

F 

g 

h 


Laplace  transform  parameter 
Servo  system  time  constant 
u = x 

Amplitude  of  v;  velocity  of  flight 

M = Lu  x y ZJT 

Element  of  coefficient  matrix  of  second  order  variables 
Second  order  variable;  forward  Cartesian  coordinate 
Element  of  coefficient  matrix  of  first  order  variables 
First  order  variable 

Element  of  coefficient  matrix  of  zero  order  variables 
Zero  order  variable 

Eigenvalue,  i.  e.  , coefficient  in  Equation  (2-5)  to 
describe  the  time  dependence  of  transient  motion. 

Note  Re(\)  = 0 indicates  neutral  stability. 

Shift  value  of  y in  eigenvalue  problem 

Shaft  rotation  of  control  surface 

Ratio  of  viscous  damping  coefficient  to  the  critical 
viscous  damping  coefficient 

Eigenvalue  of  Equation  (2-8) 

Shift  value  of  \ in  eigenvalue  problem 

Angular  frequency 
Denotes  complex  amplitude 


Actuator;  aerodynamic 
Flexible  body  motion 
Gyro 

Hinge  line 


l 


j 


E 

i * 


4 


, i 

2 


t 

I 

t 

i 
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p 

R 


x 

* 

X 

•• 

X 


y 

y 

z 

6 


Potentiometer 
Rigid  body  motion 

Corresponds  to  magnitude  of  second  order  variable 
Corresponds  to  velocity  of  second  order  variable 
Corresponds  to  acceleration  of  second  order  variable 
Corresponds  to  magnitude  of  first  order  variable 
Corresponds  to  velocity  of  first  order  variable 
Corresponds  to  magnitude  of  zero  order  variable 
Control  surface  rotation 


Matrix  Notation 


( 1 
( IT 

( rl 
I I 

L J 

r j 


Square  or  rectangular 

Transpose 

Inverse 

Column 

Row 

Diagonal 
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2.2  THEORETICAL  DERIVATION 


All  components  of  an  aero-servo-elastlc  system  can  be  regarded 
as  being  composed  of  elements  whose  characteristics  are  described  by 
second  order  (or  lower)  differential  equations.  Let  |x|  denote  the  set  of 
second  order  variables,  jyj  the  set  of  first  order  variables,  and  jsj  the 
set  of  zero  order  variables.  The  aero-servo-elastlc  equations  of  motion 
can  be  written  In  the  general  form 

raw + raw  + pyw 

♦ + Dgiw  + [zjw  * m <2'1> 


where  |F|  denotes  a forcing  function.  The  mathematical  formulation 
of  the  stability  problem  requires  first  order  differential  equations.  Wc 


By  setting  the  forcing  function  to  zero,  and  letting 

M = M exP  WO  (2-5) 

we  obtain  the  eigenvalue  formulation  of  the  stability  problem 

(y[A|  + [BJ ) {V}  = 0 (2-6) 

Instability  occurs  when  the  velocity  and/or  system  gains  are  such  that 
the  real  part  of  y is  positive. 

Since  [A]  is  singular  and  [B]  may  be  singular,  we  let 

Y * Y0  - i/U-V  tt-7) 

where  y and  X are  arbitrarily  chosen  complex  numbers.  Then 
Equation  (2-6)  can  be  rewritten  in  the  canonical  form  of  the  eigenvalue 
problem  as 

MVI  * <Y0[A]  +[B]f ' |A]  + \a  (vo[A]  +[B])]{V[  (2-8) 

where  the  new  eigenvalue  is 

x = Xo  - l/( v - Y0)  tf-9) 

The  eigenvalue  extraction  leads  to  convergence  to  the  eigenvalues  X 


corresponding  to  the  value  of  y closest  to  y . The  non-Hermitian 
matrices  yield  roots  that  appear  separately,  as  close  pairs,  or  as 
complex  conjugates. 

The  frequency  and  damping  for  each  mode  are  found  from  the  real 
and  imaginary  parts  of  each  eigenvalue.  Letting  a denote  the  decay 
rate  coefficient,  p the  damped  frequency,  ^ twe  viscous  damping 
coefficient,  and  w the  undamped  frequency  of  each  mode,  we  may  write 


y a o + ip 

(2- 10a) 

= uk+iA.  ^ J 

(2- 10b) 

from  which  we  find 

w a A2  ♦ p2 

(2-11) 
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and 


l ■ o/w 


(2-12) 


The  equations  of  motion  of  the  aeroelastlc  system  including  the 
shaft  rotation  of  the  control  surface  (flipper)  appear  as 

[M](h|  +[C]|ti|  * [K](|h|  - [h6]W)  = {F|  (Z-13) 

where  |h|  denotes  the  displacement  of  the  control  point  masses  In  the 

lumped- parameter  aeroelastlc  system,  [M]  is  the  mass  matrix,  [C] 

is  the  viscous  damping  matrix,  [K]  is  the  stiffness  matrix,  [h^]  is  the 

flipper  displacement  matrix  (the  elements  of  [he]  are  zero  for  points 

th  u 

off  the  flipper  and  equal  (x^  - x^)  for  the  i—  control  point  where  x^  is 
the  hinge  line  coordinate),  and  |F|  is  the  external  force  matrix.  For  the 
purposes  of  closed-loop  stability  analysis,  the  external  force  of  interest 
is  the  aerodynamic  force  induced  by  the  motion.  A survey  of  unsteady 
aerodynamic  influence  coefficients  (AICs)  has  been  given  in  Reference  3. 
For  present  purposes  it  is  sufficient  to  write  the  aerodynamic  force  as 

- [S.]  (4  <2-|4> 

where  [Q  ] is  a complex  matrix  of  oscillatory  AICs  valid  only  for  harmonic 
motion.  The  AICs  [(5.1  are  dependent  on  the  planform,  the  altitude,  the 
flight  Mach  number,  and  the  reduced  frequency  (Strouhal  number)  of  the 
motion  k = ubr/V  where  to  is  the  frequency,  is  a reference  dimension, 
and  V is  the  velocity  of  flight.  This  limitation  to  harmonic  motion  re- 
flects the  state-of-the-art  of  unsteady  aerodynamic  theory,  viz.  , consider- 
ably more  solutions  have  been  found  for  harmonic  motion  than  for  arbitrary 
transient  motion.  For  this  reason  flutter  analysis  has  traditionally  re- 
quired a trial-and- error  solution  to  find  the  velocity  and  frequency  for 
neutrally  stable  oscillations.  To  the  same  extent  the  aero-servo-elastic 
stability  analysis  must  be  carried  out  by  trial  and  error.  Equations  (2-13) 
and  (2-14)  may  be  combined  to  appear  as 


t [C]  |h|  + [K](jh|  - [hj  |M)  = 0 
whore  the  mass  matrix  now  includes  the  aerodynamics. 

m * csj 


(2-15) 


(2-16) 


The  inputs  to  the  hydraulic  actuators  come  from  rate  gyro  and 
accelerometer  feedback  loops.  The  angular  velocity  of  a rate  gyro 
may  be  found  by  numerical  differentiation  of  the  displacement  velocities. 

I ft*  = lDJ  {h}  (2-17) 

where  [Dj  is  a differentiation  row  matrix.  Methods  of  numerical 
differentiation  have  been  discussed  thoroughly,  e.  g. , by  Milne*  If  the 
differentiation  is  carried  out  locally,  i.  e.  , "in-the-small",  then  the 
elements  of  [Dj  will  only  be  nonzero  for  control  points  surrounding  the 
gyro  location,  e.  g.  , parabolic  differentiation  would  involve  only  the 
three  control  points  closest  to  the  gyro.  The  acceleration  at  an 
accelerometer  is  found  by  numerical  interpolation  of  the  displacement 
accelerations. 

h’=lHJM  (2-18) 

where  (Hj  is  an  interpolation  row  matrix.  Methods  of  numerical  inter- 

4 

polation  have  also  been  discussed  by  Milne  . If  the  interpolation  is  also 
carried  out  locally,  the  elements  of  [HJ  will  also  only  be  nonzero  for 
control  points  surrounding  the  accelerometer  location. 

Since  the  aeroelastic  system  has  a large  number  of  degrees 
of  freedom,  it  is  desirable  to  reduce  them  by  a series  or  modal  method. 
We  write  the  deflections  in  terms  of  a series  of  vibration  modes,  rigid 
-body  modes,  and  the  shaft  rotation(s)  of  the  flipper(s). 

|h(  = [hF](aF[  * [hR]j»B}  + [h,]W  (2-19) 

where  [hp]*  is  a matrix  of  restrained  and/or  free-free  vibration  modes, 
japjis  the  corresponding  set  of  generalized  coordinates  of  the  vibration 
modes,  [hp]  is  the  matrix  of  rigid  body  deflection  modes,  and  |a^|  is  the 
set  of  amplitudes  of  the  rigid  body  motions.  Substituting  Equation  (2-19) 
into  Equation  (2-15)  leads  to  a truncation  error  in  the  solution  since  the 
series,  Equation  (2-19),  is  finite.  The  error  matrix  may  be  interpreted 
physically  as  a distributed  force  on  the  system.  If  we  impose  the  con- 
dition that  this  error  does  no  work  in  any  of  the  flexible  or  rigid  modes 
(N.  B.  , this  is  the  method  of  Galerkin)  then  the  generalized  equations 
of  motion  are  found  by  premultiplying  Equation  (2-15)  by  [hp]  and  setting 
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to  zero,  and  also  by  premultiplying  by  and  setting  to  zero.  Sub- 

stituting Equation  (2-19)  into  Equation  (2-15)  and  premultiplying  by  [hp] 


leads  to 


[mF3{aFl  +[mFR]{aRl  +C"F5]^ 

+ [cf3)AfI  +^feFjlaFf  = 0 


(2-20) 


where  we  have  noted  that  rigid  body  displacements  produce  no  internal 
damping  or  structural  forces,  and 


py c 

Ffr!  = [hF]TMM 

[*F6>  [hF]TM[h5] 

OJ  * 


fbFJ  = [Hf3  MD'fI 


(2-21) 


(2-22) 


(2-23) 


(2-24) 


(2-25) 


We  also  have  noted  that  the  generalized  stiffness  matrix  has  a diagonal 

form  by  virtue  of  the  orthogonality  of  vibration  modes,  whether  they  be 

$ 

restrained  modes  or  free-free  modes,  and  can  be  written  in  terms  of 
the  generalized  mass  matrix  and  frequency  matrix  as 


where 


i*fj  - rvJfvJ 

• K2  "FJ 


P"fJ  =CHf3  MOfI  ’ dia8onal  elements 


(2-26a) 

(2-26b) 


(2-27) 


This  is  also  the  case  for  arbitrarily  chosen  modes,  e.  g.  , modes  given 
by  polynomial  expressions,  if  their  associated  frequencies  can  be  defined. 
However,  inclusion  of  arbitrary  modes  leads  to  a convergence  requirement 
of  a large  number  of  terms  in  Equation  (2-19)  and  therefore  will' not  be  per- 
mltted  in  this  formulation. 
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and  Wp  is  the  frequency  of  the  vibration  mode  jhpj.  The  generalized 
damping  matrix  does  not  have  a diagonal  form,  as  does  the  generalized 
stiffness,  but  may  be  assumed  so  as  an  approximate  means  of  including 
structural  damping.  If  we  represent  the  structural  damping  in  its 
equivalent  viscous  form  we  may  write 


fcri  * I*F8F^"fJ 


(2-28a) 


= fgFwFmFJ 


i 


(Z-28b) 


where  gp  is  the  structural  damping  coefficient  (N.  B.  , gp  = 2 £p  where  C 
is  the  ratio  of  viscous  damping  coefficient  to  the  critical  viscous  damp- 
ing coefficient)  corresponding  to  the  vibration  mode  jhp|. 

Substituting  Equation  (2-19)  into  Equation  (2-15)  and  premultiplying  by 
[hRf  leads  to 


where  we  again  have  noted  that  rig  id  body  dis  placements  produce  no 


internal  damping  or  structural  forces,  and 

ftp]  * [hR]T[nIhF] 

(2-30) 

C"r]  =[hR]T[M][hR] 

(2-31) 

»nd  [SRJ  -[l.R]T[ra][h6] 

(2-32) 

If  the  generalized  aerodynamic  forces  for  specific  modes  are 
considered,  we  may  rewrite  Equations  (2-21)-(23)  and  (30)-(32)  as 


[ mp]  = |mp]  + [Op] 

(2-33) 

[mFld  “ 0"fr]  + Pfr] 

(2-34) 

["f*]  = [mF6]  + PF6] 

(2-35) 
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where 


(FrWr]*Pr]  <2-^> 

CmRF]  =CWFR]T  +[QrfJ  (2-3?) 

^R6^  = t"W+^R6^  (2-38) 

[mFR]=[kjO  (2-39) 

[mFfi]  = C^F^  MW  (2-40) 

CWr]  = Or]  M0J  (2-41) 

and  C"rJ-WtMW  <2-42> 


Note  that^m^]  vanishes  if  free -free  modes  are  used  throughout.  Also, 

Or]  *s  t*ie  rigid  body  mass  matrix  consisting  of  total  mass,  static  un- 
balances about  the  coordinate  origin  (not  necessarily  the  centroid),  and 
moments  and  products  of  inertia  about  the  coordinate  origin.  The  generalized 
aerodynamic  forces  are  found  (e.  g.  , from  Volumes  I or  II  of  this  report)  for 
the  appropriate  vibration,  rigid  body,  and  flipper  modes,  including  their 
coupling. 

The  matrix  partitions  in  Equation  (2-3)  may  now  be  defined  as  follows: 


m n c s 


m 

'ey 

C"fr] 

CmF&] 

0 

A 

aF 

n 

PR] 

FrJ 

0 

*R 

lxj]  i*|  * 

feej 

> 
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0 
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(csij 
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ns 

(FSAJ 
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(RSAJ 

[SC] 

[SS2] 

e2 

V J 
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c o 

e? 

fc 

'V 
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0 0 0 0 

*R 

[X  J w * 

c 

i 

0 0 fccj  (CS2) 

> 

• 

6 

ns 

(FSGJ  (RSGj  [SC]  (SS2]_ 

*.e2.> 

m 

m n c s 

ffepj  0 0 0 

Ta 

aF 

n 

0 0 0 0 

aR 

[xj  1*1  ■ 

C 

< 

0 0 fccj  [CS2] 

> 

6 

ns 

0 0 [SC]  [SS2] 

(2-44) 


(2-45) 


where  e2  refers  to  the  set  of  second  order  servo  elements  and 


CC 

CS2 


FSA 

FSG 

RSA 

RSG 

SC 

SS2 


Control  surface  coefficients. 

Coefficients  for  second  order  servo  element  from  which 
there  is  input  to  a control  surface. 

Elements  of  the  matrix  -KRlHj[hp]  for  the  accelerometer 
Into  which  there  is  input  from  the  body  acceleration. 
Elements  of  the  matrix  -Kg[Dj[hp]  for  the  rate  gyro  into 
which  there  is  input  from  the  body  angular  rate. 
Elements  of  the  matrix  -K^Hflh^]  for  the  accelerometer 
into  which  there  is  input  from  the  body  acceleration. 
Elements  of  the  matrix  for  the  rate  gyro  into 

which  there  is  input  from  the  body  angular  rate. 
Coefficients  for  control  surface  from  which  there  is 
input  to  a servo  element. 

Coefficients  for  second  order  servo  element  from  which 
there  is  input  tP  another  servo  element. 
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The  number  of  rows  end  columns  in  each  partition  is  indicated 

m = Number  of  flexible  modes, 
n = Number  of  rigid  body  modes, 
c = Number  of  control  surfaces. 

ns  = Total  number  of  servo  elements  (zero,  first  and  second 
order),  ns  = s + q + r. 

s = Number  of  second  order  servo  elements. 

<J  = Number  of  first  order  servo  elements, 
r * Number  of  zero  order  servo  elements. 


CY*)M  = 


m 0 
n 0 
c fcS  | 
ns  jSSl] 


(2-46) 


ww 


c £si] 
ns  [SSl] 


(2-47) 


where  ei  refers  to  the  set  of  first  order  servo  elements  and 

CSl  = Coefficients  for  first  order  servo  element  from  which 
there  is  input  to  a control  surface. 

SSl  = Coefficients  for  first  order  servo  element  from  which 
there  is  input  to  another  servo  element. 


CZJH  = 


0 
0 

[cso] 

[sso] 

where  eO  refers  to  the  set  of  zero  order  servo  elements  and 


m 

n 

c 

ns 


{e0| 


(2-48) 


CSO  = Coefficients  for  zero  order  sirvo  element  from  which 
there  is  input  to  a control  surface. 

SSO  = Coefficients  for  zero  order  servo  element  from  which  there  is 
input  to  another  servo  element. 


A study  of  the  format  of  Equations  (2-43)-(48)  shows  that  each  of  the 
six  matrices  is  composed  of  partitions  that  may  be  described  as  aero- 
elastic,  servoelastic,  or*  simply  servo  terms.  PASES,  a compu.er  pro- 
gram designed  to  calculate  clcsed-loop  aero-servo-elastic  stability,  forms 
the  various  matrices  from  the  servo  partitions  which  are  input  directly  and 
the  aeroelastic  partitions  which  are  generated  from  basic  input  (e.  g.  , modes, 
frequencies,  mass,  and  aerodynamic  data). 
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APPENDIX  B 


AN  APPLICATION  OF  PROGRAM  MPASES 
TO  A 

TYPICAL  AIR-TO-AIR  MISSILE 

Introduction 

The  computer  program  MPASES  developed  in  this  report  used  the  same 
example  problem  to  demonstrate  its  usage  as  had  been  used  in  the  original 
version  of  the  program.  That  example  of  an  air-to-air  missile  was  not  a 
very  practical  one  to  the  extent  that  it  had  no  wing,  the  weight  of  the 
aft  end  of  the  fuselage  was  assumed  to  oscillate  with  the  control  surface 
(flipper),  and  there  was  only  a single  (rate-gyro)  feedback  loop.  In 
addition,  the  example  used  an  incorrect  transfer  function  for  the  actuator 
[Eq.  (78)  should  have  read  5/e3  - K /s(T  s + 1)]. 

The  new  example  studied  here  is  still  somewhat  idealized  but  is  more 
practical  to  the  extent  that  a wing,  with  its  attendant  aerodynamic  loads, 
is  added,  and  a light-weight  dynamically  balanced  flipper  with  bending  and 
torsion  degrees  of  freedom  Is  considered  separately  from  the  body.  A more 
realistic  servo  system  is  also  considered  with  three  feedback  loops,  a rate 
feedback  loop  (with  a rate-gyro  as  before),  and  attitude  and  acceleration 
feedback  loops. 

In  addition  to  being  a more  practical  problem  which  demonstrates  a 
more  general  use  of  Program  MPASES,  the  sample  also  illustrates  certain 
peripheral  aspects  of  the  problem.  One  is  the  calculation  of  a coupled 
mass  matrix  for  the  flipper.  Another  is  the  calculation  of  the  wing  and 
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flipper  aerodynamic  influence  coefficients  from  Piston  Theory  for  a more 
practical  range  of  supersonic  Mach  numbers  and  their  assembly  to  obtain 
the  proper  partitioned  format  for  the  complete  vehicle. 

The  servoelastic  stability  problem  is  considered  first,  and  then  the 
aerodynamic  influence  coefficients  are  added  next  to  consider  the  aero- 
servo-elastic  stability  analysis. 

Typical  Air-to-Air  Missile  Example 

The  typical  air-to-air  missile  is  idealized  as  shown  in  Fig.  3.  The 
ftiselage  is  idealized  as  a uniform  beam  with  ten  equally  spaced  masses. 

The  wing  is  assumed  to  be  rectangular  and  massless,  rigid  in  the  spanwise 
direction  and  to  bend  with  the  fuselage  in  the  chordwise  direction.  The 
flipper  is  also  rectangular  and  is  assumed  to  be  rigid  in  both  spanwise 
and  chordwise  directions  but  to  have  root  flexibility  in  bending  and  tor- 
sion. The  flipper  mass  distribution  is  simulated  by  eight  masses  connected 
by  rigid  massless  bars  to  the  four  degrees  of  freedom  of  the  flipper.  The 
coupled  mass  matrix  for  the  flipper  is  obtained  from  the  least  squares 
solution  in  Appendix  B of  Ref.  12.  The  various  input  matrices  to  MPASES 


are  derived  below, 


The  Mess  Matrix 


The  missile  fuselage  is  assumed  to  weigh  1000.0  lbs.  and 
each  flipper  weighs  5.0  lbs.  The  ten  equal  fuselage  weights  on 
the  half -fuselage  are  50.0  lbs.  apiece.  The  flipper  idealization 

i 

is  shown  in  Fig.  4.  The  flipper  center  of  gravity  is  assumed  to 
be  0.50  in.  forward  of  the  hinge  line  and  3.86  in.  outboard  of  the 
fuselage  side;  the  hinge  line  is  located  5.333  in.  aft  of  the  leading 
edge.  The  pitching  moment  of  inertia  about  the  hinge  line  is  50.0 
lb-in^  and  the  rolling  moment  of  inertia  about  the  fuselage  side  is 
137  lb-in^.  The  flipper  is  balanced  so  it  has  no  product  of  inertia 
with  respect  to  the  hinge  line  and  fuselage  side. 

The  input  to  Program  MASS  consists  of  the  inertial  data  and 
the  coordinates  of  the  eight  mass  points.  The  annotated  input  deck 
and  program  output  are  shown  below  in  Tables  1 and  2,  respectively. 

The  coupled  mass  matrix  for  the  four  flipper  degrees  of 


Table  1 - Input  Cards  for  Program  MASS 


C TITLE  ANO  SUBTITLE  CAROS 

LEAST  SQUARES  HASS  MATRIX 
..MASS  ELEMENTS  I* OH  TYPICAL  AIK-TO-AIR.MlSSILt 
C NUMBER  OF  STM  IPS 

l 

_C_  NUMBER  OF. MASSES. 

• 

C INERTIAL  PROPERTIES 

-.5.0  2.5 19.3 50.0 

C MASS  POINT  COORDINATES 

I. 333333333  3.0  -2.b6b666667  3.0 

_l. 333333333  6.0 I.333J33333  9.0 

•6.666066667  9.0  •2.666660667  9.0 


0.0. 


•6.666666667  3.0 
*^•666660667  64 0_ 


Table  2 - Output  from  Program  MASS 


CkAST  SUUXM6S  MAOS  MATrtlX 

NASA  616m6m7S  ►OR  TtrICai  AlK*TO-AH  MISSIL6 


MASS  MATRIX  l MRU  f UATA 
I SlMlRS 

MUMat*  OX  MASS6S  at  STAIR 


IU  II  13  13  1*  IS  Ik  IT  la  19  30 


STAIR  l 

torxL  mass  m ■ .suuoooout.oi 

X UMaALAMCk  MX  aAM  ■ .2*4000006*01 

* unHXLXMCt  mt  a.M  . . lAxuOogOt >u3 

V MOMfeNT  07  IMfcR I I A I sua  » ■ .*04904096*02 

X MOMENT  07  IMt«*IIA  I sua  X • . 13 7009446 • 0 3 

RMOUUCT  OF  Inert f lA  I sua  XT  • 0. 


A ( 

11 

■ 

. 1 JJJJJ3J6*41 

T1 

11 

• 

.340404406*01 

x< 

31 

• 

*.206666676*01 

T< 

21 

■ 

.349009996*01 

At 

3> 

■ 

••666066076*01 

T l 

31 

• 

.300040446*01 

XI 

Ai 

• 

.1  0 133396*91 

T< 

*1 

■ 

.649009446*41 

Xi 

SI 

■ 

.133334336*41 

T < 

>1 

• 

.409000446*91 

Ai 

*» 

• 

Vi 

61 

• 

.699004046*91 

Ai 

71 

• 

ri 

7) 

• 

.944004996*91 

Ai 

• 1 

■ 

*•266606676*91 

T 1 

• 1 

m 

.449090046*01 

* • » 

*N  i 

■ 

2 

BASIC  MASS  XMU  KtUUMUANl 

"ASS 

OtUMtlMT  maTMIC! 

.soaoooc*oi 

.3SOOU4t*9l 
« I9309OA*93 
_.>u004u6»0<_ 
• I 370046  *43 


. igoouoffli 

". 1333336*01 
• 3ooooi>t*0i 
. .1  77T7*»*Ol . 
.9900993*01 
.•904993*91 


.1900006*01 
*. 0660077 *01 
•9909496*01 
.••••••6*03 

.6100996*03 
■ ••OOOWOt  +0^ 


__  .1900996*01 
"-.200667**41 
•3494496*01 
..  .7111116*91 
•4990496*41 
*•0990946*41 

.1004996*01 

*•3946676*91 

.4004946*01 

.7111116*41 

•4190996*03 

*.3990466*93 


•1464466*01 
••6666676*41 
.9094946  *9 1 
..•••••6*03 
.6994046*01 
••2094096*03 


.1404996*01 
. 1333396*41 
.6494906*91 
. 1 7 7 7 706*41 
.9604996*02 
•6900946*41 


.11100096*01 
• 1333336*4 1 
.9000496*01 
.1777766*01 
.0140096*03 
.1200046*92 


.1000006*01 
*.00666  76*41 
.0004096*01 
.•♦•*•*6*02 
.3004446*42 
*.6000006*03 


MASS  MAIM1A 


STrir  l 

'.••24*7226*41 
.•16600616*42 
.104613096*41 
*.  169*72236 *41 
. 162S69**C  *4 1 
204*166 76*91 . 

•233413496*01 
*. 1 14*  73326*9 1 


As  discussed  in  Appendix  B of  Ref.  12,  the  negative  values  have  no 
particular  significance.  The  numerical  value  of  the  coupled  mass 
matrix  becomes 


4.433333  0.0010417  >0.396181 


[M]  = 


0. 970139 


symmetrical 


•0. 726042 


1.133333  >0.296181 

1.297917 


> 


Modal  Matrices  and  Frequencies 

The  modal  matrices  and  frequencies  for  the  fuselage  are  the 
same  as  those  used  in  the  MPASES  example  problems  with  the 
addition  of  the  deflections  of  the  four  flipper  mass  points.  These 
are  found  from  numerical  interpolation  and  differentiation  among 
the  last  four  fuselage  control  point  deflections,  i.  e.  , mass  points 
7,  8,  9,  and  10,  to  find  the  deflection  and  slope  at  the  flipper  hinge 
line  from  which  the  deflections  on  the  flipper  follow.  The 
Lagrangian  interpolation  formula  is 

(x-xg)(x-x9)(x-x10) 
h(x)  = h7  <x7-x8)(x7-x9)(x7-x10) 

(x-x7)(x-x9)(x-x10) 

+ h®  (*s  ' x? ) (x8  * fx8  - x J 0 ) 

(x-x7)(x-x8)(x-x10) 

+ h9  (=c,-*7)(*,-*8Kx,-*10) 

(x-x7)(x-xg)(x-x9) 

+ h10  ^i0-*7,(x10'x8)(x10’x9) 

The  hinge  line  is  at  x * 139.333  and  the  remaining  coordinates  are 
x7  = 97.5,  Xg  = 112.5,  x9  = 127.  5,  and  x^q  = 142.5,  so  the  interpolation 
formula  becomes 

hRL  « 0.049654  h?  - 0.  232235  hg  + 0.526618  h?  + 0.655963  h1Q 

By  differentiating  the  interpolation  formula,  we  obtain  the  fuselage 
slope  at  the  hinge  line,  and  the  result  is 


* 

i 


i 

t 
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h^L  = -0.009633745  h?  + 0. 048160494  hg 
-0.134086420  h?  + 0.095559671  h1Q 


From  the  hinge  line  deflection  end  elope  the  flipper  deflections  ere 


found  from 


“u  * »13  ■ hHL-  >•  "3333  1^ 


h12  * h14  * hHL+  666667  h^L 


This  generetes  the  four  flipper  deflections  required  to  complete  the 
description  of  the  three  fuselege  modes. 

The  fourth  end  fifth  vibretion  modes  ere  bending  end  torsion 
of  the  flipper  uncoupled  with  the  fuselege.  In  bending,  if  hjj  = hj4  = 

1.0,  then  hjj  = hJ2  a 0. 333333  since  the  flipper  is  eisumed  rigid, 
end  in  torsion,  if  hj2  = hj^  a 1. 0,  then  hjj  * h ^ = -0.  2,  es  cen  be 
seen  from  the  geometry  in  Fig.  4. 

The  control  surfece  deflection  mode  is  similer  to  the  flipper 
torsion  mode  except  thet  it  corresponds  to  e unit  control  surfece 
rotetion  (one  redien,  positive  treiling  edge  up),  so  thet 
1.333333  end  h12  = hu=  -6.666667. 

The  previously  essumed  frequencies  for  the  fuselege  were 

45.0,  125.4,  end  248.2  He.  The  flipper  frequencies  ere  essumed 
to  be  100.0  Hz  in  bending  end  115.0  Hz  in  torsion. 

The  damping  coefficients  in  the  three  fuselege  modes  ere 
essumed  to  be  g =0.01,  0.02,  end  0.03,  respectively,  end  g = 


0.03  is  elso  essumed  for  both  flipper  modes. 
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Servo  System 

The  servo  system  block  diagram  is  shown  in  Fig.  5.  The 
subscripts  on  the  outputs  are  the  same  as  the  servo  element 
numbers.  The  system  consists  of  three  loops  with  rate,  attitude, 
and  acceleration  feedbacks.  The  transfer  functions  and  their 
corresponding  differential  equations  and  numerical  values  are  dis- 
cussed in  the  following  sections  for  each  component  in  the  three 
loops. 


•w 


t 

! 


I 
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Control  Surface  Actuator  - The  actuator  servo  is  assumed  to  move 


the  control  surface  shaft  with  an  angular  rate  proportional  to  the 
difference  between  the  desired  and  actual  angular  displacements. 
Its  transfer  function  is 


„ _ 6 _ 1 
fs  - - s/R^+T 


and  the  corresponding  differential  equation  is 


(1/^)6  + 6 - e5  = 0 


The  numerical  value  1/Ky  =0.016  sec.  is  used  in  the  analysis. 


Servo  Element  No.  1 - The  rate  gyro  transfer  function  is 


K 


e./9  = g — ■ " ■■  ■ 

• %g  + (2Cg/«g)«  + 1 


and  the  differential  equation  is 
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The  numerical  value*  are  taken  as  l/w^  = 0.659  X 10*  sec  , 

2 f /u»  * 0.  359  x 10*  ^ sec,  and  K * 1.  0 deg/sec  per  deg/sec. 

8 8 8 

The  numerical  differentiation  matrix  to  obtain  the  rate  6 is 

baaed  on  a cubic  spline  fit  among  the  deflections  of  mass  points 

2,  3,  4,  and  S.  The  differentiation  matrix  is  (see  Ref.  LI) 


IDJ  * (1/244)1.0  +1  -27  +27  -1  0 0 0 0 0 0 0 0 0J 


where  6 = 15.0  in.,  the  distance  between  fuselage  mass  points. 

The  coefficient  1/244  =0.002777778  is  input  as  the  rate  gyro  gain, 
and  the  integer  elements  of  the  differentiation  matrix  are  input 
directly. 

Servo  Element  No.  2 - The  noise  filter  transfer  function  and 
differential  equation  are 

YRC  * e2//>el  = tTTT 

and 

t®2  + «2  * ®L  = 0 

where  t =0.001  sec. 


Servo  Element  No.  3 - For  the  variable  gain  amplifier  we  use 

e3/e2  = Kb 
or 

•3  * Kb‘2  * 0 
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where  Kg  = 0.1  deg/deg/sec  for  a Mach  number  of  3.0  at  sea  level. 
Although  the  aero- servo-elastic  stability  analysis  covers  a range 
of  supersonic  Mach  numbers,  a constant  value  of  Kg  is  used. 

Servo  Element  No.  4 - The  compensating  damping  amplifier  equations 
are 

Tds  + 1 

YDA  = *4^*3  = 27"  2 ~ , . ” 

s /wd  + (2Cd/wd)a  + 1 

and 

(1/u»d)®4  + {2Cd/wd)®4  + *4  ‘ Td*3  ‘ e3  = 0 

where  l/«d  = 0.  253  X 10”^  sec2,  2Cd/wd  = 0.  223  X 10~2  sec,  and 
Td  = 0.016  sec. 

Servo  Element  No.  5 - The  equation  for  the  rate  loop  feedback 
junction  is 


0>*** 


where  Kg  =2.0  deg/sec  per  deg/sec. 

Servo  Element  No.  7 - The  equation  for  the  attitude  loop  feedback 
junction  is 


e7  * 'll  ‘ % 


*7  + *6  * 'll  1 0 


Servo  Element  No.  8 - The  electronic  integrator  in  the  attitude  loop 
has  the  equation 

e8/.?  = 1/. 


.8  - e?  = 0 


Servo  Element  No.  9 - The  accelerometer  has  the  transfer  function 


s + (2Ca/«*>a)s  + l 


and  equation  of  motion 


(l/«a)e9+  (2C/u#)e9+  ®9  * am  = 0 


where  1/w^  = 1.013  X 10*^  sec^  and  2Ca/wa  = 0.6366  x 10  ^ sec. 

The  numerical  interpolation  matrix  to  obtain  the  acceleration 
*m  is  also  based  on  a cubic  spline  fit  among  the  deflections  of 
mass  points  2,  3,  4,  and  5.  The  interpolation  matrix  is  (see  Ref.  11) 


LHJ  = (1/16)L0  -1  +9  +9  -1  0 0 0 0 0 0 0 0 OJ 


( 


The  coefficient  1/16  = 0.06  25  is  input  as  the  accelerometer  gain, 
and  the  integer  elements  of  the  interpolation  matrix  are  input 
directly. 

Servo  Element  No.  10  - The  acceleration  loop  feedback  equation  is 

*10  = ac  ‘ *9 

which  reduces  to 

«10+  e9  = 0 

for  the  stability  problem. 

Servo  Element  No.  11  - The  acceleration  loop  gain  is 

'll/e10  * K* 
or 

ell  ’ Ka*10  = 0 

where  =4.0  deg/sec  per  g = 0.002170  sec /ft. 
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Servoelastic  Stability  Analysis 


The  input  data  cards  are  prepared  according  to  the  MPASES 
input  instructions  in  Sect.  III.  The  data  are  found  in  the  previous 
sections  of  this  appendix  with  the  exception  of  the  fuselage  vibration 
modes  which  are  found  in  the  original  MPASES  example  problem; 
the  vibration  modes  for  the  present  example  are  found  by  combining 
the  previous  fuselage  modes  with  the  additional  flipper  mass  point 
deflections  derived  above  in  this  report.  The  shift  eigenvalue 
V » +100.0  is  chosen  to  scale  [A]  to  be  of  the  same  order  of 

Q 

magnitude  as  [B]  and  with  the  same  sign.  The  input  data  deck, 
with  comments  added,  is  reproduced  in  Table  3. 

The  printed  output  is  reproduced  in  Table  4 except  that  only 
the  modes  and  deflections  corresponding  to  eigenvalues  of  practical 
interest  are  shown.  (The  modes  and  deflections  corresponding  to 
zero  or  infinite  eigenvalues  are  not  shown. ) 

The  stability  analysis  results  come  from  the  solution  of  a 
30th  order  eigenvalue  problem.  The  order  of  the  eigenvalue 


problem  is  found  from 


N a 2(m  + n + c)  + ns  + s 


where  N a order  of  the  eigenvalue  problem 
m a number  of  flexible  modes 
n a number  of  rigid  body  modes 


c « number  of  control  surfaces 


I 


; 

i 

: 

! 


I 


I 


I 


ns  * number  of  servo  elements 
s * number  of  second  order  servo  elements. 

In  this  example,  m = 5,  n * 2,  cal,  ns  = 11,  and  a a 3,  so  that 
N * 30.  The  stability  analysis  results  are  printed  for  22  modes 
since  8 of  the  eigenvalues  are  complex  conjugate  pairs.  Only 
modes  4-7,  and  16-22  are  of  practical  interest.  The  seven 
oscillatory  and  four  non-oscillatory  solutions  are  all  seen  to  be 
stable.  The  motions  involved  in  the  modes  are  summarized  as 
follows. 


Mode  No. 
4 

5 & 6 
7 

16 

17 

18 

• 19 

20 

21 

22 


Dominant  Motion 

Actuator 

Accelerometer  (critically  damped) 

Noise  filter 

Rate  gyro 

First  body  bending 

Compensating  damping  amplifier 

Flipper  bending 

Flipper  torsion 

Second  body  bending 

Third  body  bending 
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Table  3 - Input  Cards  for  Program  MPASES  for 
Servoelastic  Stability  Analysis 


C riTtF  CAN05 

rf**t c ol  Atw-rn-Arw  missile 

StWVOELASTIC  ANAL  fS  IS 
C C0UTM0L  CAHO 

2 Li 5 2 1 


C SHIFT  KlOtMVALUfc 
100.0 


JLL 


JLfi_ 


C 

c 

50,0 

HASS  MATRIX 

1ST  ROW 

0.0 

0.0 

0.0 

0.0 

0.0 

0 . o 

0.6 

oJ7^ 

o.o 

O.o 

0.0 

0.0 

0.0 

c 

2ND  Rum 

T97r 

5731! 

0,0 

o.o 

0,0 

0,0 

0,0 

0,0 

0.0 

0.0 

0,0 

0.0 

0.0 

c 

“3RD  RTJR 

50,0 

0.0 

0.00 

o.o 

0,0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

t 

"TTTR-  RQR 

50.0 

0.0 

0.0 

0.0 

0,0 

0,0 

0,0 

0.0 

0.0 

0.0 

0.0 

c 

”5TW~R0R 

50,0 

0.0 

0.0 

o.o 

0,0 

0.0 

0.0 

0.0 

0,0 

0.0 

c 

~5TFT~ROS 

50.0 

0,0 

0.0 

o.o 

0,0 

0,0 

0.0 

0.0 

0.0 

c 

"7TR“ROW 

50,0 

0.0 

0,0 

0.0 

0.0 

0.0 

0,0 

0.0 

£ 

SYH  RUN 

50,0 

0.0 

0,0 

0.0 

O.o 

0,0 

0.0 

C 

'9TR  RUR  “ 

50.0 

0.0 

0.0 

o.o 

O.o 

0.0 

C 

10TH  ROM 

50,0" 

0.8 

0,6 

0.0 

0,0 

C 

1 1 Th  ROM 

9.9333333  0.0010917 

-0,3961806 

0.0 

o.97ot3B9  o.o 
C 1 JTM  RQM 

T71 333532 — -67295TW5 

C 1STH  RQW 

1.2079166 
T 


•0.7260917 


XTJSTROf  SUPf ACT  bFFtEfVfiTN  wtlQE 

0.0  0.0  0.0  0.0 

0.0  0_#0 0.0 0.0 


~ 17333T333S3-6.WW866  7 
c oamring  coefficients,  g 


0,0  0,0 

1 ,333333333-6,666666667 


0.01  0.02  0.03 

T HQDi  C~ VITBITTHN  T PrSTJOTI  FT 

«5.0  125,«  2«8.2 


0.03 


0^03 


VIBRATION  MflOE  SMARTS 

isr"MO0F 


100,0 


115,0 


1.0  0,911195  -0.122920  -0.532567 

•0.532967  -0,122929  0.911195  1.0 

■njTS/l  620 171388111 

C i 2N0  N00E 

1.0  -0,150871  -0.9J  0892  -o',990216 


■0,756209  -0,756209 

0,821620  1,138881 


•0,919093  0,919093 
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1 - »3iLiaa4Baasga£hibtt 


1 


A irt  m sAh  -fc  I***-  - 


3a?*0315 0. *10M3  0.150671  M'.fl -0,630356  -1,361*54 

•0,6 10356  -1,381*54 

e iho  hooc 

-0,8031*5  0.645103  1,0  0.105*86  -0,*a78*tt  -0,9478*3 

0.105*86  1.0  0.645103  -0.8031*5  -0,35*3*8  -1,168304 

-0.35*3*8  -1.1683Q4 


■5 — 

4TH"W'0ftfc 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o'o 

.3333333333 

.3331333333 

or 

1.0 

c 

5TH  MQOE 

Oj.0 

0.0 

0.0 

0.0 

0.0 

0.0 

^to 

o.o 

0.0 

0.0 

-0,3. 

1.0 

-0.2 

1.0 

C 

RICIO  BQOY  MOOES 

~C 

■pcuirce“ffoo£ 

1.0 

1.0 

1.0 

r.o 

1.0 

1.0 

1.0 

1.0 

1,0 

1.0 

1,0 

1.0 

1.0 

1.0 

c 

PITCH  MUOE 

T.5 

22 - 

37.5 

52,5 

67.5 

83.5 

IT  275 

1<7,S 

14275 

IT673  ” 

'146,0  “ 

138.  < 

> 146,0 

C 

CONTROL  SURFACE 

OUTPUT 

TTO 

O70T6 

0.0 

* 

c 

l 

NUMBER  OF  INPUTS 

TO  CONTROL  surfacf. 

c 

CnFfTROl  SURF fCE“ 

TTJPuTS  FRQfl“  SEftvnS  ” 

5 

0 -1,0 

C 

OUTPUT  FROM  SERVO  ElEMFNT 

1 

^2  170 

0 .35* 

•2  0.650 

• 5 

c 

1 

NUMBER  OF  INPUTS 

TO  SERVO 

ELEMENT 

c 

t npui  FRoff-Brrwr" 

FircHRArr 

0 

-1  -2.7777778 

-3 

c 

differentiation 

MATRIX  FOR 

f rate  gyro 

0.0" 

1.0 

-2775 

27.0 

-1,0 

0,0 

0.0 

0.0 

0,0 

0.0 

0,0 

0.0 

0.0 

0,0 

c 

uUtPUTfrTTM  SERVO  ELcmCnt 

3 

1 1.0 

0,001 

c 

NUMBER  of  inputs 

TO  SERVO 

element 

1 

c 

inputs  to  servo 

ELEMENT 

l 

0 -1.0 

c OUTPirTTflUM  SERVO  ELEMENT 

3 

0 1,0 

C 

number  of  inputs 

TO  SERVO 

ELEMENT 

I 

c 

INPUTS  TO  SERVO 

ELEMENT 

2 

0 -o.t 

c 

Output f rOh servo  element 

4 

2 1.0 

0.223 

•2  0.353 

-5 

C 

NUMBER  OF  INPUTS 

TO  SERVO 

element 

1 

c 

inputs  to  servo 

ELEMENT 

3 

1 -1.0 

-0,016 

t OUT?  DTTR  OH  -SgaVO'ElTX  CTTT 

5 0 1.0 


c 

NUMBER 

OF  INPUTS  TO  SERVO  ELEMENT 

1 

e 

-r 

INPUTS 

TO  SERVO  ELEMENT 

f 

• 

0 

1.0 

—r 

o M .o  - 

c 

OUTPUT 

FROM  SERVO  ELEMENT 

6 

0 

1.0 

c 

NUMBER 

OF  INPUTS  TO  SERVO  ELEMENT 

c 

~T 

inputs 

TO  SERVO  ELEMENT 

i 

2 

0 

•2.0 

i 

c 

OUTPUT 

from  servo  element 

7 

0 

1.0 

T“ 

; 

c 

C 

INPUTS 

TO  SEHva  ELEMENT 

\\ 

1 

"nr 

0 

1.0 

j 

it 

0 

■1.0 

r; 

h 

c 

OUTPUT 

FROM  SERVO  ELEMENT 

r 

f 

O 1.0 

£ 

e 

t 

NUMBER 

OF  INPUTS  TO  SERVO  ELEMENT 

k 

V 

Twnr 

~rq  servo "eiemikt 

7 

0 

•1.0 

1 

c~ 

OUTPUT 

FROM  SERVO  tLEMEwT 

f 

2 

1.0  0.&2&*  -2  1.013  -6 

I ’ 

c 

1 

NUMBER 

OF  INPUTS  TO  SERVO  ELEMENT 

c 

& 

0 

-2  -0.062b 

J ! 

c 

Interpolation  matrix  fop  accf.lf  rune  ter 

1 

.0 

-1.0 9^0 -M 

0.0  _ 

0 

.0 

0.0  0.0  0.0  0.0 

0.0 

i 

0 

.0 

0.0 

U 

c 

OUTPUT 

FROM  SERVO  ELEMENT 

10 

0 

1.0 

c 

x 

number 

OF  INPUTS  TO  SERVO  ELEMENT 

| i 

1 : 

c 

INPUTS 

TO  SERVO  ELEMENT 

9 

0 

1.0 

V 

c 

1 , 

n 

0 

1.0 

1 

c 

1. 

NUMBER 

OF  INPUTS  TO  SERVO  ELEMENT 

c 

OUTPUT 

FROM  SERVO  ELEMENT 

10 

0 

-0.002170 
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Table  4 - Output  from  Prof-?m  MPASES  for  Servoelaatic  Stability  Analysis 
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Aerodynamic  Influence  Coefficients 


The  aerodynamic  influence  coefficients  (AIC's)  at  supersonic 
speeds  are  found  from  Piston  Theory  and  are  generated  by  Program 
PISTON  (Ref.  4).  Three  Mach  numbers  are  considered,  2.0,  3.0, 
and  4.0,  and  four  reduced  velocities  are  considered,  1000.0,  20.0, 
5.0,  and  2.0  based  on  the  wing  semichord  of  1.5  ft.  The  AIC's  for 
both  the  wing  and  the  tail  (flipper)  are  found  for  the  same  reduced 
velocities  by  using  the  wing  semichord  as  the  reference  for  both 
surfaces. 

The  wing  planform  is  shown  in  Fig.  3 and  its  motion  is  deter- 
mined by  the  deflections  of  mass  points  Nos.  7,  8,  and  9 which  are 
located  at  the  fractional  wing  chord  positions  0.152778,  0.569444, 
and  0.986111,  respectively.  One  strip  is  sufficient  to  represent  the 
wing  since  it  has  been  assumed  rigid  in  the  spanwise  direction. 

The  basic  Piston  Theory  is  used  without  thickness  (a  small  thickness 
ratio  t s 0.001  is  input).  The  input  data  cards  are  listed  below 
with  comments  and  follow  the  input  format  described  in  the  Program 
PISTON  User's  Manual  (Ref.  4).  The  input  data  are  also  echoed  in 
the  program  output  which  follows  the  input  data  card  listing  below. 

At  the  end  of  the  printed  output  the  punched  card  output  is  listed. 

The  punched  AIC's  are  input  appropriately  according  to  the  proper 
partitioning  format  into  Program  MPASES. 

The  tail  (flipper)  planform  is  also  shown  in  Fig.  4,  and  its 
motion  is  determined  by  the  deflections  of  its  four  mass  points, 

Nos.  11,  12,  13,  and  14  which  are  located  at  the  tail  quarter-  or 

163 


three-quarter  chord  positions.  The  tail  has  been  assumed  rigid  is 
both  spanwise  and  chordwise  directions  but  to  bend  and  twist  about 
the  actuator  at  the  root.  Two  equal  width  strips  are  used  to  account 
for  die  bending  motion.  Except  for  the  chordwise  and  spanwise 
locations  of  die  aerodynamic  control  points,  the  tail  is  treated  the 
same  as  the  wing  in  terms  of  input  to  Program  PISTON.  The  input 
data  cards  are  listed  below  those  for  the  wing  in  the  following  pages 
and  the  printed  output  of  AlC's  for  the  tail  follows  that  for  the  wing. 
Again,  at  the  end  of  the  printed  output  the  punched  AIC's  are  listed. 


f Table  5 - Input  Cards  for  Program  PISTON 

u . 


I 


C TITLE  ANO  SUBTITLE  CAROS 
RISTON  THEORY  ACROOYNAMICS  FOR  WINS 

NfcN  ANO  IHPROVEO  MAC  .MISSILE 

CONTROL  CAROS 

12  11 

.3 0 A A ■ A 


RLANFOHM  GEOMETRY 
l.S 

SIR  1R_6CQMETR  Y_ 
l.S 


1.0 


3.0 


3.0 


0.1S277777S  0.569AAAAA*  0.966111111  O.A 

THICKNESS  DATA 

0.001  0.0_ 0.0 

C SERIES  uF  MACH  NUMBERS 
2.0  3.0  A.O 


.15 i.m_anglC_of_a_tt ack.for_.e ach  mach_number._ 


0.0 

c series 

OF  REOUCEO 

VELOCITIES 

FOR  EACH  MACH  NUMBER 

1000.0 

20.0 

S.O 

2.0 

1000. 0 

20.0 

S.O 

2.0 

iooo. o 

20.0 

S.O 

2.0 

C TITLE 

ANO  SUBTITLE 

CAROS 

_PISTON  THEORY  AEROOYNAHICS  FOR  TAIL 
NCm  HAC  MISSILE 


COnTmol  CAROS 

.0 2 l—JL 

3 0 A A 


PLANFORM  geometry 
l.S 


JL»bL 


JWL 


o.s 


o.so 

JbJUL 


STRIP  GEOMETRY 

0.66666667  0.25 

0 . 66666667 ft*  M. 

C THICKNESS  OATA 

0.001  0.0  0.0 

£ SEM1ES_0F  MACH  NUMBERS 

2.0  3.0  A.O 

C TRIM  ANGLE  OF  ATTACK  FOR  EACH  MACH  NUMBER 


-2r0_ 


0.75 

Jdl. 


0.0 

0.0 


-Vo 


SEMI |S  OF  REOUCEO  VELOCITIES  FOR  EACH  MACH  NUMBER 

000.0  ~20.0  S.O  2.0 

1000.0  20.0  S.O  2.0 

1000.0  20.0 hi hi 


I 


rmoN  TMfiMv  uiotfNMici  ro»  king 

NCM  AMO  IMMOVCD  MAC  MIlslLf 


MKIP  MO,  ALPHA  ItHO  (DURUM 


itanovMANfC  Influence  eoirrtctiNTS  it  mi  stun  theory  »itm  camber 

oscillatory  CASE 

macm  no,  • <(ooo««t 
t/MOl  • ,1000001*0* 

I STRIRS  - 

c*t  it  me  » i ■»  j 

, IS12O«*0E*07  •,?<*))«»(»«)  •, 1 R 1 20021C*07  ,«OW»Jll*0*  ,*001*8101*01 

~,Mn«l*5t«<0 ,0i]10M|(*42 ,Tt««3401E*0O  •,*3*7o25IC*01  •,770S*««OE*0*  -,0*721**11*02 

•,UM)UII«0O  ,3*01«tl*E*01  ,*J0*2S08I«0*  -, 0*721**11*0  2 *,)lM404)(*Ot  *,»IIM77t*0> 


aerooynanic  influence  coefficients  s»  fiston  tmccot  «i  th  emu 
oscillatory  case 

. . NaCN  no,  • 2.000000 

l/«CR>  * ,2000001*02 

______  - , ,Tl,tM 

c*(  t)  stir  » j or  i 

,«1M700C.OI  _*l|«Olorl8t*0J  _;,7110*0Ml*«J  ,112070001*01  , l****2*«t*01  ,U7O)«*ot«00 

, IO0*37*OC*O1  ,tl2O7000C»0t  ,lo070100t*0l  • , 1 * 1 0S250C *02  • , 11071*7*1*01  1 7**«7**l *0 1 

•,**|Ot702E*92  , 117OJ***r»00  ,201*70011*01  •,17*0*7*01*01  •,207*00200*01  •,0121**001*01 

AtROOVNANfC  INFLUENCE  COEFF ICtCNTS  07  FISTON  THEORY  * t TH  CAMBER 
OSCILLATORY  CASE 
MACH  N(l,  * 2,000000 

l/MR)  • ,0000001*01 

...  x #TP,PS 

CH(  1)  SHE  • 1 SY  1 

,120|O| t«E»02  •,17*0*«01E*01  •,*03100021*02  ,1111*7211*00  , 12««0*11C»02  ,2**0**09t«0l 

,tO0O«O74E*42 ,1111*7211*00 ,r)kl3BI0C*0l  •,«7«oll23Ct0l  •, 10*212376*02  •,*aSOl«*»C*00 

•,27000 l tot *01  ,2**0**0*E-01  ,107211272*02  •,*«0*l**00*00  -, 12*o2SloC*02  -, I200****E«0I 

AEROOYNAmjC  INRLUENCE  COEFFICIENTS  IT  RISTON  THEORY  NlTN  CAMBER 

OSCILLATORY  CASE 

RACN  NO,  » 2,000000 

t/K(R)  * ,2000001*01 

“ t STRtRS 

CN(  1 ) SIZE  * 1 SY  1 

,32O0S7«0(*0t  •,t«*00700E*41  •, 7RO0*0O«C*41  ,112*700*1*00  ,t*««l20*E»9t  , 1 1 7ol**«E«01 

;lO0*S7«tt*0t ;H2078*(C*00 ,1»*7S1S0C*01  -.101002301*01  -,I1071*70E*01  •,17**«7*0r«00 

*,MIO«702E*00  , 117Ol«O«E*0l  ,2SlS700K*«t  -.17«**7«0k»00  •,207*00201*01  •,Sl21**SSC*00 


\ 


k 


1 
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r 


t 


[ 

t 


■ 


[ 


AOROOVNAmtC  INflUCNCC  COl'MCIONT*  It  RISTON  T*CO*Y  «UTN  CAlMllR 

■ oscillatory  cut 

NACM  *0«  ■ 1,000000 

i/fttm  ■ ,ioooooi*oo 

" i imri 

ent  i)  mi  » » by  i 


,i7tr«mi»o* 

~,iim»i(«ot 

*(minnuot 


•.lottmit'O] 

.mioitwc*** 

,309*31*00*01 


.,tPto*ojie*o7 

,itl«IMH«0t 


,*•10)7120*02 

-,90*M0**e*02 


,1110*3130*0* 

-,3173013*1*0* 

•,]<)0i01ll»>« 


*,33o330*o£*42 

-,17092ORn*03 


aorooynanic  imiucncc  coomicionts  by  riston  thcory  »ith  camor 

OSCILLATORY  CASC 

NACM  no,  • 1,00400* 

• ,2000000*02 

t #Ttln 

CMC  J)  OIK  


,390714120*4)  «,100UOSlO*0*  -,«*«271|%e*0)  ,*00*T«2*C*4* 

— *,tomooi(*oj — ,oa407o2oe*oo  “-,  127001000*02 

•,2011*0100*02  ,770002001-01  ,107101020*01  •,110111010*01 


,1111101*0*01  ,770002000-01 
>,20701OOlC*01  -, 11*111010*01 
-,110000000*03  -.1*1091210*01 


ACROOYNAhiC  INOIUONCO  COtOMCIONTS  OY  RISTON  Th00*7  RlTH  CA-RlR 

OSCILLATORY  CASC 

NACH  NO,  « 1,000000 

l/KIR)  ■ ,3000000*01 

, 

cmc  t)  mo  • i sy  i 


,210100*90,02  -,290270270*01  • 

-,**0*0770001 -,2213109*0*00 

•,101*O*l|C,0t  ,102713790-01 


1*2*70700*02  ,22191*900*00 
09*2012*0*01  -,310010210*41 
10**01900*42  -,200970910*00 


,01*711)70*01  ,1*2713730-01 
-,120107000*02'  -,20*270311*00 
-,0*29S09ot*0t  -,132*11110*00 


ACROOVNANIC  INRLUONCO  COOfriCIONTS  SY  RISTON  TNOORY  «ITM  CARMR 
“ OSCILLATORY  CASC 


NACM  NO,  • 1,40*0*0 

1/RCRJ  * ,1004000*01 

1 STRIPS 


CM(  1)  SIZO  • I SY  1 


,190710120*41  10011*910*01  . 

“,107190010*01  ,*0*0712*0-01 

•,2019*0)00  *40  ,77***2000-02 


•,*|02712*0*0l . 
,00*0*9000 *00 
,|*7)*)020*0l 


,**OOT*2U-O| 
-, 127*91*00*41 
-,114111*10*** 


,1)399*1*0 *41  ,770*0*000-02 
-,2070100)0*41  -,110)11*10*00 
-,110*40400*01  -,1*1191730*** 
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• • » ‘ V £ * .f  ,'^**< 

^yiiiaiui 


USl 


tCftOOVNAMtC  INPLUINCI  COtMICIiNTi  If  PISTON  TNtQRv  *ITN  C**B1R 

o*eiu*ro«*  cm 

»»C»  NQ,  a a, ****** 

l/MR)  ■ ,IIMM»I« 

. _ ^ aTa,^a 


cm  t>  atit  • i it  i 


,»iimru«u  •,)rtii)ii(«i) 

~2oo*io*ti*o* — , jinn  turn 
•,laiumi»ll  ,2MIHlll»«l 


•,M«ri7iu*H  ,mt*usi*o2  ,Kusim»i«  .momiifot 

,uniuu«i*  •l«ri»]in«i]  ' *, 1*79*1121,0*  •liawi)iut«i 

,l|IIM|Hfl»  sHliallUMi  •lUr»IH)tf» 


*CRQOTN*Nf C IWLUCMCt  CSCMlCIMIt  If  PISTON  ThIOHT  at TH  C*"#t* 

oactuATO«»  cm 

aaen  no,  a a ,000000 

WMD  • .lotltltfli 

t strip*  ' ' * ' 

chc  it  am  a J BY  i 


.tiiniTiiti]  .,7iiiitm.n 

-,aitnrirt»«i — ,MiM)iit««i 

•.liUMMMi  ,i7oo*9i*c.oi 


•,nioooa*t*oi 

-,?|I1UIII*H  •,*1010*011*01 

•.mmirtfii 


,194)101*1*0)  ,iT0**9l»t-01 
• ,l»I»))(t«)  •,Ht1l|l7C«ll 
•,101)122)1*0)  •,211)71)01*01 


ichootnintc  inpluinci  cocrriciiNTs  at  piston  thiort  at tn  ctmci 
osctuaWT  cm 

. N*C*  NO,  a a, 900000 

W«(R)  a ,1000001*01 

“ t strip* 

CH<  1)  SIZl  > 1 If  ) 


_,t*«7t9|9C«02_>,t0A07t*0t*01 •, 227*11701*02  ,1**170171*00  ,*27121001*01  ,10212*791*01 

, 10 1 11* 1 7C  *0 1 , 1 6*17*1 7 1 *00  .4**1721)1*01  a,239*2*llC*0t  *,**990*101*01  *, 221  l**72C*00 

>,1)71717BC«01  , l*212*79(«0l  ,702071201*01  .,2211**72£*00  *,*«1701*lt»0t  *,*1**102*1*00 


oconoTNONtc  iNPtutNCi  coeinctiNts  ov  piston  tniont  *ith  camhir 
• oactu*to«T  cm 

N*CN  NO,  * 0,000000 

!/*(*>  a ,2000001*01 

l STRIPS 

CN(  l)  SXZ1  a 1 ST  1 


,2*1110701*01  *,71220*001*00 
',002*17071*00  ~,***2R229C.or 
•,219*7*011,00  ,1700991*1*02 


*,1030900*1*01  ,***202291*01 

“•,7*0101*11*00  *, *1*10*011*00 

,121210711*01  *,0021***71*01 


,1001*01*1*01  ,1700091*1*02 
•,1111011)1*01  *,**210**71*01 
*,10)112211*01  *,211)71101*00 
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i ; 

: j 


1,000006*03 

J > t 

3 

1,312046*00' 
— T,  023*56 *05 
•1,104246*03 
2,000006*01 

- •3  1 I 

3 

3,230306*02 
—1,000546  *02 
•4,410006*01 
3,000006*00 

3 T—  1 

3 

3,201016*01' 

—l,003OOt*0l 

•2,7*0016*00 

2,000006*00 

- Sit 
3 

3,230306*00" 

- 1,000306*00 
•4,aioOoe«0l 

1.000006*03 

- 311 
3 

0,707006*05" 
-2, •70756*05 
•7,330736*04 
2,000006*01 
j , 

3 

1,307106*02- 

- 1,071506*02 
•2,033406*01 

3,000006*00 

3 - t - t 

3 

2,101046*01" 

- 4,40*606*00 
•1,034006*00 

2,000006*00 

- 3 1 1 

3 

3,307106*00. 

- 1,071306*00 
•2,033406-01 

1,000006*01 

- 1 1 1 

1 

0,506776*03. 
-2,004146*03 
•3,404676*04 
2,000006*01 
— i-  i r 

i 

2, •33316*02. 
-0,024306*01' 
•2,104756*01 
5,000006*00 

3 l-  l 

3 

1,447106*01. 

3^015346*00 

•1,371726*00 
2,000006*00 
1 1 1 

3 

2,433316*00. 
4, 024546 -01 
•2,10473C»01 


•7,403346*02 

4,423006*01 

5,401006*00 


•1,412406*04 

3,744356*03 

4,200256*05 


•'.423046*01 

•0.502436*02 

•8.072406*01 


4,000576*03 

•7,744406*03 

•3,145016*05 


1,406006*01 

1,324706*00 

l,|7fe406*01 


•7.250016*02 

1,407426*02' 

2.313706*02' 


' 1 , 744706*00« 

3,111076-01 

2.049046-02 


4,331516*01 

0,341346*00 

1.572316*01 


1.408406*00" 
1,324706-01 
1. 174406-02 


5,005336*02" 

4,030376*01 

3,034116*00 


1.001116*01" 
0, 040746-01 
7.708416-02 


2.302746*00" 

2.215106-01 

1,027146-02 


1,001116*00" 

4,400746-02 

7,706016-01 


3.741416*02" 

3,311416*01 

2,654346*00 

7,522046*00" 

4,042026-01 

5,700006-02 


1,440726*00" 

1,045716-Ot. 

1,423236-02 


7,322046-01" 

4,442026-02 

3,700046-03 


'7,230416*00 

1,107826*40 

2,315706*00 


1,210446*04 

2,407146*03 

4,104086*03 


4,642736*02 

0. 0604*6*01 

1,  *73*36*02 


3,024716*01 
4,202016*00 
1 ,044026*01' 


4.042736*00 

0,008006-01 

1,473*36*00 


'O,0O727e*03 

1,873406*03 

1,131406*03 


3.438011*02 

7,401056*01' 

1,252406*02 


2,2741je*0l 

4,001726*44' 

7,028736*00' 


3,034016*00 

7,»OJ*36-0t 

1,232406*00 


ll 324706*00 

■i;oi85l6*41" 

'1.704446*00" 


S‘.ll  1076-01 
•4.70*316*00" 
••.48*206 -0l« 


r.S2«7*e-01 
•1,014336  *00- 
• 1 .704406-01- 


*‘.430376*01 
•0.302*16*02" 
"3. 0*5306  *01" 


8.0*07*6-0 1 
■1,278346*01" 
•1.101126*00" 


2’.2l5tO6-01 
.V  10*3*6*00" 
>2.042801-01" 


1.000636*02 

3,107406*02 

2,074006*02 


1,240006*01 
l,042t?E*OI 
1 ,20*236*01 


1.000836*00 

3,107406*00 

2,074006,00 


3,334036*05 

3,173026.05 

3.450206*05 


1.335346*02 

2,070376*02 

1,300086*02 


4,347136  *00 
1,203086*01 
0, *23516*00 


0*  0*07*6-02 
•l, 2705*6*00" 
•1. 103126-01" 


]'.  331*16*01 
>0.702336  *02- 
• •.4*2016*01" 


*’,**2426-01 

■ 0 543046*00* 

■ 0. 023476-01" 


1, **37 16-0 l 
•2,30*271*00" 
•2.231476*01" 


*,*•2026-02 

■0,343066-01" 

•6,023476-42" 


1,3333*6*00 

2,070176*00 

1,140086*00 


2,304506*03 

3,070*36*05 

2,342016*03 


t, 003406*02 
1,551036*02 
1,033126*42 


0,2712*6*00 
0, *OO0*t., oo 
*,437011  .00 


1,003406*00 
1,551656.00 
1, 011126*00 


074 

074 

074 

5,081086*00074 

-4,072*06*01074 

•2,542006*02074 

074 

074 

074 

1 .174406 - 01074 
•1.704*66*00074 
•5,12*006*00070 

070 

074 

074 

2,040006-02074 
-4 ,*8420C-0 10TM 
•1 .281006*00070 
074 
074 
074 

1.174406- 02070 
-1,700486-01074 
-3,124006-01074 

074 

P74 

074 

3,83*316*00074 

-5,0*5506*01074 

-1,703276*02074 

074 

074 

074 

7.708*36-02074 

-1,103126*00074 

-3,410516*00074 

P7M 

074 

074 

1,0271*6-02074 
-2,082806-01074 
-0,324136-01074 
074 
074 
074  i 

7,708*36-01074  i 
•1,103126-01074  i 
-3,410536-01074  i 
074 
074 
074 

2, 0503 Ot *00074  : 
•4,4*2036*01074 
•1,27*006*02074  ! 

074  ; 
074  i 
074  i 

3,700006-02074  j 
•4.025076-01074' i 
•2,533736*00074  i 

074  : 
074  : 

074  . 

1,425251-02074  : 
•2.231*76-01014  : 
••,30*106-01074  ; 

074  I 
074  < 
074  I 

3,7OO«o6-O3074  . 
•0,023876-02074  < 
•2,333736-01074  i 
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*2»OOYNAFIC  (NFIUCNCI  COtFFtCItNTS  NT  Ft  STUN  7M200Y  •IThOuT  CiNU 

OSCtUAttlOY  case 


MACH  NO,  * 


2,000009 


t/K(H)  a , 1 000001*0* 

2 9701FS 

CMt  t)  S I Z»  a 2 *Y  2 


,101*423)2*4*  •,200043002*0)  •,401002)31*00  , )7  | Ofl**tt  .02 

— ;*4M0109f  *0* ;)7l 00**12.02 •,  44(401 002*00  •,240020231*03 

CM(  2)  lilt  a 2 OT  2 


,3010023)2*00  •,20004)002*03  •,4010023)2*00  , )T l 0409)2 *02 

— ,490491002*40 ;3Tl0004jt*02  “ *,444401002*00  •,234020232*03 


AERODYNAMIC  tNFUUCNCl  COIF* ICICnTs  by  FtSYON  TM20HY  >1 ThOuY  CAFBlR 
- — OSCILLATORY  CAS2 


HACK  NO,  a 


2,000000 


1/4(01  ■ ,2000002*02 

2 STRIF8 

CMC  1)  St Zfe  a 2 OY  2 


,2000404)1,01  _•, 32000*122. 01  2400404 12*03  ,7*2123012*00 

,144)00402*0)  ,74217)412*00  •,144)00402*0)  »,3 1 72*0««t »0l 

.CMC  2)  9(21  * 2 OY  2 

_j200O»04)2,0)  •,320040122*01  •,2044404)2*0)  ,70217)412*00 

,144)00402*0)  .702173412*00  '•,144)00402.0)  », 317240*42*01 


*2«ooym*mic  influence  coefficients  or  Ft sro*  tmoomy  «i7now7  ca »*2» 


OSClLLAlURT  CAS* 


MACM  NO,  • 


2,000090 


t/KCR)  > ,3000002*01 

2 S7MIFS 

CMC  I)  Size  a 2 BY  2 


,123490302*01  •,1)00213)2*01  •, 123«00302, 02  , l03S4)a32*00 

— ^12*090232*02 ,14334)432*00 4,1**000232*02  * •, 124)| 0121*01 

_ CMC  2)  3IZ2  ■ 2 OY  2 

,123*94301*42  •, 1)00213)2*01  •,123*00302.02  ,10334)432*00 

— ,124*00232*02  — ',1*334)432*00 •, 1**000*31*02  ' •,1*4)10122*01 


OtOOOVNAMIC  JNFLU4NC*  CU2FF IC 1CNTS  OY  FIS70M  7M20NY  aI7»Ou7  CANR2R 

uicruATooY  cose 


MACM  NQ,  a 


2,000000 


l/KCRi  ■ ,2000002*01 

2 9701FS 

CMC  1)  9122  a 2 OY  2 


,*00*004)2*01  *,3*000*122*00  •,2400*04)2*01  ,7*217)011-01 

— ,»44)*040**0I J7421 73*tt»flt  — •,  1*4)00402*01  •'-,317200**2*00 

_ CMC  2)  S1ZI  # 2 OY  2 

,1000444)2*01  •,320400122*40  •,2400*4*31*01  ,7*2t7)*U*0t 

— ;r*4)*o*#e*4t — ;7*ji730if*oi — .,  1043*0*02*01  “•,3ira*o*42*oo 


<*  - 


iiaoovMAMtc  C02M!Ciintj  at  ns  ton  tn»oa»  .jtnout  ci«io 

o?ciiL*ro«r  cite 

«*CH  MO,  4 J, 000009 

!/*<«)  • ,I44404(*4« 

- - , iroioi 

cut  I)  sue  liitl 


*,))17}«|}E»4* 


*,172204441*0) 


cm  a *m  • i m i 


-,iM*i*ior*oJ 


MHOtiiitu 

*,1117)4412*0* 


AtMOYt4**ic  r^iuiNd  coemcums  «t  nsto«  ?h*o«y  -ithuut  e*»»t» 
* asctustutY  cost 


*ACH  HQ,  • 0,000400 


• ,2000001*02 

2 ITOtM 


cm  n mi  ■ 2 iy  i 


,iii«7*ru.oi 

~U24*j»ir**j 


*,J«T|«0291.01 

~,«492*549I*00 


M])«7IMt«0) 


, *4520)491 •«« 
*,1**441472*41 


CHt  2)  Silt  ■ 2 »V  2 


*,f 1)47*712*41 

iiUot)i](«(]  •il)2lt)IU«0I 


*44 

>, 11**41*72*41 


A(aoevN4M|C  t«rtut4ct  cotsrjcitNrs  by  piston  tmcomy  *ithout  comic* 


oscuioroav  cast 

«€"  NO,  • 1,444444 

!/»(«)  ■ ,9444441*41 

— j |M{M 

c»<  n im  • i «v  2 


it|*>}ttll«4|  ,121414412*40  ~-,|]<S174IC»41 


• 1 U4t  Ml  t *44 
*,4414444«C*44 


W(  I)  tlZt  • 2 4Y  2 


tam«i44c*4t 

~ii4moot«or 


— i**?^***  ••4»rS42|4«*4l  , 1 2]0 1 0*  | ( *04 

,122414412*44  ••4|421TflC*4t  * >,441444442*04 


4C400YMSNU  7»7lUCNCE  COCFFICICNTS  4Y  PISTON  THIOPY  *ITHOgt  C4N4t« 

o4Ciu.*ro«Y  cor 

— - - "*CN  NO,  « 1,444440 

a ,2404441*41 
~ ~~~  2 *T«lH 

. Cm  u silt  • 2 4Y  2 


_j  11247a 7}2*4I 
,122441412*41 


_>,  147244202*40 *,111470712*41 

i<4UM»M*»I  «,  1)2441412*01 


,•*42*9492 >41 
•,144441471 *44 


....  cm  2)  stu  a 2 iy  2 


111*74 712.41 
,112449411*41 


_>il4TI«Mtt*44 

,444249492*41 


•i1914T471«»41 

•• 1 12441412*41 


,•492*9492 >41 
144441472*44 
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ACR00VNAN1C  INFLUINCE  CQtFFtetlNTS  IT  RISTON  THEORY  atTHQuT  CANSER 

OSCILLATORY  CASE 


MACH  HO.  ■ 


«, 000090 


a ,140000*. 4* 

- 2 STAIM  

CMC  1)  SIZE  a 2 IV  I 

«,I*1*4**7E»0*  i 1*9*0  71  •(  .02 

-,2*000200*. 0* , 109*071*E*02 -,3l**0200C.O*  12S**tA0C*03  ~ 

CM{  3)  SIZE  a 3 »Y  2 

,39140047**04  -,11**1901* .01  -,39140*47*. 04  , 1*9*07 1 «E *03 

— 200*0200*404 — ;i*9*0Tl*C.02 — -,240«0200CO*  -,120**1*0(401  “ 


ackootnahjc  influence  coefficients  or  a: jton  theory  hithout  c<«iii 

OSCILLATORY  CASE 


MACH  HQ,  a 


*,000040 


t/H (R)  ■ ,2040000  *03 

I 

' 2 STRIRS  - — 

CH(  I)  Size  * 2 •»  2 

,1004*l0TE»03  -,24003002*401  - , » 00*«107E .03  ,37 |S t *2*( *00 

— r**3*0*00(»02 ,37|*lo2IC*00  **3*0000(402  -,297«012tr*0l 

CH(  3)  StZl  a J BY  2 

,1004*1*71*01  -, 240*10021*01  •,t004*l*7C«4S  ,17 t* I *2*7 *00 

— ,**3*0*00C.42 ,m*r«3»g*00  --*,M140MQE*02  • ,297**32 1 E*0 1 

AEROOTHAHIC  INFLUENCE  COEFFICIENTS  *7  RISTQH  THEORY  »l THOOT  CAHOER 


OSCILLATORY  CASE 


HACK  HU,  a 


4,000000 


l/KCR)  * ,9000001*01 

2 STRIRS 

eM(  I)  SIZE  a 2 Sy  2 


,42*01147**41  *,*93079042*00  -,4**01 l»7E*0t  ,*2*939711-01 

~«*2100944(*01  ,<2*9197lE-0t  ~ -, *21409001*01  -,*4««9*42E*04 

CH(  2)  SIZE  • 2 SY  2 

.*24*11*7**41  •1*920Mg4(*40  -,4(4*11*7**41  ,*2*4197|C-0| 

,**l**904*4«t  —, *2*9197  It -01  " -,*2100900(401  -,*«**9*02C*40 


AEROOTHAHIC  INFLUENCE  COEFF  ICKhTS  IT  RISTON  THEORY  -IThOUT  CAHSIR 

OSCILLATORY  CASE 


HACH  NO,  ■ 


*,004004 


t/R(R)  a ,2044401 *01 
' 2 STRIRS 
eH(  1)  stzt  a 3 SY  2 


,100**1*7**01  -,3*00100|C*00  ••l(lO«*l(7t*41  ,17Hi»2*E-0t 

— ,491*4000C*00 ,37tSlo20(-4l  -,**1*0*001*00  '-,297*0131*400 

CN<  2)  SIZE  • 2 SY  2 

,100**1*7*. #1  -, 2*0*1002**40  -, 1 004*1071 *01  ,17|S|*20(-4t 

— ;«4S»*0**(*00  — *,17t*U30t-01 — •, «93*0000(*40  -, 397*01211 *00 


», ouox.ow, *02*5, at *04005  j'.riMHOi 

I.MMIMI  J, T| MTt  »0l •«,*•««  IF. 05*4. 5**401, 04 
I 

5,01*04*  »05.4.k00a]O04*5,01*04O05  i‘ tiootooi 

4. 00000*. 01 

•Ill 


4.00*«tl.04*f ,400*0f ,00*2,00**lC.04 

t. ••}»•* «o* 

— ( - 

I.llMlMM.IUIlMt'MIMtttlt 

“*i  80000  0 00  — ■ - 

• I I I 

I 

I .ilMICtll  >| , IM I 

1,40*00003  i,osMit-ot-i,4«*oor.oi- 

i,4«*ooooi  i ,o5s.}i*oi*t  ,4**oooo  i • 
i.ooooo* •«« 

- • I » ~ t • • - - 

a 

l. o*l*ot.oo  l.gitlilMM ,««]*oooo< 
4 

2,oo*«i(.oo*s,4oooo(*oi*4,oo*aiooo 
■I.mimmi  7,«4tT*oa*-i.«*i*aooo- 

I, 00000001 

• I I Ji 

),)*gi7IOlM,M»{i«gM,lal)rt<«l 

J. ii  t iiooi  a. *7*11001 -1,11711005. 

- i 

i, i**i7ooi*t , 7i*2*oo4*i , ]*o}7*  ,01 
l,ltMU*ai  il«7UIMI<3,3i;]H«tl< 

-4.ooooot.ot 

• ait 

a 

-n  jj*i5*  *o**i,  • t>8oooo*i,  iwiooa- 

j, o»soot»oo 

• a J I 

I 

-•,443171 ,98  1,41010001*0, 4*117000. 

I,i7it7(»gg>«,tiijg(*ai*i,373ia>aa 
• ,443171,0a  | ,23At  01*0 1 *0,4*1 17F *00, 
4,oooa*(,oo 

•13  1 

3, 314711. 08*1, *7j«4t -aw  ,11075  0 00 
UK4**(it«  • ,411*41*04*3 , ]>*««*  ,00. 

1.1147U. 00-1.  «7?8»C-8 1*1.314711, 80 
l,  14*4*1  .»•  «, .54**004-3 ,14*. .7.00' 

■'•aii  

a 

2. 51*05005*1, 10*15004*4. 51*050  05 
>,••••>1*0*  l, tS4«7C*0l*>, ••••>(*«*. 

T,  11*011. 81*i;:8«  111. 04-4, II  *811. 81 
4, •••041 .81  l, 05*07001-4, •••«>(. 01- 
4,00««0|.01 

• a i » 
a 

I ,••*•>( .«>*I,*«« 10(.80* | , 98**4f »84 
4, 41*00*. 01  l,7l01000|.*,*l*4*00l< 

a 

!,•••• at ♦•*•>, *0ll0O00*l, 00**40  04 
*4, 4 !*••{. t I },7t0l«00l**,*l*04001< 

5,0**0*014 

*411 


*,4**1404*. 
• ,)lt0U*0« 
I * 

k,44*l>*.0«. 
*.41001*. 00 
‘1.44***044 

• a i 
a 

"l,*0**4000< 

4|«tMN.(t 

>t.Ntoac*i«< 

0»*J*#0<»01 


'*,!>•  71**0 1*4,  >4*  | >1 .08 
4, >4110* •«>.*, 1100500* 

*,54(75001**, 2*01 4000 
4,4411«I*«>.»,>l«41f.*0 


4,*0*l((*«l*t  ,00»*4O«0 
i.muiHH.iutiNi 

>i,»o«m*«i*i  ,o«*«if  *«• 
J,7HI*I«04.*, 01***001 


7'.«>17«C*0t 

1,17240000 

7;«ai7»(*4t 
1, l 74*8* .08 


t;«is«j(*oi 

3,441101.00 

t‘,*55«100l 
I. 4*130000 


7, •417«l*o> 
l.l74«0t-0l 

7*.  •>!  70C*0> 
1,|74«*r*0t 


>‘.*7*11*  .01 
1,722031.04 

*'.•7*11*. 03 
3. 72203*404 


0. 412**1*03  - 

1.  •••841.00 

V,4!>k«t*tr- 

l.**»04O00 


|'.4S«ttt*01 
I',* 1 00*1-8  I 

t*,410lil*0l 
*.*10450*1  ‘ 


l‘.412**(*0t 
1, ••••»€ *01 

• '412**1.84 
].l«*«4C*0l 


310540ft .03 
i;ao*42i .04 

{',05*071  .of 
i;a*Mat.oa 


j;7i«i«i*ti 

*;»74tl*.00 

j',TI*l«O0l 
«|'l7**M  .00 


«:4**i*a«*a 
,*,**•10**01  , 

4. >41 10* *8 4 

.•'•••to*  *ot 


ll 714|***04 
>».1T401«.03 

I’.nouooi 

•2,174014  *03 


• TN  3003 

• TN  1002 
•TN  3 13 

• 71*  3 14 
•TN  3 34 

• TN  3 23 

• TN  3 24 
•TN  l 42 
•TN  4001 
•TN  2004 

0TN  2 II 

•TN  2 34 
•TN  2 14 

ON  2 21 

•TN  2 >2 
•TN  2 24 

*TN  1001 

•TN  1002 

• TN  1 U 
•TN  1 |4 
•TN  1 I* 
•TN  1 2| 
•TN  1 <> 
•TN  1 24 
•TN  •••! 
•TN  «••> 
•TN  « II 
•TN  • 14 
•TN  • |4 
•TN  • 41 
•TN  « >4 

NTN  • M 

•TN  1001 
•TN  1004 

NTN  I 11 

•TN  I 14 
•TN  J |4 

ON  3 41 

•TN  I 44 
•TN  1 24 

2TN  4001 

•TN  2004 
•TN  4 11 
~ ‘ 7TN  4 32 

•TN  4 12 
•TN  4 21 

NTN  2 42 

•TN  4 42 

•TN  loot 
•TN  1004 

•tn  l ii 

•TN  1 32 
- * NTH  1 34 

•TN  1 21 
•TN  1 42 

•TN  1 *4 

•TN  *001 
•TN  «00 4 

NTN  • l| 

•TN  • 14 
•TN  • |4 

ON  • 21 

•TN  • 44 
•TN  • 42 
•TN  3001 
•TN  1004 
•TN  1 II 
•TN  l 14 
•TN  l 14 

ON  3 4| 

•TN  l 42 
•TN  I 44 
•TN  4001 
•TN  2004 
•TN  4 It 

»TN  4 34 

•TN  2 14 
•Tn  4 21 

• TN  4 » 
ON  2 44 

•TN  loot 

•TN  1004 

’ •»»  1 tl 

•TN  1 14 
•TN  1 14 

•TN  j ai 

•TN  1 44 
•Tn  1 42 
•TN  *003 
•TN  •••> 
•TN  • tl 
•TN  0 14 
•TN  0 12 
•TN  0 21 
' ON  0 24 
•TN  O 44 


Aero-Servo-Elastic  Stability  Analysis 

The  input  data  for  the  aero-servo-elastic  stability  analysis  are 
obtained  by  adding  the  aerodynamic  influence  coefficients  (AIC's)  to 
the  servoelastic  stability  analysis  data  deck  and  by  adding  the 
necessary  aerodynamic  control  cards  and  related  data.  The  complete 
I input  data  deck  is  listed  in  Table  7.  The  AIC's  were  based  on  the 

three  Mach  numbers  M =2.0,  3.0,  and  4.0.  The  stability  analysis 
is  carried  out  at  sea  level  for  which  the  density  is  0.002378 
slugs/cu.  ft.,  and  the  speed  of  sound  is  1100.0  fps,  and  at  the 
three  Mach  numbers  2.5,  3.0,  and  3.5;  the  eigenvectors  are 
requested  for  M = 3.0. 

The  results  of  the  analyses  at  the  three  Mach  numbers  are 
shown  in  Table  8.  Table  8 shows  an  abbreviated  output:  it  shows 
the  first  page  which  echoes  the  aerodynamic  control  data,  a typical 
printout  of  the  AIC's  for  M = 2.0  and  k =0.001,  the  three  sets  of 
stability  results,  and  the  eigenvectors  for  M =3.0.  There  are  some 
differences  from  the  servoelastic  results  in  Table  4.  Both  the  servo- 
elastic  analysis  and  the  aero-servo-elastic  analysis  have  22  roots.  How- 
ever, the  15th  root,  which  was  a complex  conjugate  "zero"  in  the  servo- 
elastic analysis,  now  has  become  the  short  period  root  in  the  presence  of 
the  alrstreanu  In  addition,  an  unstable  root,  Mo.  3,  has  also  appeared. 

A perusal  of  the  3rd  eigenvector  and  structural  deflections  indicates  it 
to  be  some  kind  of  divergence  in  altitude.  This  instability  is  a subject 
for  further  research. 

( 


176 


( 


The  notions  Involved  in  the  eight  oscillatory  and  five  non-oscillatory 
solutions  are  summarized  below. 


I 


Mode  No. 

3 

4 

5 & 6 
7 

15 

16 

17 

18 

19 

20 
21 
22 


Dominant  Motion 
Altitude  divergence 
Actuator 
Accelerometer 
Noise  filter 
Short  period 
First  body  bending 
Rate  gyro 

Compensating  damping  amplifier 

Flipper  bending 

Second  body  bending 

Flipper  torsion 

Third  body  bending 


f: 

13 

h 


i ! 

I \ 

j 

I I 


i 
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Table  7 - Input  Cards  for  Program  MPASES  for 
Aero-Servo-Elastic  Stability  Analysis 

C TITLE  CAROS  __  

TYRICAL  AI«-T<j-AIR*Ml*S?TL?. 

AfiRO-SERVO-ELASTIC  ANALYSIS  USING  RISTON  THEORY 
C CONTROL  CARO  ‘ 

' 2 "l*  "'S'"  f 1 ll  5 3 5 I 50 

C SHIFT  EIGENVALUE 

100.0  

c aero" Constants 

1.0  3.0  3.0  0.05 

AIC  OATA  

maCh  NUMdk'S'S 

2.0  3.0  4.0 

REDUCED  FREQUENCIES  

0.001  O'.OS**  * 0.20‘  0.50 

AERO  OATA 

0.002370  1100.0  

numhEr  "of'  veloCTTTT5~ 

3 

MASS  MATH  IX. 

1ST  ROW 


, v v « 

C 10TH  ROM 
SO.O  0.0 


4.S333333  0,0010#l7  -0,3961806  0.0 

C 12TH  RUM 

~ir.  9 miw — d75 wrrniiai ; 

C 13TH  ROM 

1.1333332  -0.2961006 

T r#l  R-ROM* 

1.2979166 

C CONTROL  SUHFACE  OEFLECTION  MORE 


.0  0. 

0,0  0,0  0.0  0,0 

1,333333333-6,666666667 


,U  V. 

1,333333333-6,666666667 


t. 


F 


C DAMPING  COEFFICIENTS,  g 

0.01 

0,02 

0,03 

0.03 

0,03 

C NODAL  VIBRATION 

FREQUENCIES 

45,0 

125,4 

248.2 

100,0 

115,0 

C VIBRATION  MQOE  SHAPES 

C 1ST 

MOOE 

t.o 

0.411195 

•0,122424 

-0.532567 

-0,756204 

-0,756204 

•0.532467 

•0,122424 

0,411195 

1.0 

0,821620 

1.136661 

0.621620 

1.136881 

C 2N0 

MODE 

1.0 

-0.150871 

•0.910842 

-0.990216 

-0,ai90ai 

0,4190*3 

0.99Q215 

0.910842 

0.150871 

-1.0 

-0.630256 

•1.281954 

•0.630256 

-1.281954 

C 3HO 

MOOE 

-0.803195 

0.645103 

1.0 

0.105986 

-0.947894 

-0.947893 

0 . 1 0 5 9 #5" 

1.0 

0,645103 

-0.803195 

•0,259198 

-1,188204 

-0.259298 

-1,188204 

C 4TH 

MOOE 

0,0 

0,0 

0.0 

0.0 

0,0 

0.0 

0.0 

0,0 

0,0 

0.0 

,3333333333 

,3333333333 

1.0 

1.0 

TT^  5 TFT 

0.0 

0.0 

0.0 

0.0 

0,0 

0,0 

0.0 

0.0 

0.0 

0.0 

-0.2 

1.0 

"-5.2 

1.0 

C RIGID  BOOY  MOOES 

C PLUNGE  MQOE 

1.0 

170 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

C PITCH  MODE 

7.5 

22.5 

37.5 

52.5 

67.5 

82,  S 

97.5 

112.5 

127.5 

142.5 

138.0 

146.0 

"T387o 

146,0 

C control  SURFACE 

OUTPUT 

1.0 

0.016 

0.0 

~C~±  NUMBER  QF""InPUTS  to  CONTHtU,  Adtl 


l 

C CONTROL  surface  INPUTS  from  servos 

5 0 mTo 

C OUTPUT  FROM  SERVO  ELEMENT 

1 2 1.0  0.359  -2  0.659  *5 

Z WUNUER'OF  lUPuTS  TO  SERVO  6'L?m£nt 

1 

C input  FROM  BOOY  PITCH  rate 

o — ^—tTrrrrrrtt 

C DIFFERENTIATION  MATRIX  FOR  RATE  GYRO 

0.0 UO -27.0 27.0 *1.0 0.0 

TT3  O oTo  oTo  ofo  ovo 

0.0  0,0 

C OUTPUT  FROM  SERVO  ELEMENT 

2 I T75  0,001 

C NUMBER  OF  INPUTS  TO  SERVO  ELEMENt 
t 

Z iMpUtS  TO  SERvrt  ELEMENT 

1 0 -1.0 

C OUTPUT  FROM  SERVO  ELEMENT 

3 3 TT5 

C NUMBER  OF  INPUTS  TO  SERVO  CLEMENT 
1 

Z rMMUTs  Tti  serV<S  Element 

2 0 -0,1 

C OUTPUT  FROM  SERVO  ELEMENT 

5 2 T75  57723  -2  0.253  -*» 

C NUMBER  OF  INPUTS  TO  SERVO  ELEMENT 

1 
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INPUTS  TQ  3fcftvn  ELEmEn t 

t -1.0  -0,016 

OUTPUT  FROM  SERVO  ELUENT 

5 0 1,0 

NUMBER  OF  INPUTS  TO  SfRVO  FLUENT 
i 

: IWUrncflEUvo  ”¥CEmE'n  t 

4 0 1,0 

_ 0 ^1,0 

. ou'tpiTt  from  SERVO  ELEMENT 

6 0 1.0 

NUMBER  of  inputs  TO  SERVO  ELEMENT 


1 

C INPUTS  TO  SERVO  ELEMENT 

2 0 -2.0 

*C *00  T P U T~F  R OH~5  EffY0~ ETE’RPTT 

T 0 1.0 

C NUMBER  OF  INPUTS  TO  SFRVO  FLEMENT 


INPUTS  TO  SERVO  ELEMENT 

0 0 1.0 

*11  <5  n 3 

C OUTPUT  FROM  SERVO  ELEMENT 

^ S-jj  Ipj-  0 • 0 r~3  1 *4 

1 * " * 

C INPUTS  TO  SERVO  ELEMENT 


C OUTPUT  FROM  SERVO  ELEMENT 

9 2 1.0  , 0.6366 

1 

C INPUT  FROM  BOOT  ACCELERATION 


[■rrti 


-2  1.013 


INTERPOLATION  MATRIX  FOh  accelerometer 
-1.0  9.0  9,0 


OUTPUT  FROM  SERVO  ELEMENT 


u l .u 

NUMBER  OF  INPUTS  TO  SERVO  ELEMENT 


I k f 1 ft  3 : k*i  tM  3 n ! I :Uf  J 


9 0 1.0 

OUTPUT  FROM  SERVO  ELEMENT 


w i » u 

NUMBER  Or  INPUTS  TO  SERVO  ELEMENT 

1 

T>UTP0r'FH0KT'SEH70"  EtEMENl 

10  0 -0.002170 

AICS 

i sr~K,  rsT"H 

AIC  PARTITION  CONTROL 
6 1 

*icF*RTnTON-cufrnrot 

3 0 

AIC  PARTITION 


I ♦06-T^8RTJ9fT*02-l  ,fll260E*06  6,62394Et01  4,99957F*05  5,8»t9flfc*(»0PTN  1 12 

«,02395t*05  6 , 62 5*)4E *01  3.7a«55Er05-9.S9263F*O2-7,768«9Et0S-8.9?2aoE*0lPTN  1 13 
•l.!0tt24t»05  S ,8fl 19SF+00  6,2892SE*05-8.972O0Et0 t-5, 1850tF*0S-2,56200t*02PTN  1 14 

X * IC'FARTITIQN'CflFTTRoL 

1 1 

C AIC  PARTITION  CONTROL 


c aic  partition 

t,Oi602F*05-2,6004lEA02-5,Ol60?E*OS  3.7108TE*01 

180 


PTN  l 12 


I 


4.9auoig»o5  3.7ioft7g*oi-a.gaaoi£»os«2.5afa2of *02 ptn  i 12 

1.  a:c  partition  control 

2 0 

C AIC  PARTITION 

5. 0t602E*05-2. 600436*02-5. 016026*05  3.71 087E*0t  PTN  1 22 

4,984016*05  3. 710876*01-4. 984016*05-2. 586206*02  PTN  1 22 

X 2NO  K , fS  T M 

C AIC  PARTITION  CONTROL 

6 1 

"C  Al'C  "Pa* rTTtTon  control 

3 0 

C AIC  PARTITION 

57250Sl»6*02-l. 498686*01-7. 25041E*02  1.324796*00  1,999836*02  1 , 1 76aoT^O 10TN  2 12 
1,609586*02  | ,324796*00  t ,497826*02-1 .918536*01-3, 107U0E*02-i .794486+00PTN  2 13 
.4,416986*01  1.17640E-01  2,5 1570F *02- 1 . 794486 *00-2 , 074006*02-5 . 1 2400E  + 0OPTN  2 l« 

X AIC  PaR rTT't On  cwthoc  

11 

C AIC  PARTITION  CONTROL 

z 0 — 

C AIC  PARTITION 

2. 006416*02-5. 200866*00-2, 0064lF*02  7.4217aF-0t  PIN  2 12 

f. 99  Voor*02  "7;'a2n«EH) l-l  ,9~9jboFT02-5 , lT2l0cT57  Ptn  2 12 

C AIC  PARTITION  CONTROL 

2 0 

x atc~  pi  RTrrraR 

2. 006416*02-5, 200866+00*?, 0064lE»02  7*,«2t7«E-01  PTN  2 22 

1.993606*02  7,421746-01-1,993606*02-5.172406*00  PTN  2 22 

C 3RO*  KnirH 

C AIC  PARTITION  CONTROL 

6 1 

x a re  p a R'rTTrnw~emTTror  ~ 

3 0 

C AIC  PARTITION 

3".28i6lF>Ol=3.Tff6>06*06-«,531$l£*01  3.3119f£-01  i#2«989E*0t  2.9tt099£-d2PTN  3 12 
1.005996*01  3,31 1976-01  9. 361366*00-4. 796316*00-1, 942126*01-4, a8620E-01PTN  3 13 
"2 , 7606 IE >00  2.9U099E-02  l , 5723 1 6 *0 1 -4 . 486206-0 1-1 ,296256*0 1 -1 , 28 1 OOE*OOPTN  3 14 
C A I C~P  A RTTT  rOfi'ClTOT-RoL 


1 1 

C Ate  PARTITION  control 

2 o 1 

c AIC  PARTITION 

1,25401* *01-l,30022E*00-l,2540lF*01  f.85501E-Oi  PTN  3 12 

f,  246006 *01  1785 5 « 3~E-0l-l  .246006*01-1  .293 10E*00  #>TN  3 12 

C AIC  PARTITION  CONTROL 

2 0 

x rre-pATrrmaN 

l, 254016*01-1, 300226*00-1, 254016*01  1.8554JE-01  PTN  3 22 

l,2«600t*0t  1.85543E-01-I ,2a800E*Ol-l.29318E*00  PTN  3 22 

C 4 TH“KJ  rST"N 

C AIC  PARTITION  CONTROL 

8 1 

X AIC  PARTITION  CUNTRoL 

3 0 

C AIC  PARTITION 

“57257586*00-1. 49888e*oo-7, "2504 iE*00  1.32479F-01  1,999836+dO  1 , IT^oE-OUPTn  FT? 

1,609586*00  1.32479F-01  1, 4978 2E *00- 1 ,918536*00-3, 10 7406*00-1  , 7944 8E-01PTN  a 13 
•4,416986-01  1, 176406-02  2,515706*00-1. 79448E-0 1 -2, 07400E+00-5 . 124006-0 1 PTN  a 14 

x a i c~p~Airrrnroa~cuRTRtiL 1 

1 i 

C AIC  partition  control 

2 ~ o 

C AIC  PARTITION 

2, 006416*00-5, 200866-01-2,006416*00  7‘.42174E-02  PTN  4 12 

i7993606*0<TT, 42 17VE-02- 1,993606 *00-5, 1?2«06-01  PTN  4 12 

C AIC  PARTITION  CONTROL 

2 0 


I 

i 

J 

j 

1 


I 


i 


V..  JL 


m 


181 


d 


M 


1 


c aic  part i t ion  

2,0B64lE *00 -5.200 86E-01-2,006aiE*00  7,42t74E-02 
i.993*0£*0 0 7,'»217«£-02-1.'»<»360E*00-5.172«OE-01 
C 1ST  K , 2ND  h 

C AiC  partition  control 

6 1 


TTfT 

PTN 


1 


“Ait  "Partition  TOntroL" 
1 

A 1C  PARTITION  control 


Jf—  Q 

C AIC  PARTITION 

3.34937E*05-t ,73624E*02-3.3a937F*0S  2.47633E*0l 
“3V31735EV0S  P . a'7 S 33EV0  f^fTSTTTK  ♦ 6 b - 1 , 1 2 2 0 1 i ♦ OT 
C AIC  PARTITION  CONTROL 

2 0 


PTN 

■>TN 


3.34937E*05-t.7362aF*02-3.3«93?F*05  ?.47633E*0t 
3,31 735E  *05  ? . 47633E *0 1 -3 . 3 1 735E ♦OS- 1 . 7220  t E*02 
TT — — 


2N0*‘X7~ 2NCT 

AIC  PARTITION  CONTROL 
6 1 


0 

AIC  PARTITION 

I"  50  Tl  RETD  2 -r.D  0 ! 1 1 E ♦ 51  - a , fl  lf2  7 TfTO? 


1 , 0 7 1 50E  *02  «,R607aE-0! 
-2.93549E*0t  7.70863E-02 

~C ATC- PAR  T TT  TON 

1 1 

C AIC  PARTITION  CONTROL 


8 , 8A074E  «0  r 
9.98866E *01-1 ,2785afc*0l« 
t .A7363E *02- 1.193 12E*00> 
cuwnitji — 


1 

2,07037E*02« 

l,3800flE*02- 


2 g 

c aic  partition 

l,33975£*02-3.472«8E»00-l , 
“1 . 3 2 1>  9 a FVO  ? “079  S266'E^n -T7 
c AIC  PARTITION  control 

2 0 

•c atc-pvrt  rrruN 


3397SE+02  4,95 2 bfeE-01 
269<JE*02-3. 44^02  EtoO 


t,3397St*02-3,a7248E*00-l. 
1 , 3269OE*02  U.9S266E-01-1 , 

X JRO  K# 2N0~H 

C AIC  PARTITION  CONTROL 

6 1 

X” 


33975£t02  4.95268E-01 
32694E*02-3.44402EtOO 


TrC~P’AHT  I ITON-CUR THOIT 


3 0 

C AIC  PARTITION 

~7.*1  R 1 99ETD  l-2^50275ET<n^T7 


6,69<,88E*0Q  2,215196-01  8, 
-1  , 83O68E+00  1 ,927 16E-02  l, 

-c AnrFAHTiTnur-coNTRiTu 

1 1 


•B  28  7 T C*f\  ~2T2T5 1 '9E~=TT 
?429iE»00-3. |963aE*00 
04602F.*01-2,982BO£-01 


o.ia/i  if-.  ♦00 

■t,  29398E *0l« 
•6t6255l£»00« 


AIC  PARTITION  CONTROL 


2 9 

C AIC  partition 

8,37342E*00-8,68l20E-01-8. 
~<r,2R337EA<7<r'r;23f  I'AE-0 1-8 , 
c aic  partition  control 

2 0 

x— 


37342F+00  1.23818E-01 

’2R  3T7  E TirffVJTSTmFXyr 


-ITfPARrTTrON- 


8,37342E»00-8,98120E-01-«, 
8,2R337E*00  1 .23816E-01-8, 


37342r*00  1.238UE-01 
29337F *00-8, 6100 5F -01 


182 


IT 

22 


*c rrcrpiPTHTOfr-cwTRoT; 

3 0 

C AIC  PARTITION 

“BT76798E  »05-5T0055"TETo’2- 1 ,21  068E*06  4 . UTO'fTE  *0  i J.33<J85E*05  3,85«  3 1 E+OOHn  1 12 
2.67875E+05  4,43037E*0t  2 , 497 1 6E*05-6 , 39268E *02-5 , 1 7592F *05-5 , 96559E  + 0 1PTN  l 13 
-7,33873E*0«  3.8S431E+00  4. 184Q8C »0S-5 ,96SSQfc ♦01-3.45020F*OS-1 , 70S?7£»02PTN  1 la 


12 


12 


PTN  1 22 
PTN  l 22 


7TT5B63t-o2PTN"2"l2 

1 , 1 93 1 2E*00PTN  2 13 
3,ai0S3E»00PTN  2 14 


PTN  2 12 

■PT7TZ17 


PTN  2 22 
PTN  2 22 


1.92716E-02PTN-  3"I2 
2,982«or-01PTN  3 13 
8.S2633E-01PTN  3 14 


PTN  3 12 

■PT>TT~17 


PTN  3 22 
PTN  3 22 


M 


C RTH  K , 2N0  M 

X Tl C 'PA'RTYTfON  CUNtROL 

6 1 

C 4 IC  PARTITION  CONTROL 

3 0 

C A I C PARTITION 

~TT5$7i9E»00-l70'5l  1 l£700*4,8a273E  + 00  8.86074E-02  l,33$5at*00  7 , 70863E-03PTN  4 12 
I.071S0EA00  8.86074E-02  9.P«fl6»F-0 l - \ .278S4E+00-2 , 07037F»00« t , 1R3! 2E-0 lPTN  4 13 
«2,R3549F.-0t  7.70863E-03  1 *67363E»00- 1 . 1 R3 1 2E-0 1 M * 3«008E*00-3 . 4 1 0S3E-0 1 PTN  4 14 

*C  A IC  PARTITION  control 

1 1 

C A IC  partition  control 

3 0 

C A I c partition 

1 ,33R75E ♦00*3.47 248E-0 1*1  *33R7SFa00  a'.R5266E-02  PTN  4 12 

~ 173 2 SR  4 FVO  0~a7RT2  6 S E^tf2^T752F9~4ET0  6-3'.  aTPrrerTO FTN'V  12 

C AIC  PARTITION  CONTROL 

2 0 

X STC~PXfrrTTTt37J 

1.33R75E+00-3.47248E-01-1 .33975F+00  4.R5266E-02  PTN  4 22 

1.32694E*00  4.95266E-02-1  , 32694E  ♦00-3 . 44402E-0 1 PTN  4 22 

C rST-KT  3R0~ M" 

C AIC  PARTITION  CONTROL 

6 l 

x atc~part  mafncafirtfnc 

3 0 

C AIC  PARTITION 

~ 5758 8 7 7 ET8 5^37T5T aSETo  2 -97ff^T2  7 FT 3 5 T." J'3  i a i E ♦ 6 1 2 . 5o8*>0f  Ao5  2 .SsbsoFVooP'fN" T " 12 
2. 006146*05  3*33141 £iOt  1 .«7344£*05-a.792SlF*02-3,87<>63t +0b-a.«h293E+01PTN  l 13 
-5,48<»87E*09  2.850S0E+00  3 . 1 3 1 49£*05*a . 46293E +0 1 -2,5828 1 £ *05- 1 .27088E t02P TN  1 la 

X A I C~P  A R T lTnJNTUtfTWoC 

1 1 

C AIC  PARTITION  CUNTRnL 


C AIC  PARTITION 

2.51605£*05*1.30415E>02*2.5160SEa0S  » '.^Ron »0 1 

C * AIC  PARTITION  CONTROL 
2 0 


PTN  l 12 


2,51605E+05-1,30415E+02-2.51605E^05  1.85R07Ea01  PTN  l 22 

2,48402E«05  1 .8*»R07E*01-2.4A402F*0S-1  *28RR2E»02  PTN  l 22 

C 2'nO  k,  3ro  N 

C AIC  PARTITION  CONTROL 

X nC^PAWt  I T IOnTOnTBol 

3 0 

c aic  partition 

Z7S3551  rr<T2^?7527«Sr;(J3-J76'«RiP*d?''6'.6h282E~6"l"l,b!i3^0f^2  5,7'oo'99E-o2PTn  2 12 
8,02a58E*01  6.66282E-01  7 ,«R3R5F >0 1 -R ,58S06E^00-l ,55 1 85E *02-8 . «2587E-0 l PTN  2 u 
-2,lR475f*0l  5.700R9E-02  1 .25260E+02-8.R2587E-01-1 ,03312fct02-2.55375E+00PTN  2 la 

x rrtrp  ATnTTTO^cajrr  roc 

i i 

C AIC  PARTITION  CONTROL  


C AIC  PARTITION 

1 .006426 ♦02*2. 608 30E ♦00-1*0064 ?F *02  3.7l8l4E*01 
-R7RTso€r*oT7rrrmxr-(n^7RW8ETaT^?75TRffinw 
C AIC  PARTITION  CONTROL 

2 0 


PTN  2 12 
PIN  2 12 


1, 00442E ♦02-2* 608 30E *00-1 .00642F+02  3.71M4F-01 
R,R360BEa01  3.7l814E-0t-R,R3608Et0l-2.S7R83C*00 

X 3RD  K,  3R1TR 

C AIC  PARTITION  CONTROL 

* 1 


PTN  2 22 
PTN  2 22 


c iic  partition  control 

3 o 

C AIC  PARTITION 

I,6a7l9£*01-1. 880726*00-2,274326*01  1.665716-01  6,271266*00  1 .425256-02PTN  3 12 
5.015366*00  i. 665716-01  4. 681726*00-2. 396276*00-9,699086 *00-2, 231O7E-01PTN  3 13 
•l , 371726*00  1,425256-02  7, 828736*00-2. 23UTE-01-6,a5?0iE*00-6,38a3B£-0lPTN  3 la 
"C  n C P a If  T I TTOff~C  ON  T RoC 

1 1 

C Ate  PARTITION  CONTROL 

2 o : 

C AIC  PARTITION 

_ 6. 290126*00-6, 520756-01-6. 290126  *00  9.295366-02  PTN  3 12 

6, 2 1 0056 tOO  9, 29 $ 3 6E- 02-6721 0056*0 0-6, aa958f^Jl  PTN  3 12 

C AIC  PARTITION  CONTROL 

2 0 


6. 290 126 ♦00-8,520756-0 1-6, 290 126 *00  9.295366-02 
6,210056*00  9. 295366-02-6, 2 1005F ♦00-6.409586 -Ot 


PTN  3 22 
PTN  3 22 


~C~ * a'TH'K,  JR1TM  

C AIC  PARTITION  CONTROL 

6 1 

*c nc~PrAW!Tro7n:oNTRor 

3 0 

C AIC  PARTITION  . 

276355 1 ETO 0-?; 522*66-0 1 -3, 6389 16  *00  6. 662826-92  1, 00  JaoCtOO  5,760996-0  JPTN  « i2 
8.024S8P-01  6.662826-02  7 , 493956-0 1-9.585066-01-1,55 1 856 *00-8 ,9258 76*0 2PTN  a 13 
-2. 190756-01  5,700996-03  1 .252606 *00-8. 925876-02-1 . 033126 *00-2 , 553756 -0 lPTN  a 14 

C A I C~  PAR  TTTTO*rtOtfTRHC “ 

1 1 

C AIC  PARTITION  CONTROL 


C AIC  PARTITION 

1 ,006426 *00-2.608306-0 1-1  ,006a26*00  3'. 7 18 1 46 -02 

“97936  08r- ori77  in  r4eift2-TJ.9360«ETOTV275798'3r=n 

C AIC  PARTITION  CONTROL 

2 0 

*c rrmRTTTTuw 

l, 006426 *00-2, 608106-0 1-1, 006426*00  3,718146-02 

9.936086-01  3,718146-02-9,936086-01-2,579836-01 

*2  MACH  nunbEMS  PON  vEulClTies  ANO  Eiv»6NV6CHJH  REUUtSl 

2.5  -3.0  3.b 


PTN  a 12 


PTN  4 22 
PTN  4 22 


*fl  - 


»n»um  mlnil 


I,bU  8 ' °^t£ss^ss^r *■ 


i 


i 


I 


l 


i 


! 


i 


•s 


185 


fc r . ^ 


i 

r 


I 


«|C  MTIII  foa  MACK  MIlMCca  ■ }(M»M«|<  AM  AfOUCCD  fftfMC «CV  ■ I. ••••••(-•! 


t 


ItMlUIV  MMLVSIS  KM.lt 


• • » I I I * I * I I I I I I I I I I *1  • I 

'tZltXXKKAILtXZX'ZX* 

♦ /i  * * * 

7 * >»>»••  a * • -•  # 


#¥**S8!y2 

zmsis: 

• •••»••• 


? T?  ? ??  T ?1 
**;X**!X£ 
s:sflSt-A 

AjIMQt  /»#•-•• 

»««>!«  a«* 


{««••••  ^2  5 5 <5  * 


#* 828**2 ******** 8**28* 

*4**Aj«**  — — - » * <M  * A 


1 •••••••• 

*++£*#2a 


I !•>•••••• 

— a 44«f 4A 

— • a * a * m m 

4 ^ • Am  m* 


tvi*  f«MN^4i» 


1^4- 

i 

m • -*  -4 


88*2**282* ****82*28*8* 

A m m > a— m.N  *••.»••-•*  444  ft*** 
m m m a m • a «$•••••**<*»£  m«-« 

m»N<e*m*o*#aaaaa.*m.*,#»ap».* 


m m m A m » a «$•••*•**<*»£  ^*«* 
AtN#4n$A|««M4l4^4  4t«'»f 

r v r ' !“  vT,i 

-*m  <n#44haf  •«Nni4i»  a 4 • «(v 

1 , -K---.-.-.-.---.N  Nrs4 


«#44f4f 

••••••• 

• I I I I • I 

a****** 

*4  4*»N4  « 
Nf*  \ A/N- 

4 m • A • 4 m 
*»  m a n «•  k*  c 
. •»  • m * ♦ m m 

4 J*1*<a  a v 

I I 1. 1 I II I 


44*1444  ♦ 
•••••VO 

I I • I I 4 I 

2****** 
(M  • • • • * a 
m >*.  a * **  >i  * 
4 • • 4 « 4 N 

-f  f»  • *»  <*l  4 /> 

A 4 • A a m • 

— A— «4  MO?  — 

• I I li  I 


ava*»a-»a* 
»••••• a • 

iitiiiii 

2*8*8**8 

4tNC44fl» 

4»^4l4«t 

■ m **»«*  + *> 

35S32S*? 
•«««•««« 
^ a*  m ai  a «»  m m 
i.  (<ii 


i i i I • i ( ( 

8****8?* 
n a « m <*  c a « 
m \ a ^ *'  a 

• a 4 **  a a a • 

• • m # **  m •*  f 
m m a >*  «•  m * > 
aomMaa^A 
«»*!••»• 

III  III 


188 


1 


1, 


o 


! J J 


I 1 I 

:::ssi£ 

i • i • i ti 


* I ; I I 


ztzstzz 

• ••III! 


<09^00 

> • » 9 • *• 

► » 0 ♦ « ^ 


2 ! 


• i • 1 1 « • 

# # — * * 
> • * * * - 

• ♦ - - ~ •» 

T • i 

\ 
i 


» • l 

2 i • ^ • 2 Z 

I 1. 1 I I I • I 

s&ISSISf 

S51-  SS?  N ? 

1 I ‘»  1 I I * 


I f I 

sxsstsss 

t I'  I I I I I I 

$2**1$;:* 

$ iE 

- - - , -Jl# 

* ^ ♦ 

^ « 41  il  *> 

» JM  1+ 

I 


t 

Ite. 

I 

s*  -f  j3iJ 

4 «■**  <n  I 

v f*  a*  » *i  r - 

4 • • • • 

(*i**0  — « 


« I 


I • ' 


{iiiiiti 

iH2*3§2 

I?  3 • 3 a a • 

|SS»S£|| 

‘41  — n « aim  « 

• 1 ,i» 

I 

i 


! • 1 • 1 I • I 

xztxxzz 

#ji«na«i4 

9 0 • .n  <*  a 4 

* » « iM  » N 

:•  ■»  • — **  <+  ^ 
• * -•  • 41  41  • 4 
• •«•••* 
O m *•  — 0 m 0 
l l l 


I I 


1 1 

issss^ssi 

; 1 t 1 • 1 • 1 1 

!**X*$**!X 

*09  9 9 9+1+9 

"*isc*"*  ♦ ♦ * 

! * ± 9 9 5 f 2 « 
.uni****** 

♦ il-'*  0 0 9 9 


)*>  99  — 0001+ 

■ 1 1 • » 1 • 1 • 


imim  * 

0 9 — 00  m ±0  • u 

i « 

$3  5 

-IfJ  ! — 


*>  !**+  — 0 * 0 0 
9999909+ 
w(Vn  9 9m 


— — 4I«.**<M- 

r r 


•»••§•••* 
^ yi  ^ m 

• + f *•  tst  m 1+ 

liiSSSS 

jirssSa 

9 0 <C<*1  — 0 — 
• % * 


9 9 — A O 9 •+ 


I I I I I I I 

990099+ 
**  + 9 + O'*  0 
\ m**.« > 0 9 m 
nn-f*  • il  *n 
0 *>  9 9 9 0 1+ 

0 — A 0 0 9- 
• •••••* 


1 1 • r 1 1 i.  1 


ssaus*** 

-00—00—4 

*>**99  0 99*9 
9 9 00  0+*+  M 


Oj  - - — (N  <>  - ,v 

• 1 i • 


I 

SmmmmmZm 

1 • 1 i«  1 it 

I'M**  JUXX 


1 


t 


,«•»  41  4*41 

f I I I I I I I 

1XK2X??* 

,9  0 0 09  4** 

i*  1+  0 + 0 r “* 

ifl  **m< 


\*t  ■ 


IAN 


«•  (M  4 M»NOT^ 

I t I I I 

0+00+00 
\0  mm  0 0 m 0 

1 1 1 t • • t 

<2S!5  X XS  Sf 

10—0  1+0  0 — 

'S?SS52S 

!«•  # Ai  « t*  ** 

l*m*~0 m 

t+Otmrnm  0 0 — 

i 'j 

t»*  i+0+lO 

'mm  m mm  m 


• t 


: + +4 

0 1+  0 0 0 0 

•*  M 0 <+  0*  0 

.'*•  »M  9 —**  — <* 
• II  I t I,  I 

\m  90  m m 0 0 

• 1 • 1 1 • 1 

X mS  « 2 £ X 

• j • — T «♦  jj 

iiNaoi 
A n»4N^  A 

;«••  • + > 0 <+ 

\0  + -*A#  0 0 — 
t I I I 


4 « 0 9 9 0 
, 00009000 

I t I I I I • I I 

!***!***« 
•5  2#I3J{|? 

.9^909909 
\9  0 0 99  9 O — *+ 


O 1 
• • • • 


! 1 1 1 1 • 1 t 1 

\X* 

'**9090990 
x *+m  0 0 0 m 4 
m 9 0 + — • 0 

9 9 m >•  0 * m <+ 

***->9*0 

NNNN4lf^4*» 


0 01+9  + 


•■*J 

:S£33 


I I I I I I I I I 

|!XXsyS!X!K.tl 

sSsus:£3 


\nm 


• a 

_ I<  A*  N 

•<  % ill 

*JS5«T^ 


I I 


I I I I I 


'NlSgSXXS! 

• 41^  ♦ AN 

1 Nl»*.  « > 

121*2 
I A^ 


»♦  a - — * - 

Al  XNX  ■ 

*322W«. 


• • I 


;•  Ai- 


>9  0 0 *• 
it 


itzzzit 

I II  t I I tt 
imih 

222 32 3 -i I 

4 + m m 

Unix 


4*  — *•  # — m 

I 1 •<  1 J 

222222s 

i • l it  t T I* 

si  sx*— ssjf 

—2  #*22  • 4 

■•*4  4 * * 4 4 

alsxs2s 

4 4*  w*44  4 

r?l 

— *»  f*  — — IN 

44  4 • 4 4 

II  I I I l> 

3;  *123 

2$  5*52 
Si  2mU\ 

m*4+  4 — — « 
ll  ii  , 

5555555 

♦ •••III 

sxaxsss 

N 

/!•••>•• 

;i2£4£S 

23222SS 

imaI^-ZaJa!  a 

• • .ii. 


! ::::ss2 

» « * 1 1 * 1 1 

: sxxxatsx 

• 4 4^  #44  — 

4 4**  A — « 4 
A*  4 * f 4-4 
4 *.  — **  — <•*  — 

; U5XSSS 

4 — # — — — — 

I «,«  | 

♦ A 4 4 **  *»  4 

m mm mrnmm 

iiiii*i 


• * ,ii. 

i ? ! 

55555555 

I I 1 I I I I I 

*1151111 


•I--#  - 4*4 

•;  1 ; 1 j 

55555555 

I*  I I I I I I*  I 

S|5!?J?S 

fit 

Zzztstzz 

M >l  I - <M  - > 


unm 


55  5555 

♦ I III*, 

**  x ; 
«5  5525  , 

4 • . A - <M 

\ * « in « ; 

■»  A A . A # • 

«*  > »»»*.•  S 

AN«»«*lJ  ; 

T .,  J 

5555555  ! 

■ i • i > i . 

i*nm 

m A m # A « •*  ; 
— 4 — 4 41  • 

♦ iittni«  • 

♦ 4 •/»  4>  ♦ 4 : 

♦ n 4 ♦ A 4 — | 

A 4 — — « — 4 * 
I I II  I 


55555555 

• I I I I I I I 


— — — ^ — 44*  4 

i rl 

55555555 

I'  I 1 I I • l»  I 

mnm 


4 — #*  4 4 — — m 
II  I;  I I j 

55555555 

I • • I I I I Ij 

H £?**«* 

iiisiljs 

« 

i II  i i » | 


f 

! 55555555 

• 11*11  I*  I 

! *35f ]SS! 

a 55555555 

I SJiJiiil 

3 ♦♦A#NN 

* •(  • , 


55555555 

. • » . i.  i . i 

£*2.3335 

■ 5?i»  8 . 2 5 
5525t5S5 
12555555 

“7T 

5555555. 

. I* 


l J 

t*  # 4 ««  4 4 

• • m mm  m m 

; l l • •'  I l I 

<n  4 > — m 4 — 

4 — 4 *»r  4 ♦ 
I ? t Mt  4 ♦ 

5555555 

• 4*  — m 

J 9 A #C  .*3  — a* 

i • * i • 

555*555 

iiisni 

$ J 2 2 J 32 

— — 4 a 4 4** 

;il 

; — 4 — *•  — 4 — 


fkfH  ♦♦»!••• 

• 4 ••• • 

II  I I I I 


r*  tssx 

— n *.  ,v  ♦ 4 
4 4 t!S|tt 

5 a ♦♦♦ t 

4 ^ m ♦ + 4 
4 4 V A 4 • 
♦ > • > « A| 

^ <*  • r-  — m 
, ••• 

22242  2*2 

i i • T l l l 

3H2213 


22222222 
• i * i i i i i 


<*  4l 

inzstzs 

4-4445144 

41  .«  • a a « 

~ ft# 

4 4 

— — 

4 4 

*^4 

— 4 — — 1 ~ # 

8555S551 

— 4 

— — *n  • a m — 4 4 

r r j i 


,222  222  2 2 

: l I I I I 1 t I 

*¥*3$*s** 

— ••/y  m "• 
^ N -*!♦•  « • 

» t/14 

: • ♦ « <n—  -•  > ii 
* ♦ > <n -•  ai  <m 
l • • • i • 


•••j  ■ 

555S5555 

355tjt3(*3« 

• 55555 555  . 

I n-r**2  * 

! 5555555.  | 


42222 2 2 

* i • • • i ii 

w *iwi 

<M  ^ ♦ ♦ « « * 
*»  » ^ • ♦ IV  *■ 

S224221 

— 

#2 

•*  • i | 

♦ 4 •*•*«» 

••♦•• « * 

I I I I I I I 

i3ft*SK¥ 

iSiisn 

tit 

. * • . 


f ^ ♦>  m m 
mm  mrnmZ 

ii  i • i i 

4 t S • 25 

\ ^ tm  > # 

MAIf 
*»  « 4*A|  A ^ 

55  5555 

**•1111 

*3««4I** 
» I *0i44* 

3255555 

A|K«mH4^ 

ii  I.  I « i: 


i mommmmm 

• mm— — —— • 

i • • i i-  • i 

: xsyys*:*  ; 

I • # P*  rm  — m mm  j 

« + m 41  !■%  * — 

I .m  - m 4*  4 Ai 

• 4*4  » — •- 

1 <C  #*  - 

I I 

. . — .5-5 

j • «■  I I V • I 

I 

, ••  f *• 

5855551 

• !;555;a  i 

♦ t ^ *i  » 4 • 

<**♦  /i-»^  # 

J 1*1  H*  l 


S 8514*3*5'*!  S 

W (J 

Ii*  4"*f^*»t  A 4 
» ♦>*.**♦••♦  » 
X C 

5 5 

-•  •A»i*4t-i*l«  -• 

4 I III  4 


i i i i:  • i i • 


X !9 

- 4 *-41 ,5 

4*  . I I I I '«i 

I I ■ I 


2 2 2222  3 3 

• • * * • * i i 

*22*23t* 

.*  • A m a m * * 

♦ *>  ♦ 4 .* 

S2S22222 

> | 4 

4 4"»«-/l-  ♦ 

4/IN  2—  — # 4 

• • j II 

4 4 41  m 4 4 4- 

• 4 4 ——44  4 

I I I I f • II 

*^*Ss!*it* 

"*•»  />/»!•  44 
f 

A**-44»* 
> /I  ♦ «V  9 •»  * — 
-*/***•%• 

— 4I4*  — 

•t  • | • 4 

i ! i 

55SS.55. 

I I I I I I I • 

X#*S11!.!S88 

S5S555SJ 

^5J*S8«3 

♦ — — Al  A.  ♦ 

A*  • » n — « a 

• ••••••• 

4Ai<4*»4  4 41  — 

• 1,1  I 

55555558 

I I I I I I I I 

JHHS5S 

A — AI4  — 4 4 4 
— 4^  — A — > — 

SliiUt: 

4 f *-♦ ♦ 44 
III  , I • 


! 55  5555  ! 

II*  I l I l 

: mi  SsfS  ! 

— 4 »»/*• 

24  232? 

44  4 »»  4 N 

A <*  4 4 — 4 | 

4 — 4 4*  A#*  — ; 

t 1 \ * 

4444444 
4444—44  • 

I • • I I I I 

2X3X 5?S*  ; 

4 — 4 — 4 4 "•  * 

4 A m > 4 4 - , 
4 4 4 — • * 

. • |l  4 • * • • 1 

• 4 4 4 4 4 <*  4 ; 

m A m ••  4 > a . 

(4  «.  — <v  4 *•  — ? 

. I*  I • t 

i i ! i 

! 55555555 

• till  III' 

■ A ♦ 

^ — « > — »«.  M 
4*14444## 

• ->•*  •*/> 

i 4 # 4 — H 
’ *»4  il«444^ 

| — il4«^  — — 44 

• I Ij  I I ’ 

! 55555855 

* I I * I I I * Ij 

: tut***** 

+ +*+*»e%  + Am 
4—44^4—4 
**3  444#4^' 

-•  4 > — 4*  » - 4 
I 4 — 44#4/l4l 

I 8^#4*<«4^ 

• ill  I)  I . 


3 22  2 2 » 2S 

ii*ii*i* 

SJSJ52;: 

4 ,v  4 4 — 4 in  a 

?S 222242 

4 — 4 4 4***  4 4 


?#<  ^ 4 — — 4 4 

4-44** *4 

4 4»—  — *4  4 — 

7 (1  j 8 : 

i 22222222 

• • • i • i • • 

t n«sji? 


m • i .1 


1 


190 


ztxtiit 

1 a i < • • i,i 

ismwi 

illllii 


1 1 J 

X»*2XXStf 

3332222: 

sh=S£® 


TT?T?*T 

Hp 

iii* "v  I 


!S3  i:i: 

ii  iiii 


•|  • V • I 

fgfjfg 

lisflif 

n i 

IS  £238 

II  IIII* 

XX  «**g 


f iiii  i ii  ini 


r--r??r, 

3333333  i £233233 

• ••IIII  I 11*1111 

inini  »ih» 

itxUH  i 5|32ais 

|333l23X  : 323253S 

!’»!!• 

t*mnt\  mum  t 

12333333  i ? 

iumn  nrAni  * 

33322333  33332333 

I I I I I I • I I*  • I I I I It  I 

tltttnut  ******** 

IMhU  liUuh 

•••*••••  • ••••••• 

«-«?* 

33333322  3332=223 

i i i.i  h i • i i i • • i 1 1 1 

|!**3I*| 

— ♦—•f*  ;qi  «■»!•••*•*♦ 

j j?il:sh?  2 slslllfl 
! 32*33332  5 823=2323 

. hi  It  i i i _ i t i i i i ti 

lip  i|i||!i| 

i SwIHll  I *ls|l|l4 


-•a* 

*4 <"*  — *»-•*  «j*l 

j VT1 

8328=333 

I I I I I I tj  I 

5*5  •*  • $ 


Slf  ! 

til  ! 


§131 ! 

li«; 

• • • % j 


i • * 

823=2323 

• ••••••«  I 


| ' 
5 i 

y I 


3 != 


4=  A=«n«*«T 


I J 1 * 

323  5 ! 

I M * 

*p  * s 
ill  | 3 

?T  s ? 


!::; 

•Ill 


U 


• • • | 

• •• 


= i 
*:*.»  j 


•mi* 

X*2* 


3JS3  !3|23 
22  = 2 i 188.3 


«•  A*  ♦ • 

• • • • / 
ft-*  *A 


, « «N« 

i » m i • 


^ « 4i  «• 

» » • • a 

Mill 

12132 

•8*23 


* » • » 

2222 
^ ^ 

2222 

♦ 8-  - 
• • * « 
4 * 

1 I 


• • if 

■•  i* 

— :s 


4 • • J 

mm  m i 


f 4 4-4 

9 * m m 

IIII, 

***« 

a + m -i 
n > • o 
m # n # 
* * 4 ^ 
• • • t 

i*|  M f*.  >. 

1 • I 


| ! i i l| 

• • • o I * • • • Q i « « • '5 

*«*  I >••••  I • t 

• * * i » 


a a 


|i«x  | 

3 8 a 2 • S 

•ii 3 • * w 

• « a ••  ^ 


III*  ^ ,111* 

| **** 
»*»  + • — • • — a 

a ^ • u « a a — a 

I » ^ 4»f  » 

^ ^ M a a a 


22  f |223  | | jf 2 2 

^ a S S J n a ••  a *•  V ’•* 

♦ • ••  4 If  4 — — **  /»  M 

*4  U ^ • W ^4  If 


r*r 


• i • » ^ • * % » 

•li^N  - M -A  A*  <3  4NA« 

a »♦  S . i i 


f W ^ 4 ff 

* a — a m •*•  a a n 

>N  « • ,j  f 4-4 

• I • • 4 • • • • 

♦*»  m 4 

I I ,3  IIII* 


191 


• Ill 

I 

> M — » 

1*1  A M I 

« >•  ! 


' ! i 

M ai  M — i 

a a o a ; 

• ill  i 

****  j 
sC2S3f 
2S23  ! 

— A AM  : 

««-• 

• • • • I 

^♦4-  I 

• I.  I I I 


I I.  I I 

j 

jnw- 
mm  a m s 
m 4 — — I 
m — m M I 
/»m  j»  — i 
4 — I 
• • • • 


I it  I i • 

**a* 

4«4N  1 

mm  — m >4 


— — a » 

m * a ! 

♦ MV  . 
A — — I 

• 4 A 
• • • 


»*  ■ i 

«]  M M *i 

• III 

• III 
4l4l4l4l 

• • vS; 


fiiss 

A»>  ♦ 
• • • » 


-•MAM 

I'  1 

M—  V — 

• A A A 

• • • I IJ 

**** 
ONOA 
<■»♦«* 
NrtNN 
AM  4<M 
» *•  MM» 
4*4  /I  Af 


M fINN  I — «(W»  — 

^ »•♦««  <W  UJ 


I I 

4 LLLL  ig 

V44  A A A A ;lto 

i ssii  i 


I 1515  I SSsS  I 

— 4 44f»  — A A A A « 

— m mm  * ■»  9f»  go  a 


am  a 

*•-•4  { 

• • • • 
flirtAA 

f 1 

A A A A> 

• I I •. 

ft*** 

AM  A A 


I IX  Alai 

i??5 

;ft*« 

\Z12 

iZSX 

ISS7 


ri 

*$»  d 


! • t it 

:*** 
« > — 
4#a 

S » * 

1«  »•  A 

. • • • 
O 4M 

• ' i 


.sas 

.AM  A 

. • • • 

M — A. 

n 

MMM  — 
A A A A 

• III 

M *2* 

2S2m 

4.nmwi 
•M  AAA 

• • • • 
|-44* 

* I 


I I I.  I 

i M*4t 

nli 

. 4 M A A 
inN«N 

4 M 4 A 

I • ♦ • • 
!4  AO#M 


A AM 

AMM 
' — A A 
•M  4 4 
lAX 


I ?T? 

i US* 

♦ MO 


.*^54  — 

o 

-♦  A A ! 4l 


-♦  A A ! 4* 
A «i 

• • » • Vfe 

—MM  A 9 V 

• I |3 


• MSA 
4»A  A A 
A A A A 


»>  — # — ; 


« ! 1 ? 

0 !<*««• 

*»  1 A A A A 

1 !Hii 

— iOMMA 

— • M A A • 
O • A M A A 
*4  A 4 M A 


M A A 
| — 44 

- 4 M 


’ till 

; ft*** 

***$ 
• Mam 
M A A A 
» m • 4 A 
I AM  4 A 

I 4 MM4 

1 »’ 


• • • I 

I *«** 
; -?ass 
1 2 zzz 

1 a — a » 

a # 1 1 

• • • • 
— M4  — 

I* 

1 

l 

rtfliy 
AAA 
• I*  I 

$*¥ 

* MAM 

A AAA. 

1 233 

M A M 4 
5 — AM 

9 • • • • 

C M 4*  — A 

« V 

O A ♦ «A  A 


4 4 MM 
A A A A 

I « • I: 

ft*** 

4 AMM 

*S)!S 

2ftJ. 

-•AMM 
• • • » 
M A — m 
• •! 


I I I It 

2*»S 

M 4 A A 
A 4 m a 
A 4 MM 
♦ MAM 

2*52 

MM-J 
» » I 


J M A M | 

8 sst  I 

a • • • • 
3 4 A MAI 

3 MMM  A 


I J 

,MMM 
AAA 
•ill* 
* W Ul 
4 M 4 
- M *A 
’ O M — 

■ M A A 
I 4 M A 
; 4 M A 
M MlA 
’ • • • 
;M  MM 
I • 


. M A A 
4 4 M 

A 

- ft  - 
• • • 
. 4 4 A* 

• J 


— AM 
MMM 
M -*  M 

• 2« 

• • • • * 
1 4 — M A 
. I •; 


• I I ll 
4 4 Ul  4 
— AM  A 
•A  — A A 
— 4 M A 
•M  — A A 
, 4 MM  A 
' 4 M A A 
•A  AM  A 
| • • • • 
4-41  — 

I • •! 


MMM  — 

: A A A A 

; 1 • ••  1 


I 4 4-M 

AAA 

1 •;  1 
• *i«2  w 


M 4 O 

3!*? 

A — o 

• * • • I 

— AM  A J 
•: 


A OO  A 

III* 


• 44  A 
0 4 4 A 

. m A — A ! 

. M 4 — A ; 

AM  A A ! 
M 4 O A t 
• ••<*! 
I M — M — I 

i * 1 • 


“*1**  T I 


III 

I 

$rs 

A A 4 

A M A 

>00 

A— M J 

; • 


A A A A 

till! 

!ss; 

M 4 4 A 

ssiS 

M 4 4 A 
• • • • 
— M — 

i—J 

A A A A 

I.  I III 

£«££ 
ft  M M A* 
4 4AM 
MM  M A 
A A 4 M 
M 4 Al  M 
44  > 01 
* • • ♦ 
4|AM  M 
ill 


) • It  I III! 

j***  . **** 

A A A A — — A — 

• M — «A  A — A 4 A A M 

* 4 4 M « M • M ♦ 4 -M 

4— M4  M— MA 

Irtflfll  A 44AM  'w 

>-A-  | O— A—  .5 

1AM—  t AM4M  .1 


I*  ! 1 ? 

lO  — — M — 


» I l 

| fy'K 

0 4AM^ 

— M A 4 O 

— A A M A 
U A 4 A A 
Ui  A • A A 

mJ  M 4/1  A 
««  * • • * 
^ A A — — 

1 1. 


• III—  till 

| *ft**  | **** 

9 m M 4 A o 'AM— j 


* * 
M 4 M — O 

A A A A ,M 

III!  ' 


M A — A 

— A M A 
A 4 4 A 

— A m A 
4 M A A 

• • * * 

— 4 — M 

• • [ 


A A M 4 

<=& 

— 4 A 4 

M 4 A 4 

4>-S 

— 

4 4 4 M 

VJ 

4 — 4 A 

U 

A A — A 

MAMA 

• 4 

A A AM 

m a A 41 

• <4 

• • • • 

• • • • 

4» 

4 — 4 A 

1 

* 

— 4 — M 
! 1 

.3 

• A AH 

>522 

• • • m 
— | — M 

J--J 

A A A A 

• • • tl 

2 3!  5(3!! 

MM  4 M 

A ft  MM 
— 4 O M 
*44« 

• 4 4 A| 
• • • • 

• 14V 
•5  • l 


mm 

9««  o a-*. 
» • i • f • • t 

4 **000  » 0*  ^ » 
0 * i «•  # k \ a O i 

— 4 — <n 

9r4«it  7M.  / • ! 

0000«04-*i« 


>m”1r*390004 
• I » • I I • I 
(U>‘  4l  W li«  W Ul  w u t 

>n»Nte(n«i  « f»iAf 

>0046S«0N<' 

i n n ••  4 «*» 


. i 

I 

I 

I 

— Ai  — — rsjr 


I » i l a 


-«  4 <0  2f 

— **  sa  #• 


I II  I I 
4l  ul  tel  ui  w *• 
> 4 * * o a 
-#o-a. 


; i I 

0 Jl  O ~ ~ — MiJ 

• (♦♦•••♦ 

^Mteiulte 

001*00  a*  a* 

4At9«^OIMC 

9 m-*#-  <v  4 > 4 
4 0 -•  "i  4 * * a 


telteJtelvUteiuitelaJtelW4^lteJ*ii 

a >» 

^ « «#«-*««  rtrt  000 

*1  01  * * tft  I*  0 * <V  4 * * <M  0 

• ♦NinifliMnii 

*-999>»<te»  M*l  4 4 0 0 


• -»  ««  «•  mi  rtf  rtf  (\ 


! SSS'XIiSS'S 

N«t*<1N.NnK| 

1 00  ©«*•<■>©  i% 

2 "•  «*  a an  a 

I «*4  *#  0 9 000 


i aoQociao  © 


! 

«»  ft]  * Al  A 0 0 1* 


******** 
* 0rtlH-0  9 00 
*1  4 4 0«  0 * *1  <"l 
0 4 0 0 4 — 0 N 

0 00004010 


Tim 


i a a 'a  s ns  a % t a a a $ ? a an  » a 

•00^4400000094400 #0 

l-iili5SSS!SSJ5'S#S*2 

• f 4 * 0 4 0 0 0 0 0 *»2  _<n  0 4 0 S 
»«090000099a%09  4 9044 

| 4**-*^-»«0'*-*rt#  0 -te  4 9 <* 

«llt  I l;  l l Ml  M i « n i 

. I • ! ! I I 


* 001010 

\ i 

I I 


193 


MINIMUM  KtM  COMMON  LCNOIM  HC0UIHE0  ■ f«*ft 
MASLD  ON  INPUT  DATA  AND  ANALYSES  REOOtSTtO. 


SYMBOLS 


i 


In  this  modification  to  the  original  program  PASES,  the  general 
notation  has  been  retained.  The  nomenclature  used  in  the  basic  develop- 
ment (Ref.  1)  has  been  reproduced  in  Appendix  A.  The  new  symbols  intro- 
duced in  the  present  modification  are  defined  here.  The  reader  is  referred 
to  Appendix  A for  the  definitions  of  symbols  not  appearing  below: 


A 

*0'*1 

B 

br 

C 

c 

\ 

SiDh 

Si 

Sr 

Sis 

c 

D 

E 

f 

g 

I 

Sj 

Kji,ICkj 

K 


Coefficient  in  spline  fit,  Eq.  (57);  also  see  Appendix  A 
Coefficients  in  spline  fit,  Eqs.  (67)  and  (68) 

Coefficient  in  spline  fit,  Eq.  (57);  also  see  Appendix  A 
Reference  semichord 

Element  of  discrete  damping  matrix;  coefficient  in  spline 
fit,  Eq.  (57) 

Element  of  generalized  damping  matrix 
Element  of  oscillatory  AIC  matrix 
Element  of  damping  AIC  matrix 
Imaginary  part  of  C^ 

Real  part  of  C^ 

Element  of  static  AIC  matrix 
Reference  chord 

Coefficient  in  spline  fit,  Eqs.  (57)  and  (60) 

Modulus  of  elasticity 

Cyclic  frequency 

Aeroelastic  damping  coefficient 

Element  of  unit  matrix;  spline  bending  moment  of  inertia 
Element  of  spline  interpolation  matrix,  Eqs.  (74)  and  (75) 
Spline  parameters,  Eqs.  (70)  and  (73) 

Element  of  generalized  stiffness  matrix 


1 

1 


i 

i 


I 

j 


N Element  of  generalized  mass  matrix 

P Spline  reaction  force 

Q Spline  loading  parameter,  Q-P/12EI 

q Dynamic  pressure 

S Reference  area 

s Reference  span;  Laplace  transform  variable 

1*1/2  Time  to  half  amplitude  (negative  value  denotes  time  to 
double  amplitude) 

U Element  of  matrix  in  canonical  form  for  eigenvector 

analysis,  Eqs.  (45)  and  (47) 

w Spline  loading,  Eq.  (56a) 

x Coordinate  along  spline;  also  see  Appendix  A 

y Spline  deflection;  also  see  Appendix  A 

A Eigenvalue  (*0)  in  eigenvector  analysis,  Eqs.  (45)  and  (46) 

U Decay  rate  coefficient 

p Atmospheric  density 

(")  Denotes  complex  amplitude;  denotes  inclusion  of  aerodynamic 

damping  and  stiffness  in  C and  k,  respectively. 

Subscripts : 

r Oscillatory  mode  number  in  frequency  lining-up  process, 

Eq.  (53) 


s 


Number  of  mode  being  lined-up,  Eqs.  (53)  and  (55) 
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